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PREFACE 


The topics discussed in this fourth Aero-Astrodynamics 
Research Review cover a variety of subjects. Included 
are Fluid Mechanics, Aerothermodynamics, Structural 
Dynamics, Instrumentation, Selenography, Space En- 
vironment ,^and Orbit Theoiy and Prediction. It should 
be noted that no single review treats all subject areas 
being investigated by the Laboratory. Further, these 
reviews are not progress reports; rather it is attempt- 
ed here to publish works of a research nature which 
have reached a milestone worthy of note. Other sub- 
jects will be discussed in forthcoming reviews. 

It is hoped that these reviews will be interesting and 
helpful to other organizations engaged in space flight 
research and related efforts. Criticisms of this re- 
view and discussions concerning individual papers 
with respective authors are invited. 

E. D. Geissler 

Director, Aero-Astrodynamics Laboratory 
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I. AEROTHERMODYNAMICS 



CONVECTIVE HEAT TRANSFER IN TURBULENT, SUPERSONIC BASE FLOW 


By 

ifc ^ 

F. R. Krause, W. K. Dahm, R. E. Larson , and A. R. Hanson 


SUMMARY 




Total hot-wire bridge voltage, + e^ 


Turbulent, *convective base heating has been 
studied in a Mach 3 wind tunnel utilizing a two- 
dimensional model with a blunt trailing edge. Meas- 
urements of heat transfer and base pressure, along 
with probings of total pressure, total temperature, 
mean and fluctuating velocities, and flow directions 
were performed inside the recirculation area behind 
the base. The results show that the base approaching 
recirculated flow establishes a ,T base boundary layer" 
which provides 86 percent of the heat flow resistance 
and therefore dominates convective base heating. 
Furthermore, the plumes between clustered rocket 
exhausts produce extremely high fluctuations through 
impingement shock and free shear layer interaction. 
These fluctuations are then convected or radiated 
backwards towards the vehicle subjecting the base 
boundary layer to turbulence levels exceeding 100 
percent. These levels suggest that worse problems 
may be on the upper stages and that turbulence 
parameters may have to be considered when scaling 
base heating results to full scale conditions. 

LIST OF SYMBOLS 

Symbol Definition 

A Amplitude parameter or ordinate inter- 

cept in King's law 

a Local sound speed 

c Specific heat 

P 

dA Differential area 

dV Differential volume 

E, Mean hot-wire bridge voltage 

b 

’ Fluctuating hot-wire bridge voltage 

1 /2 

e_. RMS of fluctuating voltage , ( e, t2 ) 

RMS b 


f Frequency 

H One -half model base height 

H(f) Transfer function 

h Heat transfer coefficient or probe face 

height 

k Thermal conductivity 

M Mach number, U/a 

p Pressure 

Q Total heat flow rate 

q Heat flow rate per unit area 

Re Reynolds number 

St Stanton number , h/c p U 

P 

S Model length 

T Temperature 

AT Temperature difference; e.g. , AT^ w 

= T - T 
w e 

t Time 

U Mean total velocity 

U* Fluctuating component of total velocity 

IT ^ RMS of fluctuating velocity , U 

RMS KMb 



, ? 1/2 

U Total velocity, U = U+U ! -(u z + v ) 

u, v Mean velocity components in x, y direc- 

tions 


* Litton Systems, Inc. , Applied Science Division, St. Paul, Minnesota 
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LEST OF SYMBOLS (Concluded) 


l. INTRODUCTION 


Symbol 

Definition 

u T , v* 

Fluctuating velocity components 

u, V 

Total velocity components, u = u + u f , 
V = V + v f 

x, y, z 

Coordinates 

Jt 

Boundary layer or shear layer thickness 

P 

Dynamic viscosity 

V 

Kinematic viscosity 

P 

Density 

a 

Spreading factor 

CO 

Specific weight 

SUBSCRIPTS 

A 

Area -averaged value 

B 

Recirculation zone value 

b 

Base value (usually at center of base 
plate) 

e 

Equilibrium (recovery) temperature of 
hot-wire 

H.W. 

Hot-wire value 

i 

Uncorrected value 

r 

Recovery value 

s 

Wetted length, shoulder value 

t 

Tunnel total ( stagnation) value 

u 

Upper frequency limit 

w 

Average hot-wire heated temperature or 
wall value 

1 

Outer edge of forebody boundary layer 

2 

Outer edge of free shear layer 


The bases of launch vehicles and reentry bodies 
are exposed to hot gases recirculated from the bound- 
aries of adjacent burning or ablating flows. The at- 
tendant high base heating rates have been the cause 
of mission failures. Accordingly, during the last 
five years we have conducted some simple base heat- 
ing experiments designed to afford direct insight into 
the convective transport processes. 

The need for a simplified physical model of base 
heat transport was established in early attempts to 
scale the results from tests of small models to flight 
conditions. Today base heating predictions still uti- 
lize small scale tests since neither the high flow en- 
thalpies nor the low pressures that cause recircula- 
tion through plume impingement shocks can be real- 
ized on a large scale. 

Any analysis of convective transport in recircu- 
lation zones, such as base heating, requires a knowl- 
edge of the flow properties inside the recirculation 
region. Most published experimental results concen- 
trate on the adjacent jet boundary areas instead of 
the base approach flow, since this is sufficient for 
the accepted base pressure theories in which the 
recirculation zone is replaced with a "dead air" 
region. Our experiments are an attempt to elucidate 
features of the flow within the recirculation region. 
The flow toward the base is compared with a two- 
dimensional stagnation point flow, hi analysis of 
possible effects of the high turbulence levels , the 
Hiemenz -Sutera [1-3] solutions of the Navier -Stokes 
equations were used as a laminar reference. 

Mean velocities, root -mean -square velocity 
fluctuations, total temperatures, base pressure dis- 
tributions , and heating rates have been measured in 
the recirculation zone downstream of the blunt trail- 
ing edge of a two-dimensional model in a Mach 3 
airstream (Fig. 1). 

The measurements proved to be extremely dif- 
ficult because of the low densities (0.01 of atmos- 
pheric) and high turbulence levels (in excess of 100 
percent). Moreover, the geometry is different from 
that of actual engine clusters or reentry bodies. 
Therefore, the results cannot be used quantitatively. 
However, they do allow interesting qualitative con- 
clusions about heat transfer characteristics of the 
base boundary layer and the occurrence of heat 
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studies of supersonic separation bubbles, flow over 
protuberances, wake flows, and cavity flows. % 


II. EXPERIMENTAL EQUIPMENT AND 
TECHNIQUES 

A. TEST FACILITY 


a) Model Installation in Wind Tunnel 




^38/1 CA.Lie.EK TANGENT 061 ye 




b) Forebody Side View 


The two-dimensional base flow experiments 
were conducted at the University of Minnesota Aero- 
Hypersonic Laboratory in a continuous -flow wind 
tunnel [4, 5]. Air at room temperature (T fc = 509° R) 

can be delivered at a maximum continuous rate of 
23. 2 lb/sec at stagnation pressures up to 100 psig. 
The test channel in which the experiments were 
carried out consists of a 54 -inch diameter stagnation 
chamber with provisions for flow straightening and 
reduction of turbulent velocity fluctuations. Mach 3 
steel nozzle blocks produced uniform flow within a 
6 x 9 -inch test section. 

B. WIND TUNNEL MODEL 



c) Schlieren Photograph of Flow Field 

FIGURE 1. WIND TUNNEL MODEL AND 
FLOW FIELD 

transfer effects arising from intense turbulence. In 
addition, the flow field is a special case of supersonic 
flow separation and reattachment, and our results 
might therefore be utilized in other fundamental 


Figure la shows the model installed in the 
wind tunnel. The model is two-dimensional, spans 
the test section, has a forward portion whose cross 
section is a tangent ogive, and has an 8. 19 -inch chord 
(Fig. lb). Separate bases and afterbody components 
were used for the pressure, heat transfer, and prob- 
ing runs. The base plates for all configurations have 
dimensions of 1.5x6 inches. 

The base of the pressure model is a polished 
nickel plate, 0.018 inch thick, instrumented with 
eleven pressure taps and thirteen thermocouples 
(Fig. 2a). 

The heat transfer base is made of pure nickel 
0.018 inch thick. Copper -constantan thermocouples 
and pressure taps are imbedded in the base at se- 
lected locations. Figure 2b illustrates the complete 
instrumentation on the base. The base and coolant 
passages are housed in a micarta afterbody to isolate 
the base thermally from the model and improve the 
accuracy of the heat transfer data. 

Although the micarta afterbody effectively iso- 
lates the base from the forebody, the possibility of 
conduction effects on the base itself had to be con- 
sidered. Bench tests with a small replica of the 
base plate showed that slots 14 mils wide, filled with 
Armstrong cement, formed an effective barrier to 
conduction heat transfer. Following these tests, the 
thermocouples on the base plate were insulated by 
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Ol) PRESSURE 5*3>E PtATE 



- I.5D0 



0.t>!8 Tw \CK niCkE! 
Base. Plate 


TYP 1C AU 
THERMOCOUPLE 
INSTALLATION 


b) HEAT TRANSFER bAS£ PLATE 



FIGURE 2. BASE PLATE CONFIGURATIONS 




means of the cement -filled slots shown in the figure. 
The vertical slots were added to minimize horizontal 
transfer of heat across the base plate by conduction 
through the brass pressure-tap extensions. 

The base of the probing model is of highly pol- 
ished stainless steel, 0.062 inch thick (Fig. 2c). 

In addition to the nine probing openings there are 
fourteen static pressure taps. These consist of a 
vertical centerline array of ten taps, with four ad- 
ditional taps on the horizontal centerline. Included 
also are several thermocouples attached to the inner 
surface of the plate. 


combined with Equation (i) gives 

dV^b 
fa)C p dA at 

T r’ T b ' 


( 3 ) 


The heat transfer coefficient h can then be obtained 

by measuring the instantaneous wall temperature and 

its rate of change. The wall recovery temperature 

T can either be calculated or measured as a function 
r 

of the tunnel total temperature. 

D. HOT-WIRE EQUIPMENT 


C. HEAT TRANSFER TECHNIQUE 

The transient technique was used to obtain 
base heat transfer values. In this method the model 
is initially cooled below the equilibrium wall temper- 
ature to rise toward the equilibrium value. 

Cooling lines were installed inside the model to 
bring the base down to 140 °R at the start of a run. 

A shield was placed behind the model to retard aero- 
dynamic heating until proper flow had been established 
in the tunnel. 

If heat conduction through the model skin and 
radiation effects are neglected, the rate at which heat 
is transferred from the boundary layer at any point 
on the model surface may be equated to the rate at 
which heat is stored at the same location, or 

3T 

q dA = cj c dV -r-^- (1) 

p 3t 

where the notation is the following: 

q = Rate of heat flow (Btu/ft 2 hr) 

dA = Differential area element of base (ft 2 ) 

co = Specific weight of base material (lb/ ft 3 ) 

c = Specific heat of base material (Btu/lb°F) 

P 

dV = Differential volume element of base 
material (ft 3 ) 

T = Base temperature (°R). 
b 

The Newtonian heat flow equation 
q = h(T r - T b ) (2) 


Probe calibrations and mean velocity meas- 
urements in the base recirculation zone were initially 
performed with a constant -temperature hot-wire 
anemometer [6] which was manually balanced. 

Fluctuating velocities were measured with a 
DISA type 55A01 constant -temperature anemometer 
without linear izer. A simplified circuit diagram of 
this anemometer and associated apparatus is shown 
in Figure 3a. 

The sensing elements were made either of 
0.0 0025 -inch diameter platinum -iridium (80 Pt. 20Ir) 
or of platinum plated tungsten wire with a 0. 000 2 -inch 
diameter (Fig. 3b). They were always aligned in the 
spanwise flow direction; that is, parallel to the z axis. 

The use of hot-wire anemometry proved to be 
difficult because of the low density (0. 01 of atmos- 
pheric) and the very high turbulence levels (« 100 
percent) . For our flow conditions the calibration 
results are believed to be representative of what can 
be obtained with commercial constant -temperature 
hot-wire sets at the present time. 

The presence of a wall increases the heat loss 
from the wire through conductive and other transport 
phenomena. Close to the wall the corrections are very 
large, as the heat loss to the wall overpowers the con- 
vective loss to the air. Furthermore, the usual pro- 
cedure of measuring the wall effect in still air as the 
difference between the readings, [E^ (x) ] 2 (close to 

the wall) and [E b (x— *>)] 2 (far from the wall) , 

yielded inconsistent results. Therefore, we estimated 
the wall effect from the difference between the indi- 
cated mean velocity profile and an extrapolated mean 
velocity profile. Figure 4 illustrates the procedure 
and shows typical results. The following adjustment 
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(parallel shift) was used for all static hot-wire cali- 
bration .curves: 

AA = [E b (x)] 2 - [E fa <x-«)] 2 

= 1. 3 x 10 4 (V 2 -in. )/x (in. ) (4) 




FIGURE 4. WALL EFFECT CORRECTION 
PROCEDURE 

The wall effect correction was applied at all 
probing stations; it was also used to correct the RMS 
fluctuation readings , even though the wall effect had 
been obtained from mean values. 

The low density in our tests not only produced 
an unusually large wall effect, but it also limited 
the frequency response severely. This was revealed 
by dynamic calibration tests. The hot-wire was 
placed in the smooth flow of a calibration duct and 
excited by a square or sine wave of heating current 
to simulate fluctuating velocity effects in the ane- 
mometer circuit. The hot-wire voltage transfer 
function was obtained over a frequency range utiliz- 
ing the circuit and equations described by Janssen, 
et al f 7 ] and Fingerson [8] for the sine wave meth- 
od. The circuit and equation employed are shown 
in Figure 5b (See following page). The voltage 
transfer function was measured over a range of den- 
sities, mean velocities, and wire overheats. The 

best estimate for the reference density of p, « 0. 01 

b 

atmospheric and zero velocity is given in Figure 5a, 
which shows the sine and square wave results for a 
platinum-iridium probe. Included also are square 


wave and estimate sine wave responses for a 
tungsten probe. Additional sine wave measurements * 
for tungsten probes are now in progress. 

A detailed square wave response over a density 
range of zero velocity for a tungsten probe is shown 
in Figure 6a, and Figure 6b describes the measure- 
ment procedure. 


At [fi sec) - Ax (cm) • Sweep Rote 
u 2 r At 




(a) Zero Velocity Cutoff Frequency (3<Jbdo*n) 
ot Various Densities 


(c) Idealised Wave Form 


FIGURE 6. SQUARE WAVE TEST 

Tests performed over a wide velocity range 
(up to 800 ft/sec) show little improvement in fre- 
quency response at our low test density. This is in 
contrast to the marked improvement obtained at 
ambient densities. 

The differences (up to 100 percent) between the 
sine and square wave tests indicate that the measured 
power transfer function may be correct in order of 
magnitude only. Further comparisons of these two 
methods are required. However, it indicates clearly 
that the range of usable frequencies extends only to 
approximately 4 kc/sec, whereas the wind tunnel test 
data show fluctuations above 30 kc/sec. This indi- 
cates that caution must be shown in interpreting low- 
density hot-wire measurements of velocity fluctua- 
tions containing frequencies beyond 4 kc/sec. 


in. BASE APPROACH FLOW 


The base approach flow has been studied with 
and without base cooling over a range of forebody 
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Reynolds numbers. Hot-wire, pressure, and tern- disturbs the flow field; therefore, later meas- 

perature measurements were performed at various urements were performed by insertion from the base 

locations in the shear layers and recirculation zone. plate. Some probings extended deep into the jet shear 

The location of these various probings are shown in layers and the f1 subsonic hole" between the plume 

Figure 7. The early measurements were performed impingement shocks. All hot-wire traverses have 

by insertion of probes either through the wake throat been restricted to our reference "test conditions"; 

or across the shear layers. However, this method that is, to the uncooled base at Re = 4 x 10®. 
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FIGURE 7. LOCATIONS OF PROBING MEASUREMENTS 


A. HOT-WIRE MEASUREMENTS 

Mean velocities and RMS velocity fluctuations 
along the centerline are shown in Figure 8. Con- 
trary to a two-dimensional stagnation point flow, the 
mean velocities are not constant outside the boundary 
layer but show a broad maximum at about 0.35 in. 
from the base plate. The RMS velocities are com- 
parable to the mean velocities; thus, we have an ex- 
tremely high local turbulence level, which makes the 
hot-wire measurements questionable. However, the 
absolute magnitude of the velocity fluctuations along 


the centerline are still considerably smaller than the 
velocity fluctuations in the free shear layers , which 
are estimated at 20 percent of the free stream ve- 
locity, or 420 ft/sec. 

Oscilloscope traces of typical hot-wire output 

voltages e, 1 are shown in Figure 9, where narrow 
b 

band components are evident. The amplitudes of the 
peaks are of the same order as the mean bridge 
voltages, indicating the possibility of flow reversal 
or negative velocity across the hot-wire, although 
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FIGURE 8 . VELOCITY VARIATION ALONG 
CENTER LINE 

the anemometer cannot discriminate between velocity 
directions. Consequently, all mean velocities are 
probably overestimated while RMS velocities are un- 
derestimated, even if one would correct for frequency 
response. 

In spite of these difficulties, we attempted a 
spectral analysis of the hot-wire output fluctuations 
e^* to obtain a confirmation of the narrow band com- 
ponents shown by the oscilloscope traces. To obtain 
longer records of the fluctuations and allow accurate 
calculations of pcwer spectra, a traverse through the 
approach flow region and base boundary layer was 
recorded on a Hewlett -Packard /Sanborn Model 3900 
magnetic tape data recording system. The data were 
analyzed with a third octave filter and yielded the 
power spectra shown in Figure 10. An approximate 
correction of the amplitudes for the limited frequency 
response was performed utilizing Figure 5a. The 
power spectrum reveals five distinct peaks at f = 0 . 2 , 
0.4, 1, 6 , and 30 kc/sec. The higher frequency 
peaks do not normally occur in low velocity recir- 
culation areas. We surmise, therefore, that they 
reflect acoustic transport of high frequency 
phenomena from the adjacent high speed points. 

Although somewhat speculative, there appear to 
be two possible sound radiation phenomena which 
could account for the narrow band frequencies ob- 
served. The first is connected with reverberating 
waves. Subsonic velocities have been indicated by 
pressure probings, interferometer evaluations, and 
wind vanes between the roots of the two plume im- 
pingement shocks. The associated region of sub- 
sonic flow extends downstream of the stagnation point 
that terminates the dividing streamlines and provides 


a f, subsonic hole. " Sound may be radiated backwards 
through this hole toward the base. These sound waves 
will then be reflected at the base and might support 
waves reverberating between the base and the "sub- 
sonic hole. " The occurrence of standing acoustic 
waves would then lead to isolated peaks of the pres- 
sure power spectrum which should be centered ap- 
proximately at the resonant frequencies of the 
corresponding "open pipe," 

a, 

f n " «T (2n + 1) (5) 

where n = 0 , 1 , 2 . . . denotes the fundamental and 
higher harmonic frequencies, L is the distance be- 
tween the base and downstream stagnation point, and 
a^ is the speed of sound in the base region. Substi- 
tuting the experimental values 

L = 1.3 in. 

a^ = 1070 ft/sec 

we obtain 

fo = 2.5 kc/sec 

f t = 7.4 kc/sec. 

The second harmonic f t is fairly close to the middle 
peak which was observed at 6 kc/sec. This frequency 
was found throughout the recirculation zone, which 
also substantiates the speculation of standing sound 
waves. 

The second type of sound radiation would origi- 
nate in the free shear layer. For example, Ffowcs- 
Williams and Maidanik’s description of Mach wave 
sound emission [9] might indicate a mechanism by 
which the high frequency fluctuations of the free 
shear layer could propagate into the recirculation 
zone. In free shear layers a broad maximum in the 
power spectrum is centered around a Strouhal number 
of 0. 2. This would lead to a frequency of 
U 2 

f « 0.2 — ( 6 ) 

<52 

where U 2 is the velocity at the outer edge of the free 
shear layer and 62 is the average geometrical shear 
layer thickness. Substituting the measured free 
stream velocity, U 2 = 2100 ft/sec, and the stream- 
wise average of the geometrical jet shear layer 
thickness, 62 = 0. 2 in. , one obtains f « 25 kc/sec. 
This agrees roughly with the highest observed peak 
in the spectrum of f = 30 kc/sec. 
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a) x = 0. 1 5 5 in. 


Re 


4 x 10 


T b ^ T r 

(At) h.w = 750 ° F 




b) x — 0.078 in. Probing Base Plate at 
Position 9. See Figure 
2c and Figure 7 



c) x - 0. 019 in. Vertical Sensitivity, 0. 05 V/cm 

_4 

Sweep Speed, 5x 10 sec /cm 


d) x = 0. 004 in. 


FIGURE 9. OSCILLOSCOPE DISPLAYS OF NARROW BAND VELOCITY FLUCTUATIONS 


If these speculations are borne out by future ex- 
periments, we may have serious noise problems on 
the upper stages, partially produced by plume im- 
pingement effects (subsonic hole) such that they can 
not be studied by full-scale static firings or single 
jet experiments. 

The low-frequency peaks cannot be explained by 
the sound radiation sources just discussed. Other 
simple or plausible explanations, such as the re- 
circulation of a "frozen" flow pattern, were not 
fruitful. 


B. TOTAL TEMPERATURE MEASUREMENTS 

The total temperature in the approach flow 
was measured with a thermocouple rake inserted 
from the tunnel side wall. The wires were aligned 
in the spanwise direction (z) to minimize conduction 
errors. 

The total temperature distribution in the approach 
flow was measured for several base cooling rates, and 
shown in Figure 11. Even for large amounts of cool- 
ing the base thermal boundary layer thickness is less 
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fower Spectrol Density, (V 2 - sec) 
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Tungsten Hot-Wire Probe 


(A T ) hw = 392® F 

Re = 4 x 10 s 
s 

x = 0.0194 in. 

Uncorrected 

Corrected for frequency 

response (Fig. 5) 

Data affected by filter settings of 

anemometer and tape recorder 
response 

FIGURE 10. VARIATION OF POWER SPECTRAL 
DENSITY WITH FREQUENCY 

than 0. 10 in. In spite of this, the approach flow 
senses the depression in base temperature to some 
extent, resulting in the creation of an approach flow 
which has a constant temperature below the free 
stream recovery value ( « 0.90 T^) . This is due 

partially to the recirculating flow motion which trans- 
fers cool fluid particles from the immediate vicinity 
of the base plate to other portions of the recirculation 
zone. We have, thus, the unusual situation that the 
effects of wall temperature are felt outside the ther- 
mal boundary layer. This is not a heat conduction 


effect, since the temperature gradients are zero 
outside the base boundary layer. Rather, it must 
result from a turbulent convection process that mix- 
es the "cooled" particles homogeneously within the 
recirculation zone, thereby creating a region of 
constant total temperature. 



0 0.02 0.04 0J 0 2 0 3 0 4 0.3 0.6 0.7 0.8 


FIGURE 11. TEMPERATURE IN APPROACH 
FLOW AND BASE BOUNDARY LAYER 

To show the interdependence of base plate tem- 
perature T^ and recirculation temperature T fi = 

T B (*,y, z ) more clearly, the latter has been replotted 

as a function of the base temperature (cooling rate) . 
Three curves are shown in Figure 12, with position 
as a parameter. The data which average out to the 
upper curve were measured throughout the recircu- 
lation zone. The fact that they fall on one curve 
demonstrates once again that the total temperature 
is constant throughout the base approach flow. The 
second curve shows a location close to the wall and 
here the scatter of the points is probably caused by 
a slight misalignment of the total temperature rake. 
The point of intersection at the three curves indicates 
the recirculation zone and base recovery temperature , 
which is approximately 95 percent of the free stream 
stagnation temperature T^. 

The coupling between the base temperature and 
the "remote" approach flow is one effect which can- 
not be treated by a stagnation point analysis where 
the upstream conditions are unaffected by wall cool- 
ing. Furthermore, the tT base boundary layer" is 
affected at the shoulders , since the cooled particles 
must be turned back to be recirculated. However, 
the hot-wire and total temperature traverses show 
that both a velocity and a thermal boundary layer 
exist at all probing stations. 
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Mi 
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A Temperature probe 

FIGURE 12. TEMPERATURE DISTRIBUTION 
IN APPROACH FLOW 


deviations from a two-dimensional behavior. The 
decrease of the wall pressures along the two- % 

dimensional "stagnation point line" y = 0, x = 0 in- 
dicates that the deceleration of the base approach flow 
resembles that of an axisymmetric flow instead of a 
two-dimensional flow. The following discussion of 
the pressure measurements is limited to the plane of 
symmetry z = 0. 

Falkner and Skan [10] have shown that in bound- 
ary layer theory, the existence of an outer "potential" 
flow, which follows from the wall pressures, can be 
extended to stagnation point flows if one compares 
the calculated velocities with the tangential velocity 
component. The wedge solutions of Falkner and Skan 
for stagnation flow coincide with Hiemenz* exact so- 
lution of the Navier -Stokes equations. In laminar 
stagnation point flows, the wall pressure distribu- 
tions can thus be used to calculate the tangential ve- 
locity components at the outer edge of the base bound- 
ary layer. This has been done assuming the following: 

a) The stagnation pressure along the edge is 
constant and equal to the wall pressure at 
y = 0. 

b) The Bernoulli equation may be used along 
the outer edge in spite of the high local 
turbulence leveljs. 

Within these assumptions the mean tangential 
(vertical) velocity component at the outer edge 
follows from 

v(0,y,0) = ~ tP(0,0,0) -p(0,y,0)]i(7) 


The results of this calculation over a Reynolds num- 
ber range are shown in Figure 14. For Reynolds 
numbers of 2.78 x 10 6 and larger a linear velocity 
variation exists up to y =0.3 in. This is slightly 
below position number 9 (probing station 4P) shown 
in Figure 7. The scatter at the lower Reynolds 
numbers is attributed to pressure measuring errors 
and transitional flow conditions. 


IV. THE BASE BOUNDARY LAYER 


A. VELOCITY BOUNDARY LAYER 

The measured wall pressures are summa- 
rized in Figure 13. The existence of a stagnation 
point is shown by the pressure maximum around the 
centerline. However, the pressure also varies in 
the horizontal z direction. Accordingly, we observe 


The consistent trend of increase in the stagnation - 
point velocity gradient with Reynolds number is also 
shown in Figure 14. A comparison can be made 
between these data and values obtained from cen- 
terline hot-wire probings, utilizing the continuity 

equation. From Figure 8, at Re = 4 x 10 6 the 

s 

centerline velocity gradient BU/3x is approxi- 
mately equal to 6 x 10 3 sec "**, which compares 
with 7. 4 x 10 3 sec shown in Figure 14. 
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FIGURE 13. RASE PLATE PRESSURE DISTRIBUTION 


9U/dx is approximately equal to 6 x 10 3 sec 

which compares with 7.4 x 10 3 sec shown in Figure 

14. 

Comparing the "outer edge" velocities deter- 
mined from the wall pressure distribution with mean 
velocities measured by the hot-wire, one apparently 
finds little agreement. However, if one interprets 
the non-zero value at y = 0 as resulting from the 
axial velocity component u, reasonable agreement 
exists up to y = 0. 160 in. Consequently, the usual 
boundary layer assumptions seem approximately 
justified up to this position. The strong disagree- 
ment at higher stations cannot be presently ex- 
plained. However, because of the extremely high 
local turbulence levels on the base approach flow, 


one might expect the base boundary layer to exhibit 
some new or unexpected characteristics. 

Some qualitative features of the adjustment of 
this boundary layer to the high turbulence levees of 
the base approach flow can be inferred from hot-wire 
traverses. RMS velocity profiles are shown in Fig- 
ure 15. The absolute values show the usual IT hump" 
attributed to turbulence amplification near a solid 
wall. Sternberg [11] has analyzed this amplification 
in terms of "image" vortices which replace the wall 
in potential flow theoiy; whereas Sutera, et al. [2,3], 
have studied the growth of vorticity which attends the 
selective stretching of vortex lines in a two- 
dimensional stagnation point flow. This vortex 
stretching hypothesis is supported by the measure - 


15 






Vertical Coordinate , y (in) 



• =Re = 4 x 10 6 Hot-Wire Measurements 
s 

Re x 10- 6 
s 

0 1. 6? O 5. 02 

□ 2. 19 Q 5. 40 

02.78 <>6.02 

<*3.35 C 6.43 

L3. 78 k 7. 02 

b4.40 No Cooling 

FIGURE 14. BASE PLATE TANGENTIAL 
VELOCITY COMPONENT AND STAGNATION 
POINT VELOCITY GRADIENT 


ments of Kuethe, et al. [12] , in the stagnation region 
of a sphere, and also by traverses of the flat plate 
boundary layer (Townsend [13] and Sternberg [11]). 
However, the turbulence levels in the base boundary 
layer are much higher. In flat plate boundary layers 
they reach about 60 percent, whereas our measured 
values go up to 200 percent. The actual levels are 
probably even higher owing to insufficient frequency 
response and the hot-wire f s inability to distinguish 
flow direction. 

We have tried to assess the possible effects of 
turbulence by comparing hot-wire and total temper- 
ature traverses with the laminar stagnation point 
solutions [1-3], These solutions depend on an em- 
pirical parameter, the velocity gradient 3tr/3y at the 
outer edge of the stagnation point boundary layer. 

This gradient was obtained from the wall pressure 


measurements as shown in Figure 13. A strong 
Reynolds number dependence is evident, leading to 
the suspicion that it is also dependent on the turbu- 
lence level. 

The chosen velocity gradients therefore reflect 
the turbulence of the actual flow such that the lami- 
nar stagnation point solutions are automatically ad- 
justed to give the full turbulent shear stress. 

The ’’adjusted” laminar mean velocity profiles 
were calculated for these traverses where the meas- 
ured streamwise gradient of the outer edge velocity 
remains approximately constant (Fig. 14). Figure 
16 shows a comparison between these calculated mean 
velocity profiles and hot-wire results. Included at 
positions 9 and 10 (see Figure 7, probing stations 4P 
and 6P) are calculations based on Sutera’s solution 
for various values of the amplification parameter A. 
The agreement is rather poor. 

One may tentatively conclude, therefore, that 
turbulence can appreciably change the shape of the 
non-dimensional mean velocity profiles. 

B. THERMAL BOUNDARY LAYER 

Total temperatures were measured with a 
miniature thermocouple probe (0.004-in. face height) . 
The results are shown in Figure 17, which includes 
the adjusted two-dimensional stagnation point solu- 
tions from the theory of Livingood and Donoughe [14] . 

The theoretical profiles assume laminar flow and 
apply the same empirical velocity gradient already 
used in the calculations for Figure 16. The wall gra- 
dients are less than the measured profiles, indicating 
that the high turbulence level changes the shape of the 
non-dimensional temperature profile. The dashed 
curves were obtained from measured heat transfer 
values, showing good agreement with the experimental 
points. The temperature gradients at the wall are 
higher than those predicted by the ’’adjusted” laminar 
stagnation -point analysis, even though the full velo- 
city gradient (and shear stress) of the turbulent flow 
was taken into account. 

The theoretical and experimental temperature 
gradients are replotted in Figure 18 as a function of 
the temperature difference between the base approach 
flow and the base. Both gradients increase faster 
than linearly with this temperature difference; how- 
ever, the measured gradients exceed the estimated 
ones by a factor of 2. The same is true for the heat 
transfer rates which are shown in the lower part of 
Figure 18. They were obtained from the temperature 
gradients by multiplication with the thermal conduc - 
tivity k of air. The difference in shape of the heat 
transfer and temperature gradient curves is caused 
by the temperature dependence of k which was evalu- 
ated a t the base plate temperature. 
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FIGURE 15. RMS VELOCITY PROFILES IN TWO-DIMENSIONAL BASE FLOW 
INCLUDING CORRECTION FOR WALL EFFECT 


VI. HEAT TRANSFER 


Heat transfer rates were measured as a function 
of location (coordinates 0,y,z) , the base plate tem- 
perature T^, and the forebody Reynolds number Re g . 

All measurements utilized the transient technique 
discussed in Section II C. 

A typical area distribution of the heat transfer 
rate is shown in Figure 19. The measured values 
vary between 1900 and 2400 Btu/ft 2 hr, or within 
±12 percent of the mean. The heat transfer rates 


reach a maximum at the stagnation point (thermo- 
couple 3) , and fall off in both the y and z directions , 
except for the vertical traverse through thermocouple 
positions 7 and 9. This behavior is to be expected 
because of the deviations from two-dimensional stag- 
nation point flow, although the adjusted laminar 
stagnation -point calculations predict a constant heat 
transfer rate if a linear velocity variation exists. 


The total heat transfer coefficient 
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FIGURE 18. THE EFFECT OF TURBULENCE LEVEL ON BASE HEAT TRANSFER 
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FIGURE 19. DISTRIBUTION OF BASE HEATING RATES 
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is shown in Figure 20. We have tried to suppress 
* dependence on position by using the area -averaged 
heat transfer rate q A only. The effects of base plate 

temperature are eliminated by replacing the heat 

transfer rate with the heat transfer coefficient h, - 

tot 

If the heat transfer rate q increases linearly with 

A 


increasing temperature difference ( T - T) , h 

r b tot 

should be independent of the cooling rate. Actually, 
the variation with base plate temperature is 
slightly nonlinear as shown in Figure 18, which 
explains the cluster of points at each Reynolds 
number. 


h !0t 




V T b 


T. * 


T. » 


60° F 
520° R 

.95 T, (°R) 
494 °R 


Heat Tronsfer Averaged Over 
Instrumented Area 


Smooth 

Surfoce 


Rough 

Surface 


T b (°R) 
360 



FIGURE 20. VARIATION OF TOTAL HEAT TRANSFER COEFFICIENT WITH REYNOLDS NUMBER 


A. HEAT FLOW RESISTANCE OF FREE SHEAR 
LAYER AND BASE BOUNDARY LAYER 

All of the heat flowing into the base plate 
must come from adjacent jets. The transport of this 
heat is by conduction through the free shear layer , 
conduction through the base boundary layer, and con- 
vection in the recirculation process. However, the 
conductive transport can be analyzed separately if we 
let the dividing streamlines define the boundaries of 


a control surface for heat transfer ’’through" the jet 
shear layer. Since the net mass transfer across 
each surface element is zero, there is no net con- 
vective transport. Consequently, the convective re- 
circulation process operates only internally to create 
a base recirculation zone of constant temperature; 
that is, it establishes a region of a very high "appa- 
rent" heat conductivity. The heat which flows into 
the base must also be conducted through the two di- 
viding streamlines. We can draw an analogy to 
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electrical current flow by assuming the jet shear 
layer and the base boundary layer act as two heat 
flow "resistances" in series, connected by a region 
of high heat conductivity whose resistance may be 
neglected. To estimate the magnitude of these two 
resistances > we define new heat transfer coefficients , 
h gL and h^, for the jet shear layer and the base 

plate boundary layer. Then 


= h (T - T ) = h (T - T ) 
SL 1 r B ; BP 1 B V 


= h 


tot 


<T ’ T J 
r b 


(9) 


where 


for Re = 4. 5 x 10 6 . Substitution into the previous 
s 

equation gives 



9T b 

=0.857. 

b T = T 
B r 


(13) 


Accordingly, 86 percent of the total heat flow re- 
sistance can be attributed to the thermal base bound- 
ary layer, which clearly indicates its dominance in 
the base heating process. The remaining heat flow 
resistance is provided by the jet shear layer, such 
that 


T = recirculation zone and base recovery 
temperature 

T = recirculation zone temperature 

-D 

T^ = base plate temperature. 


h. , h , 

tot _ tot 

h SL = “ h BP 


dT 

( 8T 2_) =0.143,(14) 

b T = T 
B r 


The foregoing results have been obtained for one 

reference Reynolds number, Re = 4.5 x 10 6 . To 

s 


This concept of two heat flow resistances in 
series has also been used recently by Page and 
Dixon [18] , although they refer to the second re- 
sistance as that provided by the "recirculating flow. " 
In a mathematical description, the term "resistance" 
is now identified with the reciprocal of the heat trans- 
fer coefficient ("conductivity" ) . 

‘'No, - 1/h SL* 1/h BP • “»> 

which uses the additive property of series resistances. 
A comparison between the resistance of the base 
boundary layer and the total resistance is then given 
by 


show the heat transfer coefficient of the base bound- 
ary layer and the free shear layer as a function of 
Reynolds number, we have repeated this calculation 
for other test conditions, converting to the non- 
dimensional Stanton number 


St - 


tot 


P2 U 2 C 


p2 


(15) 


where the density, velocity and specific heat (p 2 , 

U 2 , c ) have been evaluated at the outer (high speed) 

pz 

edge of the jet shear layer. This enables us to com- 
pare our free shear layer heat transfer coefficient 
with those of Korst, et al. [ 19] . 



To obtain an average value, independent of the base 
cooling rate, we shall approximate the measured de- 
pendence of T fi on T fe (Fig. 12) with the straight line 


0T 


T = T + 

B r 9T, 


B 


(T h “ 
b r 


c T + 0. 143 (T - T ) 
r b r 


(12) 


The expansion process of the forebody boundary 
layer is known to be essentially isentropic [20,21] , 
at least in the supersonic portions; so the calculations 
of the properties at the outer edge of the shear layer 
were based on the measured base static pressure 
ratio P^/Pi* using the assumption of constant total 

pressure p 2 = p*. This also assumed that the lip 
shock wave is of negligible strength. 

The static -to -total-pressure ratio 


5. = ^ = ( Pb 


t2 


t2 


Pi 


Pi 

)<— ) 

P tl 


(16) 


M=2. 95 
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can thus be calculated; and the Mach number M 2 , 
density p 2 and velocity U 2 at the outer edge follow 
from isentropic flow tables and the known stagnation 
temperature T^. 


The results of the free shear layer heat transfer 
data are summarized in Figure 21, which also shows 
Korst T s semiempirical results depending on the 
choice of the spreading parameter a. The plotted 



10 s 2 3 4 5 6 7 8 9 10* 2 3 4 5 6 7 8 9 10 T 

Reynolds Number, Re^ 


FIGURE 21. BASE PRESSURE VARIATION WITH REYNOLDS NUMBER 


values take account of the fact that Korst based his 

transfer coefficient h on the area of the dividing 
oL 

streamline, whereas we utilized the base area. The 
data agree fairly well for a = 12, the incompressible 
value. Korst recommends that a be calculated from 
the formula, 

a = 12.0 + 2.76 M 2 = 24.2, (17) 

suggested by Tripp [22] and based on values for jets 
exhausting into still air; but this evidently leads to 
large discrepancies. Because of the differences in 
the flow fields this value may not be applicable to 
our case. 


B. REYNOLDS NUMBER EFFECTS 

The predominant heat flow resistance of the 
base boundary layer and the extremely high^rbu- 
lence levels have been found at comparatively snbdl 
Reynolds numbers. In order to draw conclusions on 
base heating and acoustical loads in flight one would 
like to extrapolate the result to Ml scale Reynolds 
numbers. A first step in this direction is to con- 
sider the Reynolds number dependence over the tested 
Reynolds number range. Additional measurements 
at higher Reynolds numbers are presently prepared. 

The measured total , Stanton numbers appear to 
decrease slightly with an increasing Reynolds number, 
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as shewn in Figure 22. However, this trend is barely 
recognizable within the measuring accuracies and 
any extrapolation of flight Reynolds numbers is very 
questionable. 


In view of the extremely high turbulence levels, 
one might argue that an increase in Reynolds number ^ 
should not affect the heating rates very much. The 
Reynolds number independence of highly turbulent 


Stanton Number, St 



FIGURE 22. VARIATION OF TOTAL, SHEAR LAYER, AND BASE PLATE HEAT TRANSFER 

COEFFICIENTS WITH REYNOLDS NUMBER 


flows has been established for mean velocity profiles 
in true shear layers and above rough walls. In 
such cases the shear stresses are directly propor- 
tional to the dynamic pressure pU 2 . The linear 
variation of the stagnation point velocity gradient 
with Reynolds number (Fig. 14) now suggests that 
the shear stress in the base boundary vanes with the 
stream density pu rather than the dynamic pressure. 
The base boundal^y layer therefore exhibits an un- 
expected Reynolds number dependence and any ex- 
trapolation beyond the tested Reynolds number range 
remains questionable, even if one uses physical 
analogies. 

The dimensional heat transfer rates (Fig. 2) do 
show a significant increase with Reynolds number. 
This Reynolds number dependence is provided almost 


entirely by the Reynolds number dependence of the 
recirculation base density, since the Stanton number 
itself is almost independent from the Reynolds num- 
ber. For the uncooled base the density and pressure 
are directly proportional. The Reynolds number de- 
pendence of density may therefore be evaluated from 
the plot of base pressure versus Reynolds number , 
which is given in Figure 22. The measurements 
indicate the typical "jump" which is known from 
laminar turbulent transition and may be attributed to 
transition effects of the forebody boundary layers. 

In actual flight, the forebody boundary layers 
(nozzle shear layers) are far from transition. This 
is an additional reason why the heat transfer rates 
from small scale model tests should not be extra- 
polated to full scale Reynolds numbers, in spite of 
the light turbulence levels. 
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V. CONCLUSIONS 


REFERENCES 


Heat transfer processes in a two-dimensional 
recirculation zone behind a blunt trailing edge in a 
Mach 3 wind tunnel have been studied by hot-wire and 
total temperature traverses combined with wall pres- 
sure and heat transfer measurements. The results 
indicate: 

1. The recirculating flow is subjected to ex- 
tremely high local turbulence levels of « 100 percent, 
which makes a quantitative hot-wire interpretation 
difficult. 

2. Oscilloscope traces and a third-octave fre- 
quency band analysis indicate at least five narrow- 
band frequency components in the base approach flow. 
We speculate that one of these is connected with 
sound radiation emanating from the subsonic hole 
between the plume impingement shocks and the jet 
shear layers, while another arises from noise gen- 
erated by the shear layer. 

3. Cooling the base plate produces a constant 
total temperature throughout the recirculation zone, 
which is below the recovery temperature. The asso- 
ciated influence of the wall on the remote base ap- 
proach flow can not be duplicated with a stagnation 
point analysis. 

4. Turbulent reattachment produces thermal 
and velocity base boundary layers of approximately 
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boundary layer and the jet shear layer, which act 
like two heat flow resistances in series. The ther- 
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the total heat flow resistance and, therefore, domi- 
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8. The shear layer heat transfer coefficient 
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II. FLUID MECHANICS 


FLUID MECHANICS REGIMES RELATED TO PROPELLANT OSCILLATIONS 
UNDER LOW GRAVITY CONDITIONS 

By 

Harry J. Buchanan 


SUMMARY 


This paper presents a discussion of the basic 
regimes of fluid mechanics pertinent to propellant 
oscillation problems. In particular an explanation 
is given of some of the fluid phenomena observed 
under low gravity conditions in terms of certain 
dimensionless products which are used to define the 
regimes of fluid behavior: ordinary high gravity 
sloshing, ordinary low gravity sloshing, low gravity 
sloshing predominated by liquid moving along the con- 
tainer surface, and sloshing involving breaking of 
the surface by splashing or spraying. 

LIST OF SYMBOLS 

a Tank radius 

Bo Bond number 

F Force 

Fr Froude number 

g Q Acceleration due to gravity at sea level 

g Vehicle acceleration prior to main 

engine shutdown 

g Vehicle acceleration after main 

engine shutdown 

L Characteristic length 

r Radius measured from tank center- 

line 

V Velocity of fluid elements at free 

surface 

We Weber number 

P Fluid mass density 

& Fluid surface tension 


Slosh wave amplitude at tank wall 
w First mode slosh frequency 


INTRODUCTION 


In the past, problems relating to propellant 
oscillations in the tanks of launch vehicles have 
been concerned primarily with the boost phase of 
flight where body forces due to vehicle acceleration 
are fairly large. Recently, with the advent of 
liquid fueled space vehicles designed to restart after 
an extended period in a near weightless condition, 
interest in problems involving low gravity propellant 
oscillations has increased. As of yet little progress 
has been made toward solving these problems. One 
of the specific problems to which considerable 
attention has been turned involves the dynamic be- 
havior of liquids following a large reduction in 
vehicle acceleration. 

At main engine shutdown the space vehicle ex- 
periences a transition from a high forward acceler- 
ation to a low deceleration resulting from aero- 
dynamic drag. For ordinary tanks the fluid will then 
reorient to a position of minimum potential energy. 
Since the fluid will move to the other end of the tank 
when this acceleration sign change takes place, one 
design philosphy has been to prevent this by applying 
sufficient thrust to insure a small forward acceler- 
ation. Under these conditions, if the fluid were 
initially at rest, the free surface would merely ad- 
just its slope to minimize the potential energy. 

While the resulting motion might require consider- 
able time to damp out, it is not particularly violent. 
On the other hand, if the fluid is sloshing initially, 
which is certainly the most likely case, a consid- 
erable amount of energy is possessed by the fluid. 
How this energy will be divided into potential and 
kinetic energy depends upon the phasing of the main 
engine cutoff. If this event takes place when the 
kinetic energy of the fluid is large, a particularly 
violent motion can result. 
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FLUID DYNAMIC REGIMES 


In order to examine qualitatively the dynamic 
behavior of the liquid for the case described above, 
consider four basic fluid mechanics regimes for 
describing sloshing phenomena. Figure 1 illus- 
trates these regimes in terms of the Weber and 

Froude numbers. The Weber number (We = ^ ~ ) 

a 

represents the ratio of the inertial forces to the 

V 2 

capillary forces, and the Froude number (Fr = — ) 
is the ratio of inertial forces to body forces. The 

Bond number ( Bo - gL 2 p/a ) is equal to the Weber 
number divided by the Froude number. 



The resulting fluid motion in this regime has been 
thoroughly treated, both analytically and experi- 
mentally, by a number of investigators [1-4]. The 
slosh wave is essentially planar as indicated in 
Figure 2. 
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FIGURE 2. HIGH GRAVITY SLOSHING 


B. LOW GRAVITY SLOSHING 

In this regime capillary forces and body forces 
both are greater than inertial forces. In terms of 
the dimensionless products this condition can be 
described by 

Fr <1. 0 We <1.0. 

For these conditions the fluid oscillates or sloshes, 
but the resulting wave shape is entirely different 
from that for high gravity sloshing, as shown in 
Figure 3. A detailed solution for this motion is not 


FIGURE 1. SLOSHING FLUID MECHANICS 
REGIMES 

A. HIGH GRAVITY SLOSHING 

This is the familiar case experienced in tests 
conducted at 1 g and in flight tests at several g ! s. 
Here the body forces predominate over inertial 
forces, and inertial forces are very large compared 
to surface tension forces. In terms of the above 
dimensionless products, this regime is defined by: 

Fr <1. 0, We >1.0, Bo >1. 0. 
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FIGURE 3. LOW GRAVITY SLOSHING 
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yet available although a method for determining the 
natural frequency was developed [5], This region 
can be divided into two subregions by detailed con- 
sideration of the restoring mechanism. In one case 
the body forces provide the most important restoring 
force; however, capillary forces cannot be entirely 
ignored. This regime is described by 


F Inertial _ pV 2 a 2 

F o 

Body pa^g 


ag 


We 


F Inertial = pV 2 a 2 _ pV 2 a 

F _ aa a 

Capillary 


( 1 ) * 


( 2 ) 


We > Fr Bo >1. 0. 

In the second case capillary forces supply the primary 
restoring mechanism, for this situation 

We < Fr Bo <1.0. 


C. FLOW UP TANK WALL 

This is the regime which is apparent in the drop 
tests; it occurs where inertial forces predominate 
over body forces; however, capillary forces remain 
stronger than inertial forces. For this case 

Fr >i. 0 We <1.0 Bo <1.0. 


In order to evaluate the velocity, V, let us assume 
that we are interested only in the main engine cutoff 
phasing which produces the maximum kinetic energy. 
Since this corresponds to the maximum velocity, it 
will also produce the maximum inertial force. If 
the slosh wave prior to main engine cutoff can be 
assumed to be planar, the velocity of the fluid par- 
ticles at the surface becomes 

V = r/a cl) ( 3) 

w 


and, since the maximum velocity is at the wall, 


V 

max 


= t, Cl) . 


W 


(4) 


The Froude number can now be rewritten 

In this regime the fluid climbs the container walls, 

while the surface is prevented from spraying by the ^ 2 jz 

surface tension forces. Fr = — (5) 

ag 


D. SPRAYING 


and, in like manner, the Weber number 


This last regime is characterized by liquid 
spraying or splashing. This occurs when the capil- 
lary forces can no longer resist the inertial forces 
which already exceed the body forces. For this 
condition 

Fr >i. 0 We >1.0. 


p£ 2 co 2 a 

We = — . (6) 

(7 

By setting Fr = We = 1. 0 in (5) and ( 6) , it is pos- 
sible to arrive at the equations for the lines dividing 
the various regimes in Figure 1. 


The Bond number for this case can be either high or 
low. For high Bond numbers the spraying criterion 
has been derived [6] and is consistent with the 
criterion put forward here. 


From ( 5) 
Fr = i. 0 


From (6) 
We = i. 0 



EVALUATION OF DIMENSIONLESS NUMBERS 


RESULTS 


(7) 

( 8 ) 


Figure i and the suggested regimes are of little 
practical significance unless expressions for Fr and 
We can be derived which correspond to the ratios 
described. 


In order to compare the above explanation with 
the observed phenomena, consider the case shown 
in Figure 4. This sequence of photographs was 
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a. High Gravity Sloshing 
Prior to Drop 


b. Slosh Wave at Drop 
-Maximum Kinetic 
Energy 



c. Flow up Tank Wall d. Flow around Tank 

Completed 


FIGURE 4. FLOW UP WALL 
TANK DIAMETER: 3 INCHES 
FLUID: CARBON TETRACHLORIDE 
FACILITY: x.ORTH AMERICAN 
AVIATION DROP TOWER 


taken from a drop tower test conducted by North 
American Aviation and shows a model of one of the 
Apollo Service Module propellant tanks. The test 
procedure was to establish an initial slosh wave at 
1 g and then to release the package to produce the 
maximum kinetic energy. Since this tank is nearly 
cylindrical, the natural frequency prior to drop is 
given [ 1] by 


Cl ) 2 = ( IS) 


g/go 


(9) 


Then from (7) . Fr = L 0 


w 


= .74 • 


g/g o 

g /g o 


( 10 ) 


and from ( 8) , We = 1.0 



(ID 


Using equations (10) and ( 11) it is possible to de- 
fine the four regimes in terms of ^ and g. Figure 

5 was constructed using the dimensions and fluid 
properties for the model shown in Figure 4. The 
shaded region indicates the acceleration level ex- 
perienced in the North American test facility. From 
the photographs it appears that the wave height is 
sufficient to initiate spraying. However, these 
boundaries define only the onset of each of these 
regimes. Actually all the fluid participating in the 
motion must have a velocity equal to or greater than 
that associated with this wave height. Therefore, the 
wave height required to produce the motion pictured 
in Figure 4 would probably be greater than that indi- 
cated by Figure 5. It appears that the value of the 
Bond number determines which flow regimes will be 
possible regardless of wave height. Thus, according 
to Figure 5, flow up the wall is not possible if Bo> 1.0. 


The effect of tank size on the phenomena is 
demonstrated by considering a much larger tank 
such as the Saturn S-IVB liquid hydrogen (Fig. 6). 
The current plan for propellant control requires 

—3 2 

that an acceleration level of 5 x 10 m./sec be 
applied to the vehicle for 77 seconds following main 
engine cutoff [7]. Figure 6 indicates that the 
phenomena involved here include spraying or splash- 
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Wave Height Prior To Drop (mj Height Prior to S-IZB Cutoff (m) 




FIGURE 5. SLOSH REGIMES FOR APOLLO 
SERVICE MODULE PROPELLANT 
TANK MODEL 

ing but not flow up the tank wall as pictured in Figure 
4. Therefore.it would seem that in large tanks (large 
Bond numbers) the spraying or splashing phenomena 
would be more prevalent than that described as flow 
up the tank wall. 


From the discussion it seems possible to explain 
for at least some cases the mechanism which pro- 
duces the flow up the tank wall. Although the method 
for evaluating the dimensionless numbers is based 
on several simplifying assumptions, it does offer a 
means of estimating the possible flow regimes and 
of scaling experiments. Some subjects for future 
investigation include the effect of viscous forces and 
the influence of contact angle or adhesion between the 
fluid and wall. 


FIGURE 6. SLOSH REGIMES FOR SATURN S-IVB 
LIQUID HYDROGEN TANK 
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III. INSTRUMENTATION 



REMOTE SENSING WITH OPTICAL CROSS CORRELATION METHODS 


By 

7. R. Krause and M. J. Fisher 


SUMMARY 


Our basic concept is to employ the cross corre- 
lation of two narrow light beams for remote sensing. 
The problem with most optical methods such as schlie- 
ren, shadowgraph systems or spectrometers is that 
the received signal is always integrated along the line 
of sight. However, we want information for a partic- 
ular spot along the line of sight. To this purpose we 
adjust a second beam to intersect the first at the point 
of interest. The integration along the two beams is 
then eliminated by a cross correlation of the detected 
fluctuations. This paper indicates how this new con- 
cept of combining standard spectroscopy with a statis- 
tical cross correlation analysis might be used for 
turbulence investigations in model tests and in the 
atmosphere. 

The concept has been demonstrated success- 
fully by comparing hot wire measurements with opti- 
cal cross correlation results. Since remote sensing 
methods are not limited by scale, AAP applications 
are conceivable. The main problems in atmospheric 
and AAP experiments are formulated in terms of 
questions, such that the background information be- 
comes apparent which is needed to establish the 
feasiblity of particular experiments. 

I. INTRODUCTION 


The cross correlation of optical signals has been 
used for remote sensing in jet shear layers [ 1] . Con- 
centration fluctuations of a water vapor fog could be 
resolved with good accuracy in space and time [2] . 

In principle, the method should also work on an at- 
mospheric scale. Experimental investigations about 
the generation and motion of atmospheric constituents 
such as clouds, nuclear debris, ozone, radiation belts, 
etc. , offer an exciting challenge. We will now give a 
simplified version of our new remote sensing concept 
and outline some of the problems, which one would 
have to overcome in atmospheric applications. 


H. THE CROSSED BEAM CONCEPT 


Remote sensing of distant light sources like 
planets is based on a fundamental theorem of geo- 
metrical optics, which states that the spectral radi- 
ance, I, of a source and its image are equal, if one 
neglects light extinction and emission along the line 
of sight. The spectral radiance of a distant source 
is therefore accessible through the spectral radiance 
of its image on the exit slit of a monochromator, as 
shown in Figure 1. Source area, AA, and wave 
length interval, AX, are set by the optical filter. The 
solid angle, A -a, into which the exit slit radiates, is 
set by the photodetector. Therefore, the received 
radiative power is directly proportional to the spec- 
tral radiance, I, of the distant source: 



FIGURE i. REMOTE SENSING OF DISTANT LIGHT 
SOURCES 


I(x — " - , t , X) = I g (t, X) 

_ w (x = T m 

AA A -a AX d * ; 

However, if the line of sight goes through the atmos- 
phere, the detected spectral radiance ^(t, X) is 
changed due to the light extinction or emission in the 
detector T s field of view. In the two limiting cases of 
dominant extinction or dominant emission, the change 
in spectral radiance 
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Ij(t, X) = I(5T — + «>, t, X) 

+« 

-f K(x, t, X) dx 

= I(t,X)e _0 ° (2) 

may be expressed by integrating the extinction or 
emission coefficient 

K(x, t, X)= J =K<r, X) + k(r, t, X). (3) 


In astronomy and optical communication, these ef- 
fects are a nuisance. From an atmospheric science 
point of view, they can be used to study the thermo- 
dynamic properties of the atmosphere, since the ex- 
tinction coefficient K is uniquely determined by the 
thermodynamic properties. 

The problem with local measurement of atmos- 
pheric properties is how to retrieve the local infor- 
mation that is buried in the integrated signal. A so- 
lution is possible only by repeating the remote sensing 
as a function of the remaining independent variables , 
such that some mathematical transformation can be 
used to extract the local extinction coefficient. 

Several authors f3 - 8] used the remaining space 
variables y and z and had to assume axial sym- 
metry. This assumption can be made in the study of 
single rocket exhaust and flames but not in the atmos- 
phere. Another group [9 - 15] used the dependence 
on optical wave length X. They subdivide the line of 
sight into as many zones as there are isolated mole- 
cular bands , and need some advance knowledge about 
the spatial distribution of temperatures and concen- 
trations. The method worked on a package of con- 
trolled burners [16]. In the atmosphere, a larger 
number of isolated bands which are free from scat- 
tering and the advance knowledge are hard to find. 
However, the "zonal approximation" should be con- 
sidered as a backup of our cross correlation method, 
since it does use the same test arrangement. 

We use the time aspect to retrieve the local in- 
formation. Consider the experimental arrangement 
on Figure 2. Two lines of sight intersect in an at- 
mosphere-fixed point. These lines define a plane 
which we shall call "turbulence wave front" and a 
normal in the intersection point (x, y, z) which is 
called the wave normal. We then measure the spec- 
tral radiance I(t, X) (radiative power/solid angle, 
area, and wave length interval) along both fields of 
view, split off the time average signal, I, cross 
correlate the fluctuations, i(t) , and calculate the 
quantity, G, as described in Figure 2. Let us assume 


that the two detecting telescopes have been aligned to 
the y and z axis of a Cartesian, atmosphere -fixed 
coordinate system. Each line of sight is denoted by 
a y or z subscript and each point on the line of sight- 
by the distances 77 or £ relative to the intersection 
point. The quantity G is then equal to 


_ i (t, X) i (t, X) 

G(x , X) = — 2 


Detector I 


White Light 

The Crossed Seom Technique Arec Integrates the Cross 
Correlation Function of Light Extinction Coefficients, 
K=K + k 



Detector 2 

Iy = (Iy) +1 y(t) 


<i z (t) iy(t+r)) 

<I Z > <I 2 > 


+ CD 

= G(f ,r,Xn JJk[ x.y.z + t.t.XjU (x + £,y +17, z ,t+T,X 2 ) dijdf 
*-© 


FIGURE 2. REMOTE SENSING OF LOCAL 
TURBULENCE THROUGH CROSS CORRELATION 
OF OPTICAL SIGNALS 


The interpretation of the measured quantity G 
foUows by relating it to the local fluctuations, k, of 
the extinction coefficient. We have shown that G is 
equal to an area integral over the space time corre- 
lation function f 17] , 


G(£, 7, X 

X 


+00 

= ff k(x, y, z+t , t, Xj) k(x+|, y+17, z, t+r, X 2 ) dTjdf , ( 5 ) 


provided that the foU owing assumptions are met: 
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1. The number of statistically independent at- 
mospheric disturbances is small enough to 
avoid a cancellation of the integrated fluc- 
tuation. 

2. The integrated fluctuation of the extinction 
coefficient is small enough to permit linear- 
ization. 

3. The light source fluctuations and the atmos- 
pheric fluctuations are statistically indepen- 
dent. 

4. The fluctuations of the two light sources are 
uncorrelated. 

5. Either extinction or emission dominates the 
change of the detected radiative power. 

The second assumption is not very restrictive since the 
integrated mean value would still be large. Further- 
more, the local fluctuation, k, does not have to be 
small; only the integral over statistically independent 
fluctuations is assumed to be linearized. In the very 
few cases where these weak restrictions are not met, 
one can possibly choose another wave length, X, where 
the second assumption is met. The other assumptions 
can mostly be met by a suitable experimental arrange- 
ment. Therefore, the result of equation 5 appears to 
be sufficiently general for atmospheric applications. 

Two important conclusions can be derived from 
equation 5: 

1. For a pair of statistically independent light 
sources, the source characteristics do not have to be 
known. The method would work for any source, such 
as ground beacons, the sun, or reflected sunlight. 

2. The random nature of turbulent fluctuations 
assures that integrand drops to zero over distances 
which are comparable to a typical eddy size. This 
explains, qualitatively, why the cross correlation 
partially eliminates the usual optical integration along 
a line of sight. 

According to equation 5, the crossed beam corre- 
lation works whenever the wanted signal is common to 
both lines of sight. The way in which the wanted sig- 
nal is then pulled out of the integrated signal may be 
described best by splitting each of the integrated sig- 
nals in two parts. The first describes the contribu- 
tion of the ’’correlated volume" where the integrand 
of equation 5 does not vanish. The second part de- 
scribes the rest of the integrated signal. Multiplying 
the fluctuations of the two integrated signals then leads 
to four products. Three of these change between 


positive and negative values in a random fashion. If 
we now average these products over time, their mean 
value will go to zero, whereas the deviation from the 
mean will increase with the square root of time. 

The fourth product describes the contribution of 
the correlated volume, which is common to both lines 
of sight. The associated product is always positive, 
and its summation over time should increase linearly 
with time. Therefore, it will be the dominant one, 
if one goes on adding long enough. Even small con- 
tributions at the beam intersection can be pulled out 
provided (i) the integrated signals show detectable 
fluctuations and (2) the combination of light source 
fluctuations , shot noise and instrument noise is not 
orders of magnitudes larger than the root mean square 
value of the integrated signal. 


ID. MEASURABLE TURBULENCE PARAMETERS 


The optical integration is restricted to the wave 
front. Along the wave normal, no integration takes 
place and local information inside a correlated volume 
may be obtained by repeating the crossed beam ex- 
periment for several beam separations along the nor- 
mal as shown in Figure 3. The space separation, £, 
describes the minimum beam distance, which defines 
the wave normal. Also a time separation is intro- 
duced electronically by a time delay unit. As will be 
discussed later, the measured signal, G, then ap- 
proximates a two-point product mean value, which 
would have been measured by point probes on the 
wave normal. 


qu. t, x) r _ yt) j y (t+r) <y t=0 

G(0, 0, \) T <yt) i y (t )) ?=0 (I y )^ 0 


M k(x, y, z, t, \) k(x+£, y, z, t+T , X (6) 

k(x, y, z, t, X) k(x, y, z, t, X) 

The position of these imaginary point probes is 
indicated on Figure 3 by two dots. Thus, all turbu- 
lent properties, which are commonly derived from 
two -point product mean values, may also be approx- 
imated with the crossed beam method. 

The local integral scale of turbulence is defined 
by averaging the area under the covariance curve 
G(£, t= 0) and its graphical evaluation is illustrated 
on Figure 3. 
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FIGURE 3. TURBULENCE SCALES AND 
INTENSITIES FROM CROSSED BEAM 
COVARIANCE MEASUREMENTS 

Plotting the beam correlation coefficient against 

the time separation instead of the space separation 

(Fig. 4) allows bulk convection speeds to be read. 

For a given space separation, the correlation will 

reach a maximum at the time that the turbulent eddies 

need to travel from the upstream to the downstream 

beam. The time lag, r__, which is indicated by the 
M 

maximum, is therefore a direct measure of convec- 
tion speeds. 

Power spectra and cross power spectra are de- 
fined by a Fourier transform of the measured auto- 
correlation, 4 = 0 , and cross correlation 4 / 0 , 
curves. They may be calculated using the digital 
computer programs , developed for the analysis of 
mechanical vibrations f 18] . 


A study of the envelope to the time correlations, 
r(T) provides a measure of the temporal rate of 

decay of turbulent eddies as seen by an observer mov- 
ing with the convection speed. Therefore, the decay 
of this correlation curve to 1 /e of its original value 
is used to define the average eddy lifetime. The 
Fourier transform of this envelope would indicate the 
dominant frequencies which are felt by the moving ob- 
server. 


Summarizing the graphical evaluations sketched 
in Figures 3 and 4 ; we find that the crossed beam 
method remotely senses integral scales of turbulence, 
convection speeds, eddy lifetimes and spectra for ar- 
bitrarily chosen atmosphere -fixed points at which a 
beam intersection can be achieved. 



Time Deloy r 


FIGURE 4. LOCAL POWER SPECTRA 
CONVECTION SPEEDS AND EDDY LIFETIMES 
FROM CROSSED BEAM CORRELATION 
MEASUREMENTS 
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The above turbulent properties are only approxi- 
mations , since the optical integration over the wave 
front does catch fluctuations which are not felt by the 
point probe on the wave normal. However, our ex- 
periments in a subsonic jet have shown that weighting 
of contributions close to the wave normal is suffi- 
ciently strong to yield a good spatial resolution. One 
may thus expect a good spatial resolution on a wide 
range of other turbulent motions which are of prac- 
tical interest. 

The absorption (or emission) spectrum of ex- 
tinction coefficients is commonly used to find species 
concentrations and temperatures of both stagnant gases 
and small scattering particles. Adjusting the optical 
wave length of the crossed beams to chosen atomic 
lines and/or molecular bands, one might therefore 
extend the spectroscopic methods from stagnant media 
to inhomogeneous turbulent flows and study multi - 
component or two-phase flow phenomena. 


IV. EXPERIMENTA L VERIFICA TIONS 


The most revealing and instructive test of the 
crossed beam concept is to compare the optical ap- 
proximation of point measurements with known hot- 
wire measurements. The initial portion of an axi- 
symmetric free shear layer was chosen for the ex- 
periments since (1) the turbulence is anisotropic and 
inhomogeneous and (2) convection speeds, turbulence 
scales and eddy lifetimes are well documented in the 
literature. 

All measurements were taken in a subsonic 
(M = 0.2) air jet exhausting through a one -inch 
diameter nozzle into the atmosphere (Fig. 5) . Light 
extinction was achieved by spraying a small amount 
of liquid nitrogen into the settling chamber. This 
produced small water droplets in the exhaust flow 
which attenuate the crossed beams by scattering. 
Later, the nitrogen injection was replaced with water 
spray which offered a better control of the tracer con- 
centration. 

A diagram of the optical and electronic hardware 
is also shown in Figure 5. A plane mirror and a 
spherical mirror collect the light from a powerful 
mercury arc and transmit it as a parallel beam 
through the jet. A similar mirror combination pro- 
jects the arc image on the entrance slit of a grating 
spectrograph. The aperture stop in front of the spec- 
trograph is set at a beam diameter of 2 mm. The 
photomultiplier then scans the first order spectrum 
of the grating. In this way, the wave length interval 


A A. and the wave length X are adjusted by the width 
of the monochromator slit and the rotation of the 
grating, respectively. By exchanging light sources, 
gratings, and photodetectors, the system is able to 
cover the spectrum from the vacuum ultraviolet 
(X = 1200A) to the infrared (X = 25 n ) . 


Jarrel Ash 
UV 

Monochromator 




FIGURE 5. BREAD BOARD MODEL OF CROSSED 
BEAM CORRELATOR 

The above light source and detector arrangement 
is used to generate two beams , one in a horizontal and 
one in a vertical direction (Fig. 5) . A lathe bed pro- 
vides a mechanical support which allows us to align 
the beams parallel to the nozzle exit plane. A beam 
separation in the streamwise direction gives the 
longitudinal turbulence scales and convection speeds. 
Transport of the whole system allows us to repeat the 
measurements for all axial and radial traverses of 
the jet. 

The contribution of correlated light source fluc- 
tuations, created by ripple in the power supply, was 
eliminated using electronic filters. A typical sample 
of optically measured two -point product mean values 
is shown in Figure 6. The beam intersection was ad- 
justed to the center of the free shear layer three exit 
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diameters downstream of the nozzle exit. The tur- regardless of whether nitrogen or water spray was 

bulent convection process is clearly indicated by the used as a tracei . 

displaced maxima for finite beam separation. The 

monotonic decay of the envelope shows the decay of 

the moving eddies. Measuring the time delay for 

which this curve falls to 1/e of its initial value yields 

a value of 880 microseconds. This is in excellent 

agreement with the hot-wire data of [ 19] , from which 

we estimate that the value of the time scale for our 

experimental conditions should be 920 microseconds. 



FIGURE 6 . CROSS -CORRELATION WITH 
DOWNSTREAM BEAM SEPARATION 

Figure 7 shows the same experiment under dif- 
ferent conditions. The initial intersection point of the 
beams was six diameters from the jet exit, (i. e. , 
at x/D = 6.0) instead of at x/D = 3.0. Instead of 
using liquid nitrogen to generate a water droplet fog 
tracer, a small nozzle was mounted in the jet settling 
tank which sprayed water droplets into the flow to 
form the tracer fog. The curves do appear to fall 
onto or close to an envelope as required. Unfortu- 
nately, the range of beam separations used is not 
sufficient to merit the extrapolation of the curve to 
1/e of its initial value. Therefore, the moving axis 
autocorrelation curve has been obtained from the 
hot-wire data of [ 19] and is shown as the dotted line 
superimposed on the measured cross correlation 
curve. 


The spatial resolution of the crossed beam 
method seems very encouraging, considering that the 
envelope from the hot-wire data exactly matches the 
individual curves from the crossed beam method. 

For the moving frame, the optically measured auto- 
correlation therefore agrees with the hot-wire curves, 




FIGURE 7. CROSS CORRELATION WITH 
DOWNSTREAM BEAM SEPARATION 


A quantitative comparison of individual cross 
correlation curves G(r)^ has been based on the indi- 
cated convection velocities. These velocities have 
been evaluated graphically following the procedure of 
Figure 4. To increase the accuracy, we have always 
used more than one curve. The chosen set of beam 
separations was plotted as a function of the time lag 

values t_ , which were indicated by the common tan- 
M 


gent of envelope and individual cross -correlation 
curve. The result is a straight line which is shown 
in the upper part of Figure 7. The desired convection 
speed follows from the slope of this line. 


The results of all convection velocity determina- 
tions made using the crossed beam correlation system 
to date are summarized on Figure 8, where they are 
also compared with hot-wire data from Reference 19. In 
the outer portion of the shear layer ( tj > 0) the convection 
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velocity determinations agree within experimental 
accuracy. For rj < 0, however, the crossed beam re- 
sults tend to follow the mean velocity profile (U/U ) 

o 

rather than attain the constant maximum value ob- 
served for hot-wire data. This was not an unexpected 
result. All convection velocity results presented in 
Reference 19 were taken at axial stations where the 
potential core still exists. At these positions, the large 
difference between the mean velocity and convection ve- 
locity has not been completely explained. However, 
Davies [20] has discussed the subject recently, and 
it does appear that the closer correspondence between 
mean and convection velocity observed here is to be 
expected when a measuring technique other than the 
hot-wire anemometer is being used, and/or the po- 
tential core no longer exists. 




FIGURE 8. VELOCITY PROFILES FROM HOT 
WIRE AND CROSS BEAM METHODS 


ground or on airplanes has the advantages of a larger 
field of view, quicker scanning, and avoiding vibra - 
tion problems. In addition to these benefits, which 
were discussed during the recent remote sensing 
symposium, the satellite detector can use reflected 
sunlight, whereas the ground detector is restricted 
to scattered sunlight. The orbital detector can thus 
use absorption spectroscopy for a detailed study of 
atmospheric and ionospheric constituents, whereas 
the ground detector is restricted to refractive index 
fluctuations, the interpretation of which is very dif- 
ficult. 

Sometimes the fluctuations of atmospheric con- 
stituents may be detectable with the radiative power , 
which is provided by reflected and scattered sunlight 
or by emission inside the detector’s field of view. 

In such a case, the optical hardware which would be 
needed for the crossed beam method is quite similar 
to that which is developed for the photographic sur- 
veys of the earth and the moon. Both applications 
require a diffraction-limited telescope coupled with 
an optical filter which selects the ultraviolet, visible, 
and infrared regions of the spectrum. The main dif- 
ference consists in replacing the camera with a sen- 
sitive phototube and to fly two satellites instead of 
one. 


Two orbital detectors could adjust their lines of 
sight to intersect on an atmosphere -fixed spot, by 
continuously turning the telescope during the fly-by. 
Remote sensing of disturbances around the inter- 
section might be feasible for both the lower and the 
upper atmosphere. Typical lines of sight are sketched 
on Figure 9. 

To study the feasibility of remote sensing in the 
lower atmosphere with respect to the following AAP 
experiments, we suggest: 


The optically measured results in Figure 8 were 
only 0.2 inch apart. This is much smaller than the 
geometrical shear layer thickness (~.8 T1 ) which is 
indicated by the hot-wire results. Apparently, con- 
vection speeds can be resolved inside a correlated 
volume by moving the intersection point. 

V. REMOTE SENSING FROM LOW AND HIGH 
EARTH ORBITS 

The concept of remote sensing through cross 
correlation of optical signals is not restricted by 
scale and the method should, at least in principle, 
work also over atmospheric distances. Mounting the 
two photodetectors on satellites instead of on the 


1. The generation and motion of the aurora and 
airglow using emitted radiation during the 
night cycle. 

2. The motion or generation of water vapor, 
pollutants, and trace constituents using re- 
flected sunlight ( Fig. 9) or the thermal 
radiation. 

Remote sensing in the upper atmosphere is intrigu- 
ing, if the use of ground beacons can be avoided. 

The optical path between the upper atmosphere and 
the orbital telescope is then short. Also, oxygen has 
extremely high absorption peaks in the ultraviolet, 
such that the motion of very diluted oxygen might still 
be detectable. To study the feasibility of remote 
sensing the upper atmosphere and ionosphere with re- 
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RADIATION BELT 


/ 



(C) GROUND AND SATELLITE BASED 




SUN 


EARTH 



• THE METHOD CAN USE THE WHOLE ELECTR0MA6NETIC SPECTRUM FROM 
X RAYS TO MICROWAVES 


• NATURAL RADIATION AND SUN LI6HT MAY REPLACE MICROWAVE AND 
OPTICAL BEACONS 


• A WIDE CLASS OF TEST ARRANGEMENTS IS PROVIDED BY COMBINING GROUND 
AND 0RBITIN6 STATIONS 


FIGURE 9. HOW CAN THE CROSSED BEAM METHOD BE APPLIED IN THE ATMOSPHERE? 
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spect to the following AAP experiments, we suggest: a response time of 0. 1 second along a line of 

sight which traverses the entire atmosphere? 


1. The motion or generation of ozone using di- 
rect or scattered sunlight (Fig. 9b). 

2. The motion or generation of radiation belts 
using the emitted electromagnetic radiation. 
Both vertical and horizontal motions might 
be studied using the same arrangement 
sketched for the upper atmosphere (Fig. 9b) . 


VI. PROBLEMS 


The main problems with the crossed beam 
method are as follows: 

1. To find an atmosphere constituent that pro- 
duces detectable fluctuations of the inte- 
grated signal along accessible lines of sight. 

2. To reduce light source fluctuations, shot 
noise and instrument noise to root mean 
square levels, which have the same order of 
magnitude as the integrated signal. 

3. To have the lines of sight intersect on an 
at mo sphere -fixed spot for an observation 
time, which is sufficiently long to pull the 
local signal out of the integrated signal and 
the combined light source noise, shot noise 
and instrument noise. 


2. Do the usual constituents of the upper atmos- 
phere and ionosphere such as ozone produce 
detectable fluctuations along a line of sight 
which traverses the outer portions of the up- 
per atmosphere? 

3. Do reflected or scattered sunlight and/or 
radiation belts put enough radiative power in- 
to the orbital detector's field of view to re- 
duce the shot noise below the expected signal 
levels, considering the telescopes and optical 
filters which are developed for other AAP 
missions? 

4. Does the continuous motion of the detector's 
field of view lead to fluctuations of the re- 
flected sunlight that are comparable to the 
signal fluctuations? 

5. For which observation time is it possible to 
have the narrow fields of view intersect at a 
point on the lower and the higher atmosphere 
and what are the beam diameters and pointing 
accuracies? 

6. Is it possible to extend these applications to 
the complete day and night cycle by using 
ground beacons besides reflected and scat- 
tered sunlight? 


These problems are formidable, and they may pro- 
hibit some of the applications that were outlined in 
the previous section. However, since we can pick 
the detected radiation from any region of the 
electromagnetic spectrum and since the orientation of 
the lines of sight may be chosen arbitrarily, we be- 
lieve that there is at least one interesting AAP experi- 
ment where the above problems may be overcome 
through special test arrangements. Such an experi- 
ment presents an exciting challenge, since the dynam- 
ic behavior of the atmosphere , or radiation belts 
could be added to the ’’static” mapping on oceanic, 
geographic, and atmospheric surveys. 

Several facts about atmospheric constituents and 
optical AAP hardware must be collected before we 
could decide the feasibility of particular experiments. 
These facts have been identified in the previous dis- 
cussions and may be summarized by the following 
questions: 


REFERENCES 


1. Fisher, M. J. , and Krause, F.: Local Measure- 
ments in Turbulent Flows through Cross - 
Correlation of Turbulent Fluctuations. NASA 
TM X-53295, 1 April 1965. 

2. Fisher, M. J. , and Black, T. H. : Optical Meas- 
urements with High Temporal and Spatial Reso- 
lution. IITRI Progress Report M6 114-18, Con- 
tract NAS8-11258, December 1965. 

3. Nestor, O. H. , and Olsen, H. N. : SIAM Rev. 
2,200, 1960. 

4. Freeman, M. P. and Katz, S. : J. Opt. Soc. Am. 
50,826, 1960. 

5. Bockasten, K. : J. Opt. Soc. Am. 51,943, 1961. 


1. Which constituents of the lower atmosphere, 6. Barr, W. L. : J. Opt. Soc. Am. 52,885, 1962. 
such as clouds, pollutants and trace consti- 
tuents, produce detectable fluctuations within 7. Herlitz, S. I.: Arkiv Vysik 23, 571, 1963. 


44 



REFERENCES (Concluded) 


8. Paquette, D. R. and Wiese, W. L. : Appl. Opt. 
3, 291, 1964. 

9. Howard, J. N. , Burch, D. E. , and Williams, 
D.: J. Opt. Soc. Am. 46, 186, 237, 242, 334, 
452, 1956. 

10. Burch, D. E. , and Williams, D. : Appl. Opt. 1, 
473, 1962. 

11. Burch, D. E. , Gryvnak, D. , Singleton, E. B. , 
France, W. L. , and Williams, D. : Air Force 
Cambridge Research Labs. , Report AFCRL-62- 
698, 1962. 

12. Edwards, D. K. : J. Opt. Soc. Am. 50, 617, 
1960. 

13. Plass, G. in.: J. Opt. Soc. Am. 48, 690, 1958. 

14. Plass, G. N. : J. Opt. Soc. Am. 50, 868, 1960. 


15. Plass, G. N. : J. Opt. Soc. Am. 49, 821, 1959. 

16. Tourin, H. and Krakow, B. : Appl. Opt. 4, 

237, 1965. 

17. Krause, F. R. , and Fisher, M. J. : Optical 
Integration over Correlation Areas in Turbulent 
Flows. Fifth International Congress on Acous- 
tics, Paper K65, Liege, Belgium, September 7- 
14, 1965. 

18. The Vibration Manual . NASA, Geo. C. Marshall 
Space Flight Center, 1st ed. , 1964. 

19. Davies, P. O. A. L. , Fisher, M. J. , and Barrat, 
M. T. , J. Fluid Mech. , Vol. 15, pp. 337-367. 

20. Davies, P. O. A. L. , Near -Field Generation 
by Intense Turbulence. 5th International Con- 
gress of Acoustics, Liege, Belgium, paper L53. 


45 



APPLICATION OF THE Q-BALL ANGLE OF ATTACK TRANSDUCER TO LARGE 

SPACE VEHICLES 

By 

Paul E. Ramsey 


SUMMARY 


LIST OF SYMBOLS 


The angle of attack of large space vehicles can 
be measured either indirectly or directly. Indirect 
methods consist of using accelerometers or obtaining 
rawinsonde wind data, which are used in conjunction 
with telemetered attitude angles and control param- 
eters to calculate the angle of attack. Direct meth- 
ods, which tend to be the most accurate, consist of 
using externally mounted angle of attack sensors. 
Early direct methods consisted of ’’weather cock” 
stable vanes of various types and gimbaled bodies, 
with and without wings , which were boom -mounted 
on the nose of the vehicle. Later, stable vane con- 
figurations and slotted, self-servoed, cylindrical 
probes were extended perpendicularly from the ve- 
hicle skin. These devices were limited in accuracy, 
sensitivity, and use because of their electromechan- 
ical nature, fragility, external mounting positions, 
and mode of operation. 

The Q-ball was developed to eliminate or reduce 
the severity of these limitations normally associated 
with the earlier angle of attack meters. Since the 
Q-ball is faired into the nose of the vehicle, the 
strict requirements for surface smoothness ahead of 
the sensor and local upwash calibrations do not exist. 
Also, since the Q-ball is not a fragile, protruding 
device, it is not likely to sustain damage during ve- 
hicle ground handling. The accuracy and sensitivity 
of the Q-ball are enhanced because its principle of 
operation requires no mechanical moving parts; 
rather it uses sensitive pressure transducers and 
associated tubing and electronics. 

The Q-ball provides a valuable method for de- 
termining the angle of attack of any space vehicle 
during most of the atmospheric portion of flight. 
Because of its accuracy and ability to detect slow 
angular divergence rates away from normal flight 
angle of attack, the Q-ball has been used to obtain 
data for post-flight analysis and is being developed 
for use in the manned Apollo-Saturn emergency de- 
tection system to detect abort situations. 


Symbol 

M 

q 

p 

a 


a 


4 > 

AP 

c 

AP 

AP 


Definition 

Free stream Mach number 

Free stream dynamic pressure, psia 

Pressure, psia 

Angle of attack, degrees 

Angle measured from stagnation point 
on the spherical nose, degrees 

Angle of yaw, degrees 

Pressure differential measured across 
the pair of 45 ° ports in the pitch plane , 
psia 

Pressure differential measured across 
the pair of 45° ports in the yaw plane, 
psia 

Pressure differential measured across 
the Q -ports, psia 


Q AP 

— (— )f 

da K q f \ 


AP 


or 

AP 


Slope of 


a 


versus a curves at 


a 


a q 
K 


a = 0 °; 


AP a 

_9 


deg. 


AP 


a 


max 


Pressure coefficient at stagnation point 


Subscripts 


Measured at angle 0 
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LIST OF SYMBOLS (Concluded) 
Symbol Definition 

°o Free stream conditions 

1 Upper a port 

2 Lower a port 

39 Q~port located 39° from Q-ball 

center line 

51 Q-port located 51° from Q-ball 

center line 


L INTRODUCTION 


The control systems of today's advanced space 
vehicles are a relatively complex organization of ve- 
hicle attitude sensors, electronic computers, and 
actuators needed for movement of a portion of the 
vehicle, such as a fin, an engine, or a control sur- 
face. The major function of such a system is to 
keep the space vehicle on a prescribed path as closely 
as possible. An equally important function for the 
control system, especially for large slender Saturn 
vehicles, is to keep lateral aerodynamic loads within 
the failure limits of the vehicled structure. Since 
these loads are proportional to the angle of attack 
and the flight dynamic pressure, the allowable angle 
of attack must be restricted, especially during the 
high dy na mic pressure regions of flight. Also, since 
the Saturn vehicles are generally aerodynamic ally 
unstable throughout most of the flight range, the con- 
trol system must induce artificial stability. These 
control system functions are normally accomplished 
through some means of engine thrust vector control. 

To determine the amount of thrust vectoring 
needed, the vehicle attitude and angle of attack or 
accelerations caused by angle of attack induced forces 
must be known at all times during flight. An excep- 
tion would be the attitude control system presently 
planned for Saturn V which will not require angle of 
attack inputs. This information is fed into a control 
computer which determines the required engine de- 
flection according to predetermined gain factors. 

The attitude or vehicle orientation with respect to 
the vertical can be determined by measuring the 
angle between the axis of a space -fixed gyro- 
stabilized platform and the vehicle centerline. Angle 
of attack is normally measured during flight, either 
directly or indirectly, by any one of various means. 


Direct methods, which tend to be the most accu- 
rate, use externally mounted angle of attack meters; 
indirect methods use accelerometers or rawinsonde 
wind data along with telemetered attitude angles and 
control parameters to calculate the angle of attack. 

In the past, several types of angle of attack meters 
have been used. Early devices consisted of ’'weather 
cock" stable configurations which were mounted on 
the nose of the vehicle and aligned themselves with 
the direction of the winds. Boom-mounted vanes and 
gimbaled bodies with and without wings were the most 
common examples. Later devices consisted of stable 
vane configurations and slotted, self-servoed , cyl- 
indrical probes which were extended perpendicularly 
from the vehicle skin into the flow field around the 
vehicle. 

All of these devices had limitations characteris- 
tic of their mode of operation. Some of these limi- 
tations are discussed below: 

a. Local body -mounted probe and vane sen- 
sors require an accurate calibration of the upwash 
around the vehicle in the vicinity of the angle of attack 
sensor for each sensor ^vehicle combination. 

b. The surface ahead of body -mounted sen- 
sors must be free of protuberances or sudden discon- 
tinuities which could cause unpredictable local flow 
conditions. 

c. Because local body -mounted sensors are 
electromechanical devices, good sensitivity and ac- 
curacy are limited to relatively high dynamic pres- 
sure regions of flight where the aerodynamic forces 
are sufficient to overcome bearing friction and the 
momentum of mechanical parts. 

d. Because these angle of attack devices are 
usually fragile and are mounted externally , they can 
be easily damaged during ground handling, and are 
subjected to aerodynamic heating during the ascent 
and reentry phases of flight. 

To eliminate or reduce the severity of these 
limitations , a more advanced sensor for use on large 
space vehicles was developed by adapting the movable 
nose sphere Q-ball used on the X-15 research air- 
craft. The result was a fixed-sphere Q-ball which 
indicated the angle of attack by sensing differential 
pressures on the nose sphere. The Q-ball was de- 
signed to be incorporated into the nose cones of 
ballistic missiles (without disturbing the overall 
vehicle aerodynamic characteristics) and to be 
sturdy enough to withstand aerodynamic heating dur- 
ing the ascent and to remain essentially intact, even 
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though inoperative, when exposed to the heat of at- 
mospheric reentry. 

During the development period of the Q-ball, 
angle of attack devices gave way to accelerometers 
for control purposes, but were and still are required 
for other applications. Because of an accuracy ad- 
vantage over accelerometers, angle of attack meters 
(such as the Q-ball) are used to measure precisely 
the angle of attack during flight. This information is 
telemetered back to Earth for postflight analysis of 
vehicle performance. A Q-ball is now being devel- 
oped for use as part of an abort warning system for 
the manned Apollo -Saturn flights. In this capacity, 
the Q-ball can detect slow divergence from normal 
flight angles of attack toward the abort limit more 
accurately than accelerometers. 

This paper is concerned with the Q-ball and its 
application to large space vehicles, and with the 
major problems and sensor accuracies encountered 
during actual use. 

n. DESCRIPTION AND OPERATION OF Q-BALL 

The standard O-ball is a sphericallv blunted 15° 
half-angle cone with diametrically opposite differen- 
tial pressure orifices which are located 45° from the 
stagnation point in both the pitch and yaw planes . A 
pair of differential pressure orifices, called Q-ports, 
are located 39° and 51° from the stagnation point in 
a plane rolled 45* from the pitch plane. Figure 1 
shows the Q-ball configuration and location of the 
three pairs of pressure ports. 


FIGURE 1. PORT CONFIGURATION OF 
Q-BALL 


The Q-ball works on the principle that a linear 
function exists between the differential pressures 
measured across the 45° (a and ip) ports and the 
angles of pitch and yaw. Because the Q-ball is axi- 
symmetrical, this pressure -to -angle relationship 
can be represented by a single set of calibration data 
for both the pitch and yaw planes. Therefore, for 
simplicity, all further comments will be restricted 
to the pitch plane. Calibration data are obtained ex- 
perimentally by placing a scale model of the Q-ball 
into a wind tunnel (Fig. 2) and measuring the result- 
ing AP as a function of a. Since the differential 

pressure variation with a is linear for angles of at- 
tack up to approximately 15° , the data are conven- 
iently presented as the slope of the AP^/q versus a 

curves or AP /aq as a function of free stream Mach 
a 

number. Typical data are presented in Figures 3 and 
4 for several Q-ball sizes and configurations. More 
complete data are shown by Mills and Lukesh f 1] ana 
Ramsey [2], Also, for comparison, suppersonic ex- 
perimental data for hemisphere -cylinders, along with 
some subsonic perturbation theory [3] and modified 
Newtonian Theory, are presented in Figure 4. 



FIGURE 2. TYPICAL LARGE SCALE Q-BALL 
EXPERIMENTAL MODEL SHOWN INSTALLED 
IN THE CORNELL AERONAUTICAL 
LABORATORY 8-FOOT TRANSONIC 
WIND TUNNEL 
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Sym Moch 
O 0.40 
A 0.80 
k 1.00 
d 1.25 



Angle of Attack a 

FIGURE 3. VARIATION OF DIFFERENTIAL PRESSURE ACROSS THE 45° PORTS/q 

WITH ANGLE OF ATTACK 




Sym 


Facility 


Model 

V 

AEDC 

Tunnel A 

(Ref. 5) 



A 

MSFC 

14 x 14 inch 

TWT (Ref. 6) 


o 

Cornell 

8' TWT, 2“ 

nose R., 

13" Dia. 

(Ref. 2) 

□ 

Cornell 

8' TWT, l" 

nose R., 

13" Dia. 

(Ref. 2) 

— 

Ames 

Research Center, F~I6 

Q-Ball 

(Ref. 1) 


Subsonic Perturbation Theory (Ref. 4) 
Modified Newtonian Theory 



Cornell 8 TWT, I nose R., 6.5 Dia 


Mach Number , M 


FIGURE 4. SLOPE OF DIFFERENTIAL PRESSURE ACROSS a -PORTS/q VERSUS a CURVES 
AS A FUNCTION OF FREE STREAM MACH NUMBER 



The modified Newtonian data were obtained by 
using the familiar Newtonian expression to write the 
equation for the pressure at each of the a -ports and 

a Ap 

combining them to obtain 7 ( ) for the pitch 

3 a q 

plane. Mills andLukesh [1] give a complete deriva- 
tion which includes the effects of yaw. 


The pressure on a sphere at any point 0° from 
the stagnation point is as follows: 

P = P + C q cos 2 0 . 

0 <*> p 

* m o -v 


Subtraction of Equations (1) and (2) and substituting 
the expressions for 0 A and 0 2 yields 



C cos 2 - C cos 2 0 2 

P P 

max max 


cos 2 - a) -cos 2 + o; ) 

Now, we have 

cos (t - oi) = cos 7 cos a + sin 7 sin a 
v 4 4 4 


AP 


a 


~ C 


max 


(3) 


At some angle of attack (or) , the 0 angles for ports and 

1 and 2 , shown in Figure 5, are as follows: 

7T . 7 r 

cos 7 = sin 7 . 
4 4 



Substituting we get 

cos ( j - O') = sin ^ ( cos 01 + sin 01 ) . 

Therefore, 

cos 2 (j - a) = sin 2 - (cos 2 a + 2 sin a cos a + sin 2 a) 

cos 2 ( “■ + a) = sin 2 - (cos 2 a - 2 sin a cos a + sin 2 a). 
Substitution in Equation (3) yields 
AP 

= C sin 2 7 (4 sin a cos o?) 

q P 4 


FIGURE 5. NOTATION FOR PRESSURE 
CALCULATIONS WITH Q-BALL AT AN 
ANGLE OF ATTACK 


7T 

01 = - - O 

*2 = f + « • 

Therefore, 

P t = P^ + C q cos 2 0! (1) 

P max 

p 2 = P oo + ^ Q cos 2 0 2 . (2) 

P max 


AP 


= C 


sin 2 a , 


max 


Differentiating, we obtain 


AP 


a a 


■) = 2 C 


cos 2 a . 


max 


8 


AP 


a 


(4) 


) in radian 


Equation (4), which determines Q a (' 

measure ( — ,) , was converted to degrees before being 
rad 

presented in Figure 4. 

The Q-ball measured differential pressures are 
also a function of the dynamic pressure level; conse- 
quently, a means of determining the dynamic pressure 
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(5) 


is necessary, which is the function of the Q-ports. 

. Differential pressures measured across these ports 
are directly related to free stream dynamic pressure. 
Necessary calibration data are normally presented 
as the differential pressure across the Q-ports, AP , 

Q 

ratioed to free stream dynamic pressure, q, or 
simply AP^/q as a function of Mach number. Typi- 
cal AP calibration data are presented in Figure 6 

q 

for several Q-ball configurations . Some hemisphere- 
cylinder data, along with subsonic perturbation and 
Newtonian theories, are also presented for compari- 
son with the AP calibration data. 

q 


K = 


(AP /q) 
-SL- 


AP a 
JL 


AP /oiq A P_ 
rv (X 


Then, to obtain the vehicle angle of pitch, the flight - 
measured AP^ and AP^ are ratioed and multiplied 

by K as follows 2 


AP 


ol = K ( 


a 


AP 


) 


( 6 ) 


T 

measured in flight 


111. Q-tsALL TYPES 


The AP /q and AP / ctq data can be combined to 
q a 

obtain the flight angle of pitch, a , in the following 
manner. First, for any Mach number, AP / aq and 


the corresponding AP /q value can be used to de- 

q 

termine a K factor (represented by \ [4]) as follows: 


There are three versions of the Q-ball either in 
existence or under development. The Q-ball princi- 
ple was used on early Saturn vehicles by measuring 
differential pressures across pairs of holes drilled 
into the nose cone. A complete Q-ball package known 
as the F-16 was flown on Saturn SA-4 through SA-7 
vehicles. This is a normal Q-ball which has a 2-inch 
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0.32 
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0.16 
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AEDC Tunnel A (Ref. 5) 

MSFC !4 x 14- inch TWT (Ref. 6) 

Cornell 8* TWT, 2" nose R.,13" Did. (Ref. 2) 
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Saturn S A - 6 , Flight data 
Modified Newtonian Theory 
Subsonic Perturbation Theory (Ref. 4) 

Ames Research Center (Ref. I) 

Preferred fairing for Cornell 13" Dia. (Ref. 2.) 



Cornell 8' TWT , l" nose R., 6.5" Dia. 
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FIGURE 6. THE DIFFERENTIAL PRESSURE ACROSS THE Q-PORTS/q AS A FUNCTION OF 
FREE STREAM MACH NUMBER FOR ZERO ANGLE OF ATTACK 
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nose radius and was used with the q -measurement as 
explained above. 

The measuring Q-ball, a second highly refined 
version now being developed at MSFC , has a nose 
radius of 4 inches and operates on the same principle 
as the F-16 model, but uses six pressure transducers 
instead of three. This arrangement provides a low 
and high range transducer for each of the three dif- 
ferential pressure measurements (AP , AP , and 

a $ 

AP ) . The extra set of transducers extends the 

q 

range of this Q-ball to cover a larger portion of the 
vehicle flight period. Flown on Saturn I vehicles 
SA-8, SA-9, and SA-iO, this version will be used on 
future vehicles requiring angle of attack measure- 
ments. 


The third Q-ball, which is presently of greatest 
interest, is the emergency detection system (EDS) 
Q-ball being developed for the manned Apollo-Saturn 
flights. This device has a 4 -inch nose radius and 
only Of and 0 ports. Since the differential pressures 
are a function of both the vehicle angular attitude and 
dynamic pressure, they are an effective indication of 
the aerodynamic loads on the vehicle. During flight, 
the vector sum of AP and AP , will be displayed to 

01 if) 

the astronaut. Should this sum go beyond a set limit 
or ’’red line” and if the astronaut has confirmation 
from other EDS sources, he may initiate manual 
abort. The ’’red line” will be based on vehicle struc- 
tural limits, trajectory data, and Q-ball calibration 
data. 


IV. CALIBRATION PROBLEMS AND Q-BALL 
LIMITATIONS 


The existing experimental evidence indicates that 
the AP^ measurements are relatively unaffected by 

afterbody configuration, scaling effects, or slight 
hardware inaccuracies due to manufacturing toler- 
ances. On the other hand, because of the disagree- 
ment between data obtained from various experimen- 
tal test models and flight Q-balls, the AP measure - 

q 

ment appears to be very sensitive to these same 
parameters. This is indicated in Figure 6. A natu- 
ral consequence is that any AP measurements ob- 

q 

tained experimentally will have a limited accuracy 
because the full scale flight hardware and flight aero- 
dynamic parameters cannot be exactly simulated. 


During the present discussion, the AP /q cali- 

q 

bration curves have been conveniently assumed to be 
constant with variations in a , but an examination of 
the typical AP /q versus a data, obtained from 

q 

Ramsey [2] and presented in Figure 7, indicates that 
there are significant effects which must be considered 
when the Q-ball is used. In effect, this reduces the 
precision of the angle of pitch as computed from 
Equations (5) and (6) f . To obtain an accurate a, a 
second AP^/q value, which corresponds to the first 

computed a , must be obtained from data similar to 
that in Figure 7 and used to calculate a new K factor. 
This K can then be used to compute a new a which is 
used to start the entire process again. This iteration 
should be made at least three times to obtain the 
necessary accuracy. 


The AP^/ q nonlinear variation with a apparently 


is caused by the location of the Q -ports in the plane 
rotated 45° from the pitch plane. By using spherical 
geometry to derive expressions for the angular dis- 
tances of the Q -ports from the stagnation point and 
using the Newtonian expression for the pressure dis- 
tribution, an equation that describes AP /q for vari- 

q 


ous angles of attack can be derived as follows. 


Figure 8 indicates the geometry involved in com- 
puting 6 for the Q -ports. For the 39° port, 


sin b = sin B sin c 


(7) 


tan a 


cos B . i 

“7 and — 

ctn c ctn c 


tan c. 


Therefore, 


tana = cos B tan c. (8) 

Equations (7) and (8) become 

sin b = sin 45° sin 39° (9) 

tan a = cos 45° tan 39° (10) 


Now, since we have a right spherical triangle, 

cos 0 39 = cos (a - oi) cos b. (il) 

Solving Equations (9) and (10) and substituting into 
Equation (11), we obtain 

cos 0 39 = 0. 8956 cos (29. 80 - a) . ( 12) 
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FIGURE 8. THE RELATIONSHIP BETWEEN THE 
Q -PORTS AND THE STAGNATION POINT AT 
ANGLE OF ATTACK 


Similarly , 

cos $ 5 ! = 0.8355 cos (41.13 - a). (13) 


The Newtonian expressions for the two Q -ports 

are 


% 

= P + C q cos 2 0QQ 

oO p 00 

max 

(14) 

P *51 

= P + C q cos 2 05 i . 

oo p 

max 

(15) 


Rearranging and subtracting Equations ( 14) and 
(15), we obtain 


AP 
9. 

q 


c (cos 2 0 39 - cos 2 0 51 ) . 

p 

max 


(16) 


Squaring Equations (12) and (13) and substituting 
into Equation (16) yields 



C 

p 

max 


[cos 2 (29.80 - a) 0.8019 - cos 2 (41. 13 -a) 0.6981]. 

(17) 


Data obtained from Equation ( 17) are presented in 
Figure 7 for angles of attack between ±10° and a 
Mach number of 1.25. 

Because the theory is less valid for low super- 
sonic Mach numbers, the modified Newtonian values 
are lower than the experimental data, but the non- 
linear trends are similar. The locations of the ex- 
ports cause the distance between each of the two ports 
and the stagnation point to vary differently with chang- 
ing a as is shown in Figure 8. Since the pressures 
sensed by these ports are a function of the distance 
from the stagnation point, the measured AP^ 

nonlinear with a. The difference between the two 
cos 2 0 terms in Equation ( 17) represents this con- 
clusion. 

Because of these limitations of the dynamic pres- 
sure measurement, it is sometimes desirable to use 

the Q-ball without the AP measurement and instead 

q 

to rely on the dynamic pressure obtained from other 
sources such as some type of on-board pilot-static 
probe or calculation from the telemetered trajectory 
parameters. The source which is used to obtain the 
dynamic pressure measurement should be determined 
by whether the Q-ball is used as a measurement in- 
strument or as a control input. For instance, if the 
Q-ball were being used as a measurement device, 
immediate on-board dynamic pressure information 
would not be necessary since all data would be tele- 
metered back to Earth for postflight analyses. 

A rigid mounting must be provided for the Q-ball 
to prevent errors caused by flexing of the mounting 
structure. Also, the vehicle configuration should 
provide for a nose -mounting position for the Q-ball to 
ensure that the sensor will be located in relatively 
undisturbed free stream flow. 


V. ESTIMATED Q-BALL ERRORS 


The probable errors for the Q-ball are dependent 
not only on the accuracy of the angular pressure meas- 
urements (AP and AP ), but on whether the dy- 

a ip 

namic pressure measurement is used. The 99.73 
percent probability errors for a normal Q-ball using 
both the angular and dynamic pressure measurements 
are presented in Figure 9 for a typical Saturn tra- 
jectory. These errors would be somewhat less if a 
more refined externally obtained dynamic pressure 
were used. 
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FIGURE 9. THE 99.73 PERCENT PROBABLE ERROR FOR THE Q-BALL DESCRIBED AS A 
FUNCTION OF MACH NUMBER FOR 5* AND 10° ANGLE OF ATTACK 


Other sources which contribute to the errors 
presented in Figure 5 are listed below. 

a. AP^ measurement error due to the pressure 
transducer inaccuracies. 

b. Manufacturing errors, such as port mislo- 
cations for the flight Q-ball and the test models, 
which introduce inaccuracies in the calibration curves. 

c. Slight nonlinearities in the pressure distri- 
bution over the spherical nose. 

d. Q-ball misalignment with respect to the 
space vehicle. 

e. Vehicle subsonic upwash ahead of the Q-ball 
transducer. The upwash cannot be completely cor- 
rected by wind tunnel calibrations because of the 
impracticability of placing the complete Q-ball/ve- 
hicle combination in a wind tunnel. Calculation of 
the upwash is impractical because of the range and 
complexity of possible vehicles on which the Q-ball 
will be flown. 

Since vehicle elastic bending, which causes 
Q-ball angular misalignment, can be different for 
various vehicle configurations, it has not been in- 


! 
| 
i 

eluded. In most cases, elastic bending can more 
than double the expected error in measurement of 
angle of attack. 

Figure 9 indicates that the accuracy of the Q- 
ball drops off rapidly during the early and latter part 
of flight because of the relatively low dynamic pres- 
sure levels in these regions. The Q-ball differential 
pressures drop to such low values that they approach 
the measuring inaccuracies of the pressure trans- 
ducers used in the Q-ball. The accuracy of all ex- 
ternally mounted angle of attack sensors suffers 
during these regions of flight, but because of its 
greater sensitivity, the Q-ball is less affected than 
the other devices. 


VI. COMPARISON WITH ACCURACY OF 
OTHER METHODS 


The relative accuracies of the various methods 
of measuring the angle of attack are investigated by 
Addison [4] . Saturn SA-2 through SA-4 flight data 
from the Q-ball, rawinsonde winds, and acceler- 
ometers were compared to a slotted , cylindrical 
probe sensor, which was used as a reference be- 
cause, at the time of this study, the probe device 
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had been flown successfully on a large number of 
Jupiter missiles and the early Saturns. The probe 
sensors have limitations, such as a need for an ac- 
curate calibration of the upwash field around the ve- 
hicle and possible misalignment with the vehicle 
centerline which can introduce errors. However, if 
these errors are assumed to be reasonably small, 
the various methods for determining angle of attack 
can be evaluated. 

The Q-ball data compared favorably with the 
probe sensor data for most of the vehicles studied ; 
therefore, the Q4mU sensor was considered to be an 
accurate method for measuring angle of attack. 


vn. CONCLUSIONS 


The angle of attack of large space vehicles is 
difficult to measure accurately because of the in- 
evitable limitations of any known technique. Direct 
methods, i.e. , externally mounted angle of attack 
sensors, are inherently more accurate than indirect 
methods such as accelerometers and rawinsonde 
winds. 

The limitations of various externally mounted 
sensors used in the past are summarized as follows: 

a. Local body -mounted probe or vane sensors 
require an accurate calibration of the local upwash 
around the vehicle for each sensor -vehicle combina- 
tion. 

b. The surface ahead of body -mounted sensors 
must be smooth and free of structures and protuber- 
ances which could cause unpredictable local flow 
conditions. 

c. Bearing friction and momentum of mechani- 
cal parts reduce the accuracy, sensitivity, and 
therefore the range of flight times for which the 
sensor is useful. 

d. Externally mounted sensors are easily 
damaged during vehicle ground handling and are sub- 
jected to aerodynamic heating during ascent and re- 
entry. 

The Q-ball was developed to eliminate or reduce 
the severity of the limitations of the earlier angle of 
attack meters. Because the Q-ball is faired into the 
nose ahead of the vehicle itself, the requirements 
for surface smoothness ahead of the sensor and local 
upwash calibrations no longer exist. Also, since 
the Q-ball is not a fragile device which protrudes 


from the vehicle skin, the sensor is not as likely to 
sustain damage during ground handling. Although it * 
does not maintain an operational capability , the Q- 
ball is designed to survive reentry. Accuracy and 
sensitivity of the Q-ball is enhanced because its prin- 
ciple of operation requires no mechanical moving 
parts, only sensitive pressure transducers and asso- 
ciated tubing and electronics. 

It can be concluded that the Q-ball provides a 
valuable method for determining angle of attack of 
any launch vehicle during most of the atmospheric 
portion of flight, especially the important maximum 
dynamic pressure region. Because of its accuracy 
and ability to detect slow angular divergence rates , 
the Q-ball has been used to obtain data for post -flight 
analysis and is being developed for emergency de- 
tection system purposes. 
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IV. ORB IT THEORY AND PREDICTION 



A SOLUTION OF A PROBLEM OF ORBIT DETERMINATION WITH POSSIBLE 
APPLICABILITY TO THE ABORT GUIDANCE PROBLEM 

By 


Robert J. Hill and William D. Goldsby, Jr. 


SUMMARY 


H. ANALYSIS 


This paper presents a closed solution to the prob- 
lem of orbital determination if two position vectors 
and the path angle at one of the vectors are given. The 
orbit thus determined may be considered as having 
application to the abort guidance problem where an al- 
ternate mission might be the alternative to mission 
failure. 

The procedure here developed is available in an 
operative computer program. A slight modification 
in the program could allow the user, by varying the 
specified terminal position vector, to iterate on some 
final desired geocentric coordinates. 

I. INTRODUCTION 


Consider a vehicle with a given position vector 
( P) and velocity vector (V) traveling in an isolated 
central force field. What is the elliptical path which 
the vehicle must follow so that it will arrive at a de- 
sired later position with a prescribed direction of the 
velocity vector (flight path angle, *?)? This is, in 
effect, a statement of the abort guidance problem, the 
application of which is assumed in investigating sit- 
uations in which alternate missions may become 
necessary. 

The problem is solved in closed form by use of 
some geometrical properties of the ellipse and by use 
of some standard procedures of celestial mechanics. 
The approach is first to determine the required path 
(ellipse) the vehicle must follow, and then to deter- 
mine the necessary velocity vector needed to transfer 
the vehicle from its present path to the required el- 
lipse. 


The abort guidance problem is one of insuring that 
certain specified end points of vehicle position and 
flight path angle are met after mission abort. 

For the purpose of this analysis , this problem 
statement implies the determination of the required 
ellipse to meet the specified end point and the velocity 
increment necessary for the vehicle to transfer to that 
ellipse. 

The adoption of the following notation is conven- 
ient for the required analysis: 

Pi - the position vector of the vehicle at the time 
the abort is initiated. 

P 2 - the position vector of the vehicle desired at 
some later time. 

F t - the position of the attracting body. The co- 
ordinates of F j are taken without loss of 
generality as the origin 

F j - (0, 0, 0) T 

I = (1,0,0) T 
J = (0,1,0) T 
K = (0,0, i) T 

The angle, A <p , between P t and P 2 in the plane 
determined by Pi and P 2 may be computed in the 
following manner: 
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cos cr 


N ■ I 



cos A <p = — w 

|P|I IPil 


sin A cp — 


(Pi x Pi) 
I Pi I IP 2 1 


(Pi x P 2 ) 

I Pi X P 2 ! 


sin a = (N x I) 

-i / sin a \ 
V cos cr / 


a = tan 


(N x I) 

|N x l\ 


Thejingle, i, which is the inclination of the 
plane of P t and P 2 to the X- Y plane, is given by 


Pi x P 2 

cos i = — — * K 

IPi x P 2 | 


sin i = 


(P t x P 2 ) x K 
I (Pi X P 2 ) X K | 


N 


A 


(p = tan 


-1 


( sin A <p \ 
cos A <p / 



Now, let F 2 denote the position of the other focus 
and have coordinates (X, Y,Z) to be determined. In 
order to determine these values , it is convenient to 
transform the coordinate system so that the X and Y 
axes lie in the plane of TV and P 2 and so that the three 
coordinate axes form a right-handed system. In or- 
der to perform this transformation, the following 
computations are made. Let £ be the angle meas- 
ured in the plane of P t and P 2 from P t to the inter- 
section of the X-Y plane and the plane of P* and P 2 . 

The unit vector , N, in the direction of the de- 
scending node is given by 

_ (Pi x Pi) XK" 

N = — — - 

l(P t x P 2 ) XK j 

The angle £ is then given by 
P 4 ~N~ 

COS £ = — — 

I Pi ! 


Now the transformation may ue made from the 
(X, Y, Z) system tojhe (X f , Y f , Z') system, where 
the X axis lies jdong Pi and the Z axis is perpendicular 
to the plane of P 2 and P 2 and the (X, Y, Z) system 
forms a right-handed system. 


X’ = ( 0)3 (i)i (£) 3 X 

where 


<cr) 3 


( cos a 
-sin a 
0 


sin a 
cos a 
0 




cos i 
-sin i 


sin 

cos 


and 

U>3 



sin £ 
cos £ 
0 



The transformation takes the vector JP t into the vec- 
tor P t = (Pi, 0, 0) and the vector P 2 into the vec- 
tor P 2 = ( 1 P 2 1 cos A ip , | P 2 1 sin A <p, 0) . 


(Pi x N) ^ (Pi x N) 

IPlI IPi x N| 



The angle , a, which is the angle between the 
line of nodes of the two planes and the X axis , is 
given by 


A useful property of an ellipse is now stated. 

Given the foci, Fi and F 2 , of any ellipse and any 
point P on an ellipse , then Jie angle that F t P makes 
with the tangent at P is equal to the angle that F 2 P 
makes with the tangent at P. This property is es- 
sential to the analysis and is proved in the Appendix. 

The angle that the vector F*i makes with the velo- 
city vector is 180-^j where is the path angle at P*j. 
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Then the angle, °c , between the extension of and 
F 2 is given by 

oc = 360 - 2i» t . 

If the origin is translated to P^ and a rotation 
through the angle oc is performed, the coordinates 
of F m 2 , the double prime denotes the third coordinate 
system, become (X", 0, 0). The transformation 
may be represented in general by 

X” = ( 00 ) 3 (X 1 - P T t ) where 

( cos oc sin oc 0 \ 

-sin oc cos oc 0 1 . 

0 0 1 / 

Another property of the ellipse is shown: 

Given any point P on an ellipse and the two foci, 

then 

|Fj - P | + IF 2 - P I = Constant. 

Using this property yields 

1 

Pj + X" = Pj+J^Y" - Y" 2 ) 2 + (X" - X" 2 ) : J 2 . 

But Y" = 0, Z" = 0, after the transformation. 

Solving for X" results in 

X„ = Y" 2 2 + X" 2 2 - (P, - P 2 ) 2 

2(P, - P 2 + X" 2 ) 


Then 

X' = (a ^ X” + P, 

gives the coordinates of F 2 . Since F x is at^the origin, 
the distance between the foci is given by | X f | . The 


distance between P ! t and F’ 2 is given by | P' 2 - F ! 2 | 
and the semi-major axis, a, is given by 

a = p i + IP f i - ^ 2 [ 

2 


For any ellipse, the distance between the center and 

1 

either focus is given by (a 2 - b 2 ) 2 where b is_the semi- 

I X* | 

minor axis. This distance is also given by since 

F'j is at the origin. 

Hence, (a 2 -b 2 ) 2 = ^ ^ ; 



The numerical eccentricity, e, is given by the rela- 


tion 




(a 2 -b 2 ) 2 

e = — L — (0<e<i). 

a 

The semilatus rectum, p, is given by 
b 2 

p = T • 


The angular momentum, h, is given by 


h = (/x p) 2 , h > 0 , 


where n is the gravitational constant. 

Now, the magnitude of the velocity at P’ 2 is given by 



The unit vector in the direction of peri -center U 1 is 
is given by 



I XM 


The true anomaly, <p 2 , of the position vector P’ 2 
may be obtained as follows: 

U f • P' 2 

cos <P 2 = 3 

IP’ol 


sin cp 2 = 


<p 2 = tan 


-1 


U 1 x P t 2 

|P f 2l 
sin <p 2 
cos <p 2 


K 
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The path angle, , at P T 2 may be computed 



where 90 < — 180 if <p 2 ^ 180 

0 < i> 2 < 90 if <p 2 < 180 . 


interval, the earth would have rotated so that P 2 would 
no longer be at the specified geocentric position. In 
order to use this solution, then, an iteration would be 
necessary to match the time of travel along the de- 
termined ellipse and the rotation of the earth. 


APPENDIX 


The velocity vector X' 2 may be computed by the 
relationship 


X ’ 2 = (-A <?) 3 


V 2 COS t ? 2 

V 2 sint? 2 
0 


Then the velocity vector X 2 in the original coordi- 
nate system is given by 

— T T T — 

X 2 = (£>3 ( i ) 1 (c)3 X f 2 . 

• • • 

Then the vector A V = X 2 - X where X is the 

m m 

velocity of the_vehicle_at the time of abort. The 
magnitude |AV[ = |X 2 - | and the angle £ be- 


sin £ = 


£ = tan 


. 

m 

x 2 

• X 

m 

|X 2 

m 


X X m > 

|X 2 

1 

m 

. L 

sin £ 

\ 

cos £ 


in. 


in. CONCLUSION 



Fl = (c»0) ; F 2 (-C, 0) ;a 2 = b 2 + c 2 

F = “S” + = 1 ;b 2 X 2 + aV = a 2 b 2 . 

Then = 2b 2 X + 2a 2 y ^ = 0. 
dx J dx 

Hence, the slope of the tangent at any point 
P(x,y) is given by 

= s 

dx aV 


This paper has presented a closed solution to a 
problem of orbit determination that has direct appli- 
cation to the abort guidance problem. The problem 
can be presented mathematically as follows: 


This slope is undefined aty = 0 (apogee and perigee) ; 
however, the angles j3 t and 0 2 are readily attainable 
here. 

Now consider any point P (x,y) where X ^ ±c and 
y ^o. 


Given two radius vectors and the path 
angle at one of them; determine the unique 
ellipse which satisfies these conditions. 

If this solution were to be applied to the abort 
guidance problem, it would be necessary to place the 
vector Pj at a specified geocentric position; however, 
the time it takes a vehicle to traverse from P 2 to Pj 
must be taken into account, and during this time 


The slope of R 2 = 
and 

the slope of R t = 


x+c 

c-x 


= m 2 x ^ -c 
= m t x ^ c . 
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-b 2 * y 

S -“l 1?T C-X 

•“O' ■ l* (s) (m,) ■ 


1 + 


bSc y 

aV c -x 


-bSc (c-x) + aV 
y a 2 (c-x) + b 2 * 


-b 2 xc + bV + aV 2 _ -b^c + a 2 b 2 
y a 2 c - x(a 2 -b 2 ) yc (a 2 -xc) 

b 2 (a 2 -xc) = b^_ 

= 3 1 

yc (a -xc) yc 


Hence , tan. fa = (y ^ 0 , X ^ c) 


m 2 - s 


b 2 x 


tan 0 2 = 


x+c 


1 + (s)(m 2 ) 


+ ?y^_ 


x ^ -c 
y ^ o 


y b 2 x 
x+c aV 


aV + b^x (x+c) 
y a 2 (x+c) V bht 

aV 2 + bSc 2 + bStc 

y x(a 2 -b 2 ) + a 2 c 

,2 U 2 


a 2 b 2 + b 2 xc _ bf_ (a* +xc) = b_ 
y xc 2 + a 2 c yc (xc + a 2 ) yc 


yc (xc + a ) yc 


Hence, tan fa = ~ = tan 0j 

since for y f 0, x ± ± c; 0i, fa f- ± 90 
fa = fa ~ ^ • 


In particular, let us take the case where x- c. 


m 2 is still given by 


x+c 


2c 


_b 2 x -b 2 c 
S is still given by ~j 


ary 


a 2 y 


and tan 02 - 


b^_ 

yc 


However, slope of Rj = nij is not defined for x c, 
but 

x*-c x'-c 


x -c x -u __ j_ „ 

ton ^“ ■rT" ~ y ; y * ° 


where x*is the value of x where the extension of V 
intersects the x-axis. 

Further 

y'-V -b 2 c / 

= 2 — » y/° • 

x -c a 2 y 

Now solving of x* (y '= 0) we get 

-b 2 c , , 

-y = — 2 — ( x -c) 

J a z y 

or 

a 2 y 2 

Ti + c = x . 
b 2 c 

Then tan0j = ^ x = c, y ^ 0 . 

Since 0*, 0 2 ^ ± 90 > * 

if 0! = 02 » t 11611 


Now, let us consider the cases where x = ± c. 
(This will occur at four points on the ellipse. ) 


tan 0 A = tan 0 2 
i.e. , 



i 
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t^x - b 2 c = 0. 


Since x = c 
b 2 c 2 = b 2 c 2 . 


Since the steps are reversible, ^ 

Solution may be made similarly for x = -c. 
Where y = 0, the proof is trivial. 
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V. SELENOGRAPHY 



LUNAR MARIA TERRAIN MODEL FOR USE IN LUNAR SURFACE VEHICLE MOBILITY STUDIES 


By 


OthaH. Vaughan, Jr. and Donald Rose 


Symbol 

c 

0 


K 


0 


K 

c 


n 


Ki 

k 2 


H 

max 

W 

W t 

P 

Z 

y 


LIST OF SYMBOLS 


Symbol 


Definition 


Definition 


R 


Motion resistance, lbs 


Coefficient of soil cohesion, lbs/in 2 


R. 

1 


Internal resistance, lbs 


Angle of internal friction (between R^ 

grains) , degrees 

Modulus of soil deformation due to 
frictional ingredients of the soil, 

lbs/in n+ 2 b 

Modulus of soil deformation due to 
cohesive ingredients of the soil, HP 

lbs/in n+ 1 A 

Exponent related to the shape of the 
pressure-sinkage curve, dimension- ^ 

less 

Slip parameter reflecting the degree 
of soil compactness, i/in 


Total motion resistance, lbs 

Drive motor to track efficiency, 
per cent 

Width of bearing surface, per 
contact area, in 

Power, hp 

Total area in contact with the 
ground, in 2 

Length of bearing surface per 
contact area, in 

I. INTRODUCTION 


Slip parameter reflecting the fun- 
damental shape of the soil shear 
curve, dimensionless 

Maximum traction available from 
soil characteristics, lbs 

Weight of vehicle per contact area, 
lbs 

Total weight of the vehicle, lbs 


Ground pressure per track, lbs/in 2 
Sinkage in the soil, in 
Specific soil density, lbs/ in 3 


Lunar surface mobility capabilities of a roving 
vehicle are determined by the interactions of the 
vehicle with the lunar surface and the characteristics 
of the lunar soil (bearing strength, slope distri- 
bution, obstacle distribution, and many other para- 
meters). The proposed model is an attempt to pro- 
vide a more realistic maria model than the maria 
model as proposed by the study, ” Engineering Lunar 
Surface Model” [ 1] , and is based on small scale re- 
solution data as determined from Ranger VII photo- 
graphic and photometrically generated topographic 
information. The slope data were developed from 
the ”Topographic Map of a Small Area of the Moon’s 
Surface in Mare Cognitium” (Fig. 1) prepared by the 
Branch of Astrogeology of the United States Geolog- 
ical Survey (USGS) under NASA contract. 
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FIGURE i. TOPOGRAPHIC MAP OF A SMALL AREA OF THE MOON’S SURFACE IN MARE COGNITUM 


H. DISCUSSION constructed, aU of the positive and negative slopes 

were measured and tabulated. The distribution for 
A. DISTRIBUTION OF SLOPE FOR THE the positive and negative slopes for this particular 

MARIA MODEL traverse was 71. 4 per cent and 28. 6 per cent, re- 

spectively. Only the positive slope data are pre- 

Since the previous Engineering Lunar Model Sur- sented since these are the slopes that require energy 

face (ELMS) [ 1] was developed using small scale to be expended if the vehicle is not considered to use 

maps (1:1, 000,000) and supporting photographs to the power system for braking, 

produce a lunar surface model, it now is possible to 

use the same techniques to develop a large (1:80) B. SOIL PARAMETERS FOR MOBILITY 

scale engineering lunar model surface based on 

Ranger VII data. To construct this new surface 1. Although the slope frequency distribution 

model, a typical traverse for a lunar scientific diagram gives an indication of the terrain profile, it 

mission was constructed in the map area. This is still necessary to determine what soil parameters 

traverse was 100 meters in length and consisted of should be used with the model profile. Bekker [ 2] 

of seven different legs as shown in Figure 2. The in his work with locomotion and soil studies evaluated 

typical profile for each leg of the total traverse is soils in terms of a series of constants. These con- 

shown in Figures 3 and 4. After these profiles were stants are divided Into two categories: 


67 



a. Soil Strength 
c (cohesion) 

0 (friction) 

b. Soil Deformation 

K 0 (friction modulus of deformation) 

K (cohesive modulus of deformation) 
c 

n (exponent of sinkage) 

K (slip coefficient) 



FIGURE 2. LUNAR MARIA SURFACE TRAVERSE 


Elevation (m) 


3 — Traverse I * 

, End of Traverse 



Length of Troverse (m) 


FIGURE 3. LUNAR MARIA TRAVERSE PROFILE 


Elevation (m) 





2. Values for these constants have been de- 
termined for various types of terrestrial soils, and 
values have been postulated for lunar soils based on 
the terrestrial laboratory data. Bekker [ 2] indicates 
that the lunar soil, in the absence of an atmosphere 
and the absence of water, will probably not consist 
of the saturated plastic clay types, but will probably 
be similar to sand, gravel, or dry mineral powder. 
Many other investigators of the lunar surface and its 
characteristics [ 2 - 12 ] also postulate (based on 
thermal, photometric, and radar data) that the soils 
will probably consist of gravel, sand, lava, pumice, 
or even fine powder. If the lunar soils are of the 
granular type, they will probably not exhibit the 
cohesive forces which result from moisture present 
between the grains. Laboratory data [ 13] have 
indicated that the modulus of cohesion ( K^) is equal 


FIGURE 4. LUNAR MARIA TRAVERSE PROFILE 

to zero for both pumice and quartz sand in atmos- 
phere and in vacuum. The soil properties selected 

for the ELMS model were c = 0, K = 0, and 0 = 32* 

c i 

typical values for a cohesionless dry sand type soil. 
Laboratory data [ 14] indicate that the value of n 
(exponent of sinkage) for terrestrial soils ranges 
from a minimum of 0 . 5 for loose, minimum soils to 
a maximum of 1 . 25 for hard-packed high strength 
soils. These experimental data also imply that the 

value of K (friction modulus of deformation) ranges 

0 

from a minimum of 0. 1 for a fine pumice or a fluffy, 
dust type soil to a maximum of six for a hard, dry- 
packed sand. In selecting these constants for the 
ELMS, the values for n and K for each slope were 
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selected such that the small slope material would 
have values of n = 0. 5, and K^O. 5, respectively, 
increasing to values of n = 1. 25 and K^-6 for slopes 
of 25* . Slopes larger than 25* were considered hard 
surfaces. Since the soil constants for the lunar sur- 
face are still largely unknown and since many thermal, 
radar, and photometric studies have indicated that the 
lunar surface may be covered with a fine dust or ex- 
hibit the soil characteristics of dry sand, then the 
soil constants as used in the ELMS model [ 1] are 
considered satisfactory at the present time and will 
be used with this model slope distribution. The re- 
vised lunar maria terrain model is presented in 
Table I. 

TABLE L LUNAR MARIA TERRAIN MODEL 


Slope Angle* 
(Degrees) 

Percent of 
Total Distance 

Soil Constants* * 

0-2 

44.7% 

K<p = 0. 5 

2-4 

4.1% 

n = 0. 5 

4-6 

5. 6% 

(Loose fluffy dust) 

6-10 

5.4% 

= 1. 0 
n = 0. 75 

10 - 20 

7.9 % 

= 3. 0 
n = 1. 0 
( Loose dry sand) 

20 - 30 

3.0% 

K 0 - 6. 0 
n =1.25 

(Hard compacted sand) 

30 - 35 

0. 4% 

Hard surface 


* The distribution of slopes for this model are the distribution for the positive 
slopes. Negative slopes are not considered since they do not increase the power 
requirements unless the power system is used for braking. 


** c = 0 

K = 0 

c 

0 = 32 


C. SLIPPAGE FACTORS FOR MOBILITY 

In addition to the values of c, K , K , n, and 

c <\> 

(f>, in terms of mobility, one must consider some 
additional parameters, the slippage factors K t and 
K 2 . These slippage factors, as proposed by Bekker, 
and other expressions as defined by Janosi and Ha- 
namoto [ 15] are the result of the interactions of the 
particular soil and wheel combinations. Terrestrial 
data have indicated that K A values range from 0. 3 


for average soils to approximately 1 . 0 for the 
brittle or coherent soils, while K 2 values range from 
1. 0 for the brittle or coherent soils to approximately 
3. 0 for loose dry sands. The values for K t and K 2 
for relatively smooth wheels, as selected for the 
ELMS [ 1] were K t = 0. 2 and K 2 = 1. 25. Computer 
studies [ 16] have shown that as the value of Ki 
decreases from the larger values to smaller values, 
the effect is to increase the amount of slip for a 
given soil thrust; however, changes in the value of 
K 2 does not afreet the slippage, but tends to increase 
or decrease the soil thrust for a given slip. Thus, 
by selecting a low value for K 1 and K 2f we will have 
pessimistic values of slip for our lunar soil model; 
and by a careful selection of a certain set of c , , 

K , n, and d > , it is possible to establish a model 

<t> 

which represents a "worst case" moon surface from 
the standpoint of determining mobility per- 
formance estimates. Since no soil characteristics 
data have been determined from the Ranger photo- 
graphic data, we will use the same soil slippage 
factors as recommended by ELMS [ 1] for our model 
slippage factors. 

D. ENERGY REQUIRED FOR MOBILITY 

Through the use of computer programs, in- 
corporating the equations as developed by Bekker [ 2] 
to determine mobility energy requirements and using 
this maria model, the energy requirements for a 
vehicle traversing a smooth surface can be calculated. 
The mobility energy, as determined by this maria 
model, does not include that energy which is re- 
quired to negotiate obstacles in the maria. An 
illustration on the use of the model constants and 
equations used in a mobility analysis can be found 
in Appendix A. 

E. OBSTACLE DISTRIBUTION 

The obstacle distribution for this particular 
maria has not yet been developed. Although the 
Ranger data did not show the surface littered with 
blocky debris, one must remember that lunar maria 
in some areas may still have some obstacles which 
are not recognizable due to the resolution of the 
Ranger VII data. 


m. CONCLUSIONS 


Although this model is developed from a small 
area within Mare Cognitium, it seems to represent 
a typical maria model for use in future lunar sci- 
entific mission planning studies. When better data 
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on the characteristics of the lunar surface are ob- 
tained, then many different traverses will be con- 
structed and a more realistic lunar maria surface 
model will be developed. The major difference 
between this model and the ELMS [ 1] is the slope 


distributions. This proposed maria surface model 
produces a more rugged terrain for use in scientific 
mission planning studies and appears to be more 
realistic since it is based on a topographic map pre- 
pared from one meter resolution data. 


APPENDIX A 


The soft soil mobility of an off-the-road land 
roving vehicle is determined by marry vehicle and soil 
parameters such as number of wheels or tracks, 
shape and size of ground contact area, center of grav- 
ity location, soil characteristics, vehicle speed, 
etc. The following simplified analysis will be used 
to illustrate how some of the soil constants as post- 
ulated can be used in a soft soil mobility analysis. 

I. Statement of the Problem 

A lunar roving vehicle weighing 2, 328 pounds on 
the lunar surface is equipped with tracks for mobility. 
The vehicle tracks in contact with the surface each 
have a width of 12 inches and a length of 97 inches. 
Each track supports a load of 1, 164 pounds. The 
centerline of the tracks is 12 feet apart, and there 
is no elastic deformation of the track. The soil model 
will be the loose fluffy dust soil (K^-0. 5, n = 0. 5, 

K c = 0, c = 0, and <p = 32*) . The vehicle will operate 

on level terrain, and the slippage effects between 
soil and tracks will be neglected. Since the soil 
constants are normally stated in English units, the 
problem will be solved using English units; however, 
the results will be given in English and metric units. 

IL To Be Determined 

A. The maximum tractive effort available from 
the soil. 

B. The sinkage of the vehicle in the soil. 

C. The total resistance (R ) created by the soil 

l 

and the vehicle during motion. 

D. The horsepower (HP) required to drive the 
vehicle at 5 miles/hour over a loose, fluffy dust, 
level terrain assuming a power transfer efficiency 
of 80 per cent (the drive motor to the track). 


HI. Solution 

A. To determine the maximum tractive effort 

(H ) produced by the soil, we use the Bekker 
' max' 

equation [ 14] for a tracked vehicle while recognizing 
that it neglects the effect of slippage between the 
tracks and the soil: 

H = Ac + tan <£ 

max t 

since c = 0 from the model, the equation becomes 

H = W tan 4> 

max t 

where </> is the angle of internal friction and is the 
total weight of the vehicle, 

using W t = 2, 328 lbs. and 0 = 32 # from the model, 

then H = 1,455 lbs. (659. 98 kg). 

max 

Thus in frictional soils (c = 0) , the tractive 
effort depends only on the vehicle weight and the 
angle of internal friction and not the contact area. 

Since the tracked vehicle as stated in the problem 
has a large track length with respect to the width, the 
effect of the slippage constants K t and K 2 can be ne- 
glected in this simplified analysis. However, if a 
comparison must be made of the track versus the 
wheel from an overall mobility standpoint (various 
soils, obstacles, etc. ) , then these constants must be 
considered in the mobility analysis. 

B. The depth of sinkage (Z) can be determined 
using the Bekker equation [ 2] for the sinkage of a 
track or a high deflection elastic wheel in soils: 
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using b = 12 inches, 
P = 1. 0 lbs/in 2 


where 


K= 0. 5, and 


W 

t = 2,328 

2bl (2) (12) (97) 


1. 0 lb/in 2 (70. 31 g/cm 2 ). 


Since K = 0 from the model, the equation reduces to 
c 


n= 0. 5 
then 

R^ = 32 lbs/ track (14. 52 kg/ track). 



using P = 1. 0 lbs/in 2 , n = 0. 5, and K » 0. 5 then the 

<P 

sinkage/track 

Z = 4 inches sinkage (10. 16 cm). 

The vehicle track will sink into the soil 4 inches 
(10. 16 cm) , before it will exert a bearing pressure 

of 1. 0 lb/in 2 (70. 31 g/cm 2 ). 

C. The total resistance (R^) to motion created 

by the soil and the vehicle during its motion is due to 
the soil compaction, resistance due to soil bulldozing, 
and the internal resistance due to friction in the 
bearings, links, etc. of the track. 

1. The rolling resistance due to compaction of 
the soil can be determined by using the Bekker 
equation [ 18] . 


n + 1 



The motion resistance due to both tracks is equal to 
64 pounds, (29.03 kg). 

2. The bulldozing resistance of a track or wheel 
can be determined using the equation as recommended 
by Hegedus [ 17] . 

R^ = 1/2 y bz 2 tan 2 (45° + <p/ 2 + 2cbz tan (45° + <p/ 2). 

Since c = 0, the equation reduces to 
Rk = 1/2 y bz 2 tan 2 (45° + <p/2). 

Using Z = 4 inches, b = 12 inches, <p = 32°, and with 
a loose fluffy dust we assume y - . 01 lbs/in 3 . 

Then 

= 3. 12 lbs/track (1. 42 kg/track) . 

Thus, the total bulldozing resistance due to both 
tracks is equal to 6. 24 lbs. , (2. 83 kg) . 

3. The internal resistance (R^) due to the 

friction in bearings, links, etc. , of the track can be 
determined as follows: 

R i = fw t 


where 

b = width of the bearing surface per contact area. 
Since K c = 0 from the model, the equation becomes 


n+ 1 



where 

f = coefficient of internal friction, approximately 0. 04 
for a small tracked vehicle, 

using 

W= 2,328 lbs and f = 0.04 
then 

R. = 93. 12 lbs. (42. 24 kg) for both tracks. 
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4. The total resistance (R^) is the sum of all 
these resistances: 

R = R + R + R. 
t c b 1 

R t = 163.36 lbs. (74.10 kg). 

D. The horsepower (HP) required to drive the 
vehicle at 5 mph over a loose, fluffy dust, level ter- 
rain can be determined as follows: 

(RJ (V) 

HP = 

(550) (E) 

(Total Resistance) (Velocity) 

(Drive motor to trackefficiency) (550ft-lbs/ sec-hp) 

using 

R^ = 163. 36 lbs. , V = 7. 30 ft/sec (5 mph) , and a 
drive motor to track efficiency of 80 per cent. 

Then, the required horsepower is 2. 71 hp, (2. 02 kw). 

This analysis has demonstrated how the model 
soil constants can be used to determine some of the 
capabilities of a vehicle. In a complete mobility a- 
nalysis, many other parameters must be determined 
before the vehicle mobility capabilities can be estab- 
lished. 
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VI. SPACE ENVIRONMENT 



REVISED METEOROID FLUX AND PUNCTURE MODELS 


By 

Charles C. Dalton 


SUMMARY Symbol Definition 

Interrelated cubic polynomial models are con- p 

structed for log puncture flux versus log effective 
target thickness and for log mass -cumulative inci- 
dent flux versus particle mass, which give the best t 

agreement between meteoroid measurement satellite 
data and photographic meteor data without relaxing C 

the assumption that p/d is independent of d for parti- 
cles of different velocity, where p and d are punctur- 
able thickness and particle diameter, respectively. E 

An accurate account of the expected mean hazard 
for meteoroid puncture in earth orbit is needed for c 

design and operations planning purposes for future 
programs. The model developed below is thought to 
give mean puncture flux accurate to within a factor 2 v 
probable error with target equivalent thickness be- 
tween 15 and 10 -4 centimeters of aluminum. The in 

corresponding model for mass -cumulative incident 
flux is also developed. 

0 

LIST OF SYMBOLS 


Symbol 

Definition 

P 2 and 0 3 

Respectively, the slope and intercept 
of the tangent to the curve for log 

log 

Logarithm with base ten. 


F^ versus log m. 

P 

Thickness of a just -punctur able sheet 
in centimeters. 


Logarithm of the puncture flux en- 
hancement factor due to the statisti- 
cal variation of the particle and im- 

d 

Projectile diameter in centimeters. 


pact parameters. 

V 

Meteoroid air -entry velocity in 
kilometers per second. 

x 2> X 

Angle of impact of incident and of 
puncturing meteoroids, respectively, 
with respect to the surface normal 

V 

m 

f(v) 

Meteoroid lunar impact velocity in 
kilometers per second. 

Probability density function for v. 

H t 

in radians. 

Material hardness in Brinell units. 
Material parameter , (2/3) log 

m 

Mass of a meteoroid in grams. 


(p t C t } " 1 * 162 ‘ 

F > 

Flux of meteoroids of mass equal to 
or greater than m; mean number of 

y 

Target parameter, Iq + log p. 


hits per square meteor per second 
per 211 steradians without earth 

a 0 , . . . a 3 

Constants . 


shielding. 

z 

Log 0, a 0 + a t y + a 2 y 2 + a 3 y 3 . 


Crater depth in centimeters for a 
semi -infinite target. 

Target subscript. 

Velocity of sound in kilometers per 
second in an unbounded solid. 

10 -8 x Young ! s modulus in kilograms 
per square centimeter. 

Ductility, relative elongation in 2- 
inch gauge length at fracture*. 

Poisson ! s ratio. 

Meteoroid mass in grams nominally 
sufficient to puncture. 

Puncture flux; mean number per 
second per square meter of surface 
with 2n steradian exposure. 
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LIST OF SYMBOLS (Concluded) 


f («) and III Gravity - Free Component 


Symbol 

R 

H 

t 


At 

A 

P 


Definition 

Radius of effectively shielding atmos- 
phere, 6487 km. 

( 2 ) (Perigee + Apogee) -109 km, 
effective height above atmosphere. 

(i) [1 + cos sin 1 [R/ (R+H) ] , rela- 
tive exposability of randomly oriented 
surface [ 18] . 

Second -square -meters (partially 
shielded) . 



Centimeters thickness allowance for 

non -metallic mass in capacitor FIGURE 1. PROBABILITY DENSITY FUNCTION 

structure. FOR METEOROID AIR -ENTRY VELOCITY 


M. L. Maximum likelihood value. 

G 10“® x torsion modulus in kilograms 

per square centimeter. 


F g and 0 Mean incident and puncture flux, 

respectively, per second per square 
meter of a randomly oriented sur- 
face near the earth ( £ = £ ) . 


I. DISCUSSION AND RESULTS 


It had been expected previously [1] that meteo- 
roids with air -entry velocity only slightly in excess 
of escape velocity could be relatively so much more 
abundant, especially for those of small mass, that 
their flux in cislunar space might be strongly depen- 
dent on gravitational potential. But a weighted anal- 
ysis of photographic meteor data [2] shows no signi- 
ficant statistical relation between mass and velocity, 
and gives the velocity distribution shown in Figure 1, 
The weighted mean log velocity is log 27 kilometers 
per second; and only 25 percent of the total flux is 
due to the gravity-effect factor. The remaining 
component, shown also in Figure i, can be trans- 
formed into the distribution of lunar impact velocity 
shown in Figure 2, with median coincidentally 27 
kilometers per second. In the same analysis the 
mean air -entry meteoroid flux was found to be 


F 


-14.49 -1.34 

10 m 


( 1 ) 


Ht«) 



FIGURE 2. PROBABILITY DENSITY FUNCTION OF 
LUNAR -SURFACE IMPACT VELOCITY FOR 
PRIMARY METEOROIDS 

when the particle material density is f 3] 

p = 0.44 (2) 

and the least -countable mass m is in the interval 
-113 -0 23 

10 ^ m ^ 10 ’ .In establishing Equation ( 1) 

the mass of the zero -visual -magnitude -30 -kilometer - 
per -second meteor has been increased 0. 2 order of 
magnitude over Whipple’s [3] one-gram estimate in 
consideration of the cube of the relative change in the 
85-90 kilometer pressure for Reference 4 ^versus 
Reference 5 which Whipple used (basis: Opik [6]). 
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Correspondingly, the mean lunar impact flux would 

be 10 m * ,34 # When allowance is made for 

the mean solid angle of exposure for a surface in 
orbit , the mean flux will not differ appreciably from 
the fixed -level -surface values given above. 

In a survey and multivariable analysis of available 
laboratory hypervelocity impact data, Reismann, 
Donahue and Burkitt [7] found that in the highest of 
three velocity regimes the multiple correlation co- 
efficient is nearly maximized (to 0.989) by repre- 
senting the relation between projectile and target 
parameters for impact at normal incidence as fol- 
lows: 


p o /d = 


2.5{v/C t ) 0 - 50 (E/E t ) 1 - 31 (v/^) 8 - 0 (e/ et ) 0 - 43 (3) 


The relation between puncturable target thick- 
ness p and crater depth p^ has been treated in 

considerable analytical depth by Andriankin and 
Stepanov [ 8] . Their results indicate that the ratio 
p/p Q is a function of the projectile and target ma- 
terial parameters and of the impact velocity, varying 
between somewhat more than i to somewhat less than 
2 for impact velocity above about 6 kilometers per 
second. Herrmann and Jones [9] illustrate experi- 
mental results (which they attribute to Kinard, et 
al. at NASA, Langley) and an empirical formula in- 
dicating that, as target thickness is decreased toward 
the value p, the value of p/p q approaches 

p/p Q — (i/1.3) 2 + i = 1.59. (4) 


If Equations (3) and (4) are valid for meteoroid 
impact, then p/d is independent of d, and the slope 
of log <p versus log p should be three times the slope 
of log versus log m, where 0 is the puncture flux 

corresponding to incident flux F . Therefore, if the 

slope of log versus log m for photographic meteors 

in the vicinity of 10 1 gram mass, expressed by the 
constant -1. 34 in Equation ( 1) , were applicable also 
for the particles puncturing the meteoroid measure- 
ment satellites , then the slope of log 0 versus log p 
should be -4. 02 for the satellite puncture data. But 
in a previous working paper [10] the value -1.966 
for the slope of log <p versus log p was shown to fit 
the satellite data very well. At first the author [ 10] 
interpreted this discrepancy between the -4. 02 and 
-1. 966 results as indicating that p/d is not essen- 
tially independent of d for meteoroid impact circum- 
stances; but, after extensive consultation and further 


analysis, the present interpretation is that the slope 
of log F> versus log m is not essentially independent 
of m. 

Because puncture flux 0 versus effective target 
thickness p is of more direct technological interest 
than incident flux F versus particle mass m, even 
though the opposite might be true for scientific in- 
terests, a model should be developed first for 0 ver- 
sus p by using the satellite puncture data and com- 
puted values of 0 and p corresponding to the photo- 
graphic meteor data, Equation (1). The relation 
between incident flux F and puncture flux 0 must be 
established through 

F > = iO^ 3 n/ 2 (5) 

and 

<P = l(f 3 + fti m^ 2 , (6) 

where 0 2 and 0 3 are, respectively^ the slope and 
intercept of the tangent to the curve for log F ver- 
sus log m, m is the mass of the meteoroid winch, 
with mean values of density p, velocity v, and angle 
of incidence x 2 with respect to the surface normal, 
is just sufficient to puncture the target of interest, 

and 10^ 5 is the puncture flux enhancement factor [il] 
due to the statistical distribution of the parameters 
p, v, and x 2 . 

Recently, Equation (3) had been preferred [ 1 , 

10, 11, 12] as a basis for the relative puncturability 
of targets of different materials. That preference 
was based both on the detailed involvement of vari- 
ous material parameters and on the fact that the ve- 
locity exponent had the same value as had been pre- 
ferred previously by the author [13, 14] because it 
was intermediate between the values supported by 
different authors [9, 15] at that time. The author 
[13, 14] had preferred 

P o = 10°- 434 (pm/p t H t ) 1/3 (v cosx/ (7) 

where is the material hardness in Brinell units. 
Then, by Equations (1), (2), (4), and (7), a mete - 

roid with nominal mass in = 10 gram, nominal 
density p = 0.44, nominal velocity v = 27 kilo- 
meters per second, and nominal impact angle 
x 2 = n/4 radians impacting a sheet of 2024 T-3 
aluminum with density p^ = 2.77 and Brinell hard- 
ness H =120 would just puncture the sheet of thick- 
ness p — 10 or o. 96 centimeter. For the 

same conditions except with the intermediate value 
/3 2 = -1.23, and Equation (7), the author [11] found 
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* the puncture flux enhancement factor to be l(/* s = 

10 0,70 . Then, with Equations (1) , (5) , and (6) , 

the mean puncture flux for m = 10 1 would be 
-14 45 

0 - 10 . A line between this point (log $, 

log p) = (-14.45, -0.017) and the corresponding 
point for the thickest Pegasus target would have a 
slope of -4.02 = 3 (-1.34), which would support the 
extrapolation of Equation (1) if it were not for the 
fact that the puncture data for thinner targets (both 
on Pegasus and on Explorers XVI and XXIII) defin- 
itively establish a lower slope. It does not seem 
reasonable to accept an abrupt change in slope from 
about -4 to about -2 in the immediate vicinity of the 
Pegasus data point. Then it must be assumed that 
puncturable thickness p varies more steeply with 
velocity v than would be indicated by Equation (7) . 
Therefore^the author prefers to abandon both of the 
formulations in Equations (3) and (7) for meteoroid 
puncture. 

The NASA -Ames penetration criterion (by 
Summers [16] in 1959) has been recommended re- 
cently for meteoroid puncture analysis by Dohnanyi 
[17] ; i.e. , 

1/3 

P q - 3.56{(mp/2p t 2 ) [v(cosx 2 )/C t ] 2 } . (8) 

Together, Equations (4) and (8) indicate that p/d is 
again independent of d. With C = 6.25 from [7] 
for 2024T-3 aluminum, the nominal meteoroid, with 

m = 10 and with the values for the other param- 
eters mentioned above, the puncturable thickness is 
0 229 

P = 10 or 1.69 centimeters, 76 percent thicker 
than that which was inferred by Equation (7) . 

Equations (4) and (8) lead to some revision in 

the puncture flux enhancement factor 10^ 5 and in the 
probability density function of the angle of impact 
x of a puncturing meteoroid, as derived in [ 11] . 

It seems appropriate to retain the same value 0. 52 
for the standard deviation of log puncturable mass 
instead of increasing it to 0. 57 to account for the 
stronger dependence on velocity in Equation (8). 

This accommodation is supported by a reduction in 
the standard deviation of leg material density from 
0. 44 to 0. 38. Then, by the derivation [11] , differ- 
ential puncture flux d <f> must be multiplied by 

e 4 (cos x) 2 . The results of this alteration are 
that 

Pi = 0.63 (9) 


Therefore, by Equations (5), (6), (9), and (10), 
the puncture flux enhancement factor is 4, 26 instead 
of 5. 05 and the median impact angle x for the punc- 
turing meteoroids is 34. 7 instead of 33.5 degrees 
from the normal. 


By Equations (1), (5), (6), and (9), the mean 
puncture flux which is inferred by Equations (4) and 


-1 _\a 50 

(8) for m = 10 gram is <f> = 10 * , with p = 

0 229 

10 * for 2024T-3 aluminum. However, because 
Equation (1) represents the photographic meteor data 


—1 13 —0 23 

over the interval 10 * m ^ 10 * , the nomi- 

nally puncturing mass which is represented best is 


Hr. to" 0 - 68 w, th f -..T 13 - 579 


This corresponds to the point (log <£>, log p) = 
(-12.949, 0.336) for 2024T-3 aluminum. 


By introducing a material constant k 4 and a tar- 
get parameter k 4 + log p and the values of the param- 
eters already used with Equations (4) and (8) , one 
can write 


- _ ..3ki - 0.678 3 
m = 10 p 

(11) 

= (2/3) log (p t C ) - 1.162 

(12) 

y = k 4 + log p 

(13) 

= (1/3) log m + 0.226 

(14) 

Z = log (f> 

(15) 

= a 0 + a t y + a 2 y 2 + a 3 y 3 

(16) 

where a 0 a 3 are constants to be determined be- 

low. Also, 

02 = ( 1/3) (a t + 2a 2 y + 3a 3 y 2 ) 

(17) 

0 3 = z - 0.63 - /? 2 log m . 

(18) 


Therefore, with Equations (5) and (14) through (18), 
one can determine the mass -cumulative flux F as a 
function of mass m. > 

The constant term in Equation ( 12) for the ma- 
terial constant k^ has been chosen so that the target 
parameter y in Equations (13) and (14) vanishes 
for the photographic meteor value log m = -0. 68 and 
log <p = -12.949. Therefore^ Equation (16), 


f(x) 


1.77 (cos x) 1,54 


sin 2 x. 


( 10 ) 


a = -12.95 
o 


(19) 
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and, by Equations (i) and (li) , 

ai = - 4.02 . ( 20 ) 

The values of the other constants a 2 and a 3 in Equations 
(16) and (18) must be determined by the satellite 
puncture data. 

Puncture data for five meteoroid measurement 
satellites , compiled in Table 1 , adjusted to remove 
earth -shielding according to the mean solid angle of 
exposure, are plotted in Figure 3. The values of the 
physical constants for the materials are also given 
in Table 1. The values used for p for the Pegasus 



-4 -3 - 2-10 1 


FIGURE 3. METEOROID PUNCTURE FLUX VS 
EFFECTIVE TARGET PARAMETER 

-3 

sensors include a 10 - centimeter increment in con- 

sideration of the non -metallic components of the ca- 
pacitor puncture sensors [1] . 

The values of log (j > , as shown in Table 1 and 
Figure 3, have been corrected for geometric shield- 
ing [18] by dividing the mean puncture flux by £. 
Poisson's ratio for the Explorer XVI transducer ma- 
terial was calculated from Baumeister [ 19] : 

v = (E - 2G)/2G. (21) 

Values for velocity of sound in the transducer ma- 
terials were calculated with the extended -medium 
formula [7, 20] 

i 

C =9.90 { (E/p) [ (i-v)/(i+v) (l-2i>)]} 2 (22) 

instead of from the corresponding formula for thin 
plates [21] . 


The two conditions for the satellite puncture re- 
sults from Table 1 and Figure 3 which seem most 
significant for the determination of the constants a 2 
and a 3 in Equation (16) are 1) the values -1.717 
and -7.171 for the abscissa y and the mean ordinate 
z, respectively, for the thickest Pegasus target, and 
2) the values -2. 383 and -1. 58 for the abscissa y , 
and the Hastings [27] maximum likelihood slope 
dz/dy, respectively, for the Explorer XVI target of 
intermediate thickness* Therefore, by these two 
conditions and Equations (16), (19), and (20), 

a 2 = -0.424 < 23 ) 

a 3 = -0.025 . < 24 ) 

By Equations (13), (15), (16), (19), (20), ,(23), and 
(24) , the relation between puncture flux <p and target 
parameter (Iq + log p) is 

log <p = -12. 95 -4. 02 (k 4 + log p) -0. 424 (k 4 + log p) 2 
- 0.025 (k« + log p) 3 (25) 

where the material parameter k 4 is given by Equa- 
tion (12) . The corresponding numerical results are 
superimposed on Figure 3 for comparison with the 
satellite puncture data. It should be noted that, al- 
though the Pegasus puncture flux has been adjusted 
upward to correct for the approximately 30 percent 
earth shielding, it has not been corrected upward to 
correct for any counting loss, which might be about 
14 percent [ 29] . 

The values of the mass -cumulative flux F^ and 
mass m, corresponding to any point on the curve in 
Figure 3, can be computed most readily from 

log m = 3(k4+logp) -0.68 (25) 

log F > = log $ -0.63 . (26) 

For example, the thickest Pegasus target is being 
punctured by particles with about 1. 5 micrograms 
mass. 

The results, which may be accurate to within a 
factor 2 probable error, are illustrated further in 
Figure 4 for the mass -cumulative flux F (per ef- 
fectively exposed hemisphere). For better compari- 
son with other models. Figures 5 and 6 show the 

mass -cumulative flux F and puncture flux <t> for a 
• s s 

randomly oriented vehicle at low orbital altitude from 

the earth. 
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TABLE 1. DATA FOR TWELVE SATELLITE PUNCTURE TRANSDUCERS 


ENTITIES 

DATA AND SOURCES 

Satellite 

Explorer 

XVI 

Explorer 

XVI 

Explorer 

XVI 

Explorer 

xxm 

Explorer 

xxm 

Pegasus 

I 

Pegasus 

n 

Pegasus 

H 

Pegasus 

n 

Pegasus 

m 

Pegasus 

m 

Pegasus 

m 

Launch Date 

12-16-62 

12-16-62 

12-16-62 

11-6-64 

11-6-64 

2-16-65 

5-25-65 

5-25-65 

5-25-65 

7-30-65 

7-30-65 

7-30-65 


[22] 

[22] 

[22] 

[23] 

[23] 

[24] 

[25] 

[25] 

[25] 

[26] 

[26] 

[26] 

Perigee (km) 

750 

750 

750 

463 

463 

497 

506 

506 

506 

529 

529 

529 


[22] 

[22] 

[22] 

[23] | 

[23] 

[24] 

[25] 

[25] 

[25] 

[26] 

[26] 

[26] 

Apogee (km) 

1180 

1180 

1180 

980 

980 

745 

748 

748 

748 

532 

532 

532 

H 

856 

856 

856 

612 

612 

512 

518 

518 

518 

422 

422 

422 


[27] 

[27] 

[27] 

[23] 

[23] 

[28] 

[28] 

[28] 

128] 

[28] 

[28] 

[28] 

Punctures 

44 

11 

0 

24 

25 

125 

121 

18 

58 

71 

14 

41 


[27] Fin. 

[27] Fin. 

[27] Fin. 

[23] 

[23] 

[28] 

[28] 

[28] 

[28] 

[28] 

[28] 

[28] 

Tally Date 

7-22-63 

7-22-63 

7-22-63 

2-11-65 

2-11-65 

9-5-65 

9-2-65 

9-2-65 

9-2-65 

10-8-65 

10-8-65 

10-8-65 


[27] 

[27] 

[27] 

[23] 

[23] 

(28] 

[28] 

[28] 

[28] 

[28] 

[28] 

[28] 

Log At 

7.060 

6.738 

6.512 

6.731 

7.065 

7.953 

7.737 

8.014 

9.101 

7.587 

7.938 

8.951 

t 

0.734 

0.734 

0.734 

0.703 

0.703 

0.688 

0.689 

0. 689 

0.689 

0.672 

0. 672 

0.672 


[27] M.L. 

[27] M.L. 

[27] M.L. 










Log <t> 

-5.266 

-5.696 

S. 275 

-5.197 

-5.514 

-5.693 

-5. 492 

-6.597 

-7. 176 

-5. 563 

-6.619 

-7. 165 

10 4 (p^Ap) 

[27] 

29.2 

[27] 

54.6 

[27] 

127 

[23] 

25.4 

[23] 

50.8 

[29] 
38 ±3 

[29] 

38±3 

[29] 

203*25 

]29] 

406*50 

[29] 

38*3 

[29] 

203*25 

[29] 

406*50 







[1] 

[11 



[11 

[11 

[11 

10 4 (Ap) 

0 

0 

0 

0 

0 

10 

10 



10 

10 

10 

Log P 

-2.535 

-2.263 

-1.896 

-2. 595 

-2. 294 

-2.317 

-2.317 



-2. 317 

-1.671 

-1.381 


[22] 

[22] 

[22] 

[23] 

[23] 

[1] Al. 

[ 1] Al. 

[30] Al. 

[ 30] Al. 

[ 1] Al. 

[30] Al. 

[30] Al. 

Material 

Berylco 

25 

Berylco 

25 

Berylco 

25 

AISI 

302 

AISI 

302 

1100-H14 

1100-H14 

2024-T3 

2024-T3 

1100-H14 

2024 -T3 

2024 -T3 

Condition 

[22] 

Annealed 

[22] 

Annealed 

[22] 

Annealed 

[23] 
\ hard 

[23] 
^ hard 








Rockwell No. 

[22] 
B 60 

[22] 
B 60 

[22] 
B 60 










P t 

[31] 
8. 23 

[31] 

8.23 

— QO 

[30] 

7.83 

[30] 

7.83 

m 

2.72 

[7] 

2.72 

[7] 

2.77 

[7] 

2.77 

[7] 

2.72 

[7] 

2.77 

[7] 

2.77 


[31] 

[31] 

[31] 

[19] 

[19] 

[7] 

[7] 

[7] 

[7] 

[7] 

[7] 

[7] 

E t 

1. 195 

1. 195 

1.195 

1.941 

1.941 

0.703 

0.703 

0.745 

0.745 

0.703 

0.745 

0.745 


[31] 

[31] 

[31] 

[19] 

[19] 








°t 

0.457 

0.457 

0.457 

0.745 

0. 745 












[19] 

[19] 

m 

[7] 

[7] 

[7] 

[7] 

[7] 

[7] 

v t 

0. 308 

0.308 

0.308 

0.305 

0.305 

0. 33 

0. 33 

0.33 

0.33 

0.33 

0.33 

0. 33 

C t 

4.44 

4.44 

4.44 

5.73 

5.73 

6. 12 

f 7] 
6. 12 

[7] 

6.25 

[7] 

6.25 

[7] 

6.12 

[7] 

6.25 

[7] 

6.25 


[32] 

[32] 

[32] 

[30] 

[30] 

[7] 

[7] 

[7] 

[7] 

m 

[71 

[7] 

€ t 

0.46 

0.46 

0.46 

15 

0. 15 

0. 20 

0.20 

0. 18 

0.18 

0.20 

0. 18 

0.18 


-0. 120 

-0.120 

-0. 120 

-0.061 

-0. 061 

-0.348 

-0.348 

-0.336 

-0.336 

-0. 348 

-0.336 

-0.336 

k< + log p 

-2.655 

-2. 383 

-2. 016 

-2.656 

-2. 355 

-2. 665 

-2. 665 

-2.007 

-1.717 

-2. 665 

-2.007 

-1.717 

Sensor 

Structure 

Pressure 

Can 

Pressure 

Can 

Pressure 

Can 

Pressure 

Can 

Pressure 

Can 

Capac- 

itor 

Sand- 

wich 

Capac- 

itor 

Sand- 

wich 

Capac- 

itor 

Sand- 

wich 

Capac- 

itor 

Sand- 

wich 

Capac- 

itor 

Sand- 

wich 

Capac- 

itor 

Sand- 

wich 

Capac- 

itor 

Sand- 

wich 
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FIGURE 4. METEOROID MASS CUMULATIVE 
TERRESTRIAL MEAN INFLUX 


^ s , Punctures Per Square Meter Per Sec (with 5 = 1/2) 


-WHIPPLE -1963 (3) 


-MSFC-1964 (14) 

• MSFC-1965 (12) 


PROPOSED 
MODEL ; 


-BELLCOMM -1963 (35) 


OMSF -APOLLO -1965 1 MSFC-1966 

Upper 90% Confidence (34) ( 1 0) 


BELLCOMM - 1965 (17) \\ 


F s , Impacts Per Square Meter Per Sec (with £ = 1/2) 




10'* K) 10 ‘ 10 w 

m, Meteoroid Mass in Grams 


FIGURE 5. MEAN IMPACT FLUX OF 
METEOROIDS WITH MASS > M FOR A 
RANDOMLY ORIENTED VEHICLE IN A 
NEAR EARTH ORBIT 


10’ 4 10" 1 10 ° 
p, Centimeters Wall Thickness of Aluminum 2024T-3 


FIGURE 6. MEAN PUNCTURE FLUX FOR A 
RANDOMLY ORIENTED VEHICLE IN A NEAR 
EARTH ORBIT 


II. EQUIVALENT PUNCTURABILITY 


By Equation ( li) , if sheets of different materials 
must have the same meteoroid punctur ability , then 
their target parameters must have the same value; 

i.e. , (10^ 4 p) is invariant. For example, sheets of 
2024-T3 aluminum must be 1.65 times as thick as for 
annealed Berylco 25 beryllium copper and 1. 89 times 
as for \ hard AISI 302 stainless steel. But the nec- 
essary weight of aluminum is still less in either case 
and soft aluminum (1100-H14) is nearly as effective 
as 2024T-3. 











( 31 ) 


in. interplanetary space 


Most of the puncture hazard In cislunar space, 
and in interplanetary space near or between the 
orbits of Earth and Mercury, is due to particles of 
cometary rather than asteroidal origin. This follows 
from McCoy's [36] preference for Beard's [37] 
estimate that the flux of cometary particles should 
be proportional to D -8 / 2 , where D is the heliocentric 
distance in astronomical units, and from Parldnson's 
[38] model for the distribution of asteroidal parti- 
cles with respect to mass m and distance D. For 
the flux of vehicle-puncturing asteroidal particles , 
Parkinson [38] preferred the extrapolation of 
Hawkins' [39] model for large stones, giving 
10-15.93 m -1 particles per square meter per second 
for the terrestrial influx. By this formula the flux 
of photographic meteor particles larger than 10 -1 
gram is nearly two orders of magnitude below the 
corresponding value from Figure 4. This result 
agrees with the observation of Whipple and Hawkins 
[40] and of Jacchia and Whipple [41] about the rel- 
atively great predominance of cometary particles 
in samples of photographic meteors. Then 

F > ~ D- 3 / 2 . (27) 

The 27 kilometers per second median lunar im- 
pact velocity, for the distribution illustrated in 
Figure 2, is only slightly less than the heliocentric 
circular orbital velocity for the particular helio- 
centric distance. Anyway, it would seem appro- 
priate to assume that, at any other heliocentric 
distance D, nominal closing velocity is proportional 
toD- 1 / 2 ; i.e., 

v = 27 D -1 / 2 . (28) 

Then, by Equations (4) and (8) , the thickness p of a 
target sheet nominally puncturable by a meteoroid 
of particular mass is proportional to D"* 1 / 3 ; i.e. , 

p ~ D -1 / 3 . (29) 

Also, by Equations (4) and (8) , p 3 ~ mv 2 ; i.e. , 
together with the result in Equation (28) , it follows 
that the particle mass m nominally sufficient to 
puncture a particular target is 

m ~ D. (30) 

Therefore , puncture flux <f> is changed not only by the 
factor D 3 ' 2 due to the change in the mass-cumulative 
flux ( according to the model by McCoy [ 36] and 

Beard [37]) but also by the factor D 9(1 ° g F < )/d log m 


- D *. These two factors combine to give 

<p ~ D~ 5 / 2 . 

Figures 3 and 4, with respective ordinates 
log <p and log F >? presuppose that D = 1 astronomical 

unit. Then, by Equations (31) and (27) , it is nec- 
essary to subtract 2. 5 log D and 1. 5 log D from the 
ordinates of Figures 3 and 4, respectively, for 
applic ability to other heliocentric distances D. 

At heliocentric distances very much beyond the 
orbit of the Earth, the meteoroid hazard should be 
calculated separately for cometary and asteroidal 
populations of particles. The results for the 
cometary particles, adjusted with respect to D as 
already described in this section, are applicable in 
interplanetary space near or between the orbits of 
Mercury and Jupiter, except that for close orbits 
around Jupiter it is necessary to account for the 
significant effects of the gravitational field of the 
planet. 

Because of the relative sparsity of the aster- 
oidal meteoroids in the terrestrial influx, pertinent 
information about this population is less substantial. 
The nominal closing velocity will be assumed to be 
the same as is given as a junction of D in Equation 
(28). It seems appropriate for the vehicle- 
puncturing asteroidal particles to foUow Parkinson's 
[38] extrapolation of Hawkins' [39] model for the 
terrestrial cumulative influx of stones, 10” 15 * 93 m _i 
particles per square meter per second. Parkinson 
[38] preferred the following flux enhancement 
factors for these particles with respect to helio- 
centric distance: i. 0, 2 x 10 3 , 3 x 10 4 , and 20 for 
the orbits of Earth, Mars, asteroids (at 2. 5 AU 
point of max. ) , and Jupiter, respectively. The 
Jovian and terrestrial gravitationaUy enhanced flux 
factors, for these asteroidal particles, were given 
by Parkinson [38] as 170 and 2. 1, respectively. 

Opik [42] gives the average density of meteOritic 
stone as p = 3. 4. It is expected that a somewhat less 
dispersion of density for these particles may be just 
about equally compensated by a somewhat greater 
dispersion of closing velocity, so that the dispersion 
of the mass just sufficient to puncture a given target 
sheet may be essentially the same as for the come- 
tary particles. Then the value of the puncture flux 
enhancement factor given in Equation ( 9) should 
be valid for both populations. 

IV. CONCLUSIONS 

When a meteoroid impact target of interest is 
the puncturable equivalent of a sheet p centimeters 
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thick of a material with specific gravity and sonic 

velocity C^. kilometers per second, it has a material 

constant k 4 , a target parameter k 4 + log p, and is 
just puncturable under average impact circumstances 
by a meteoroid with mass in at heliocentric distance 
D astronimical units; i. e. , 

k 4 = (2/3) log (p t C ) - 1.162 (12) 

“ = i0 3 k 4 - 0.OT8 p 3 d. (11-30) 

In cislunar space, except for Earth or moon re- 
duction of the solid angle of exposure to a celestial 
hemisphere, the mean flux of primary meteoroids 
is nearly independent of distance beyond low orbital 
altitude and is reduced by only about 25 percent out 
to the boundary of the earth's gravitational domi- 
nance. 

For a surface with 211 ste radian exposure, the 
flux F> of particles with mass exceeding min 

Equation (11) is 0. 63 order of magnitude below the 
puncture flux 0, because of the statistical variation 
of the particle material density, velocity, and angle 
of impact; i. e. , 

F = 10 - °* 63 0. (26) 

> 

The mean puncture flux 0 per second per square 
meter per 2U steradian exposure is related to the 
target parameter k 4 + log p by 

log 0 = -12. 95 -4. 02 (k 4 + log p) -0. 424(k 4 + log p) 2 

-0. 025(k 4 + log p) 3 -2. 5 log D. (25-31) 

The result for Equations (25) and (26) are 
illustrated graphically in Figure 3 and 4 for a unit 
value of D. All the values for mass-cumulative 
flux F > and puncture flux 0 by the curves in these 

two figures are thought to be accurate to within a 
factor 2 probable error (normal distribution) , and 
constitutes the primary-meteoroid flux and puncture 
models proposed for use in design and operations 
planning in MSFC programs for Earth-orbit, 
cislunar and interplanetary space near or between 
the orbits of Earth and Mercury.* 

* Example. The flux values computed from 
Equations (25) and (26) are statistically nominal 
for the mean flux. To obtain the 99. 86 percent one- 
sided confidence value for the mean flux ( 3 standard 
deviations from the nominal value) , multiply by 
3(2/0. 6745) « 9. 


In interplanetary space much beyond the orbit 
of the Earth, it is appropriate to consider that the 
puncture hazard is due separately to the cometary 
and asteroidal populations of particles. The re- 
sults for the cometary particles, adjusted with re- 
spect to D as already described in this and the pre- 
ceding sections, are applicable in interplanetary 
space near or between the orbits of Mercury and 
Jupiter, except that for close orbits around Jupiter 
it is necessary to account for the significant effects 
of the gravitational field of the planet. The best 
available description of the asteroidal particle 
environment is given in the preceding section. 
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ON RADIATION PRESSURE AND ITS EFFECTS ON SATELLITE MOTION 


By 

W. H. Heybey 


SUMMARY 


The pressure of light falling on a plane surface 
is bound up with the time -averaged stress prevailing 
in the mov ing electromagnetic field, A pr^ssiirp ex- 
pression for elemental surfaces inclined to parallel 
light is derived with the stipulation that the photoelec- 
tric effect is absent. Other losses (through trans- 
mission and absorption) are briefly discussed. Fi- 
nite surfaces are supposed to be given in analytic 
terms; a method is described for finding the light 
force acting on them, and its point of attack. The 
motion of a satellite, otherwise undisturbed, under 
the impact of light is investigated in general. An ap- 
plication to an orbiting cylinder illustrates the meth- 
ods. 


i. experimental and theoretical basis 


for surfaces slanted to light direction. It will be 
carried out here in order to support the easy gen- 
eralization usually adopted, which is open to discus- 
sion if closely looked at. Reservations to be made 
will appear at several places, as will the need to 
average. 


One may start with the complete set of Maxwell’s 
equations ( including boundary conditions) written 
down for a nonconducting medium: 


i 9 5 

c =vx « 


i 

0 9 r=- v *£ 


V • D = 47 rp 

V • B = 0 . 


( 1 ) 


As a natural phenomenon the mechanical force 
accompanying the incidence of light is observed and 
measured with the use of a torsional balance ( See [ 1] , 
pp. 556-557) . It is found to be very small; however, 
over sufficiently long periods of time, it may produce 
unwanted disturbances in the motion of satellites sub- 
ject to rigid attitude requirements as, for example, 
those of an orbiting telescope whose line of vision 
should remain fixed within 0. i of an arc of a second 
or less in current estimates. 

Theoretical expressions for the radiation pres- 
sure are derived in many treatises on optics or theo- 
retical physics. In one way or another, the argument 
is mostly from the standpoint of the electromagnetic 
theory. One notable exception is that of Drude [2] , 
pages 488-491, where the proof is based on the gen- 
eral energy principle and thus avoids the conceptual 
difficulty inherent in field -body interaction. Since the 
demonstration has to rely on light energy converted 
into heat energy by an ideally black body , it cannot 
very well be adapted to reflecting surfaces and will 
not be used in what follows. 

Most of the textbooks readily accessible to the 
author restrict the case to light striking the surface 
at a right angle which is an exceptional incidence. 

None of them presents a direct and cogent derivation 


If the medium is isotropic and homogeneous, as we 
will and can assume here, the vectors D and B are 
related to the electric (E) and magnetic (H) field 
strengths by the equalities 

D = € E, B = p H (2) 

where the dielectricity coefficient and the permeabil- 
ity are constants. 

Sunlight falling on a metallic surface will liberate 
electrons from it, especially in empty space where 
the ultraviolet part of the radiation is strong. The 
surface charge (density p) thus created will interfere 
with the electromagnetic field, and some of the light 
energy will be lost in photoelectric work. For both 
reasons the surface density will be disregarded 
(p = 0) , so that V • E is considered as zero every- 
where. 


With these limitations equation (1) can be re- 
arranged to give 


f 


a 2 E 

at 5 


= AE 


(3) 


and thus assumes the form of the wave equation. In 
a vacuum (c = p = 1) the propagation velocity of 
the wave is given by the constant c which will be 
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If, in general, 


regarded as known. After a solution is established 
to equation (3) , the associated wave motion of the 
magnetic vector H follows from the second Max- 
wellian equation where B = p H . 

If a rectangular Cartesian system is used (unit 
vectors on the x,y,z-axes: i, j, k) , a particular so- 
lution of equation (3) can be set down as 

E = A E cos 2 7T ^ i = Eji. (4) 

It satisfies the boundary condition (V • E = 0, here 
3E t /3x = 0) and represents a plane wave with the 
z-axis as its normal. In an isotropic medium this 
normal is in the direction of the energy transmission; 
i.e. , it gives the ray direction. The wave length, 

X, and period, T, are linked to the ray propagation 
velocity by the relation 



One of the quantities X and T can be chosen arbi- 
trarily. 

The second Maxwellian equation can now be in- 
tegrated to give 

s ■ a e 7T cos 2 ’(t - k) 1 ■ Hli ' <5) 

The wave form and motion are that of the E-wave; 
but the amplitude 



is equal to the electric amplitude A only if € = p, 

Xli 

in particular, if the electromagnetic wave moves in 
a vacuum. 

With this solution the electric and magnetic field 
vectors are in the x- and y -directions , respectively , 
and the ray direction is normal to both. 

Plane wave motion can be said to prevail if a ve- 
hicle travels far enough from the sun to be subject to 
virtually parallel radiation. This assumption will be 
adopted as a further limitation. The force acting on 
a surface element, dS, can then be obtained with the 
use of the integrals (4) and (5) . It can be computed 
from the electromagnetic stress in the medium or, 
rather, in the ether underlying it in Faraday ! s and 
Maxwell’s view. 


E = iEj + jE 2 + kE 3 , 


the nine components of the electric stress tensor are 
known to be 


Y 

y 

z 


z 

Y 


z 

Z 


x 


X 

y 


= ~ <-Ei 2 + E 2 2 + E 3 2 ) 

O 7T 

= ~ (Ej 2 - E 2 2 + E 3 2 ) 

o 7T 

= (E t 2 + E 2 2 - E 3 2 ) 

= Z = - E 2 e 3 

y 4 7T 

= X = - Eg Ei 

Z 4 7T 



(6) 


In the case of the plane wave (4) , all the shear 
components are zero (E 2 = E 3 = 0) , while the nor- 
mal components add up to the stress vector 


E e ■ 17 • 

Expressions (6) also define the magnetic stress 

tensor, when the E. are replaced by the corresponding 
i 

quantities H.. The evaluation in terms of the wave (5) 
for which H 3 = Hi = 0 yields 

4 ■= 


Both the field stresses P E and P R are periodic 

functions of time. The radiation pressure to be com- 
puted from them becomes comparable to the constant 
found experimentally by taking the average over a 
period. Since 


T / 

f cos 2 7r( ; 

t z\ 

r x ) 

d * ■ i • 

0 ^ 

1 A / 



(E ‘\v^ A E 2 > (H 2 \v = it A E 2 • 

Thus, the two stress vectors combine to a common 
average 
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. e- <iy ♦ %> m 

av av 


This expression holds for any values of e and /z . 


Since the energy density at a given point in the 
(time -dependent) electromagnetic field of the waves 
(4) and (5) may be written as 

€Ej 2 + fjH 2 2 

8*7T * 

its average value is found as ~ — A 2 and gives the 

O 7T E 

magnitude of the average stress vector. For short, 
we may put 



and 


( 8 ) 


P= pk. 
- o 


(9) 


The quantity A 2 is proportional to the light in- 
E 

tensity. In polychromatic sunlight it depends on the 
wave length and should be regarded, in formula (8) , 
as an average taken over the spectrum. 

Expression (9) shows that the stress in the par- 
ticular electromagnetic field chosen is parallel to the 

light rays. The quantity p has the dimension of a 

o 

pressure ( — = ^° rCe ) , and will shortly be 
volume area 

taken as a mechanical pressure. Strictly speaking, 
it cannot be fully understood in what way a stress 
existing in the ether (i.e. , the void) could give rise 
to a ponder omotive force; observation only can sug- 
gest that it does so* . 


Consider a light beam of cross-sectional area 
dS Q (Fig. 1). At any point (A or A') in it, the oppo- 
site stresses P and -^P are active. Let the light 
at A be intercepted by a cross-sectional material 
surface. An open face of the stressed ether is then 
created, one of the stresses is freed, and the force 


M. Planck [3] after commenting on this point goes 
on to say that a similar difficulty exists in trying 
to visualize energy flow through a vacuum. 



ACTING ON PLANE ELEMENTAL SURFACES 


F = p dS k 

“OO 

is attacking the surface. Since it is normal to it, p Q 

can be interpreted as the pressure on the surface ex- 
erted by the force F. 


Regarding sunlight, the radiation pressure p^ 

at the outer fringes of the atmosphere has been de- 
termined as 


1. 374 x 10 g 
c 


dyne 
cm 2 * 


( 10 ) 


This figure depends on the distance of the Earth from 
the sun and is slightly different at perihel and aphel. 
The significant coefficient varies between 1. 351 and 
1. 397 (See [4] , p. 4 ) . For satellites moving in 
near -Earth orbits, one will often take the average 
value used in formula (10) . The constant c (« 3. 10 10 ) 
is the measure of the vacuum light speed in cm/sec. 

The radiation force on the whole Earth (assumed 
as a perfectly absorbing body) equals the weight of 
5.7 x 10 7 kg while the sun ! s attracting force corre- 
sponds to 3. 6 — 10 21 kg [5] . The insignificance of 
the radiation force is apparent. On the other hand, 
in the interior of fixed stars where the light intensity 
is immense, radiation may outweigh gravitation 
( Eddington) . 

The force F, as an outcome of the stressed state 
of the ether, is characteristic for it and cannot change 
either magnitude or direction if the beam is termi- 
nated by any arbitrary cut-off surface dS (at A 1 , 

Fig. 1) . * The force on dS may be given the form 

F = pdSk (11) 


If light behaved like a Newtonian fluid, the normal 
component only of F would act on dS. 
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where p is the pressure on dS. It follows that 
dS 

p =p o ¥ ' 

For any two plane surfaces lying in parallel projec- 
tion, the product of the area with the cosine of the 
angle the surface normal makes with the projecting 

direction is constant. Thus, dS = dS cos a and 

o 

p = p cos a . (12) 

o 

It appears that if the surface is illuminated under the 
angle of incidence, a 9 the pressure on it is the frac- 
tion p^ cos oi of the basic pressure p^. 

If perfect reflection is assumed ( ideally smooth 
surface, no transmission, no transformation of en- 
ergy ) , the outgoing light will have the same intensity 
as the incoming light (p^ remains unaltered). With 

plane waves the reflection angle is equal to the angle 

(r) 

of incidence, so that the force, _F , caused by re- 
flection will point in the direction marked k* in Fig- 
ure 1. The total force combining incidence and re- 
flection, 

F. . , = F (l) + F (r) = p cos a; dS (k + k') , 

- total “ ~ o 

is in direction of the vector 

t_ - k + k 1 . 

Since k and k* are of equal length, Jt is normal to the 
surface (pointing toward the interior) . The magni- 
tude | _t | is found from the trigonometric sine- 
theorem: 


— sin ol 

Thus, with t_ = |_t I n, 

F = 2p cos 2 0! dS n. (13) 

“t o — 

* 

The basic pressure p^ is multiplied by 2 cos 2 a. 


The factor cos 2 o' connected with oblique incidence 
is mentioned in [6] (p. 616) , where reflection is 
taken into account for normal incidence, too, so 
that p^ goes into 2p^. Cos 2 o: is given for incident 

light alone in [1] , p. 557. Ditchburn considers the 
normal component of the light force only. 


Physical reality seldom offers even quasi -ideal 
circumstances. On striking a surface, light will 
partly be transformed into different energy species , 
partly be reflected, specularly or diffusely. The full 
value of p^ will rarely prevail even at first incidence 

(as soon as there is measurable transparency). The 
dependency of transmission and absorption on wave 
length and angle of incidence is studied in [4] , where 
applications to a plane surface and to convex surfaces 
symmetric to incidence are also given (even to sur- 
faces containing zones of different material) . The 
purpose of the present paper does not permit limita- 
tion to preferred light directions; however, in view of 
the complexity of the corrections otherwise necessary , 
it seems advisable first to consider perfectly absorb- 
ing or perfectly reflecting surfaces only. Pertinent 
corrections to the pressure p^ can afterwards be in- 
corporated into the expressions; the integration over 
the surface then will have to be stepwise as a rule. 
When scattering was to be considered, a more in- 
volved mathematical approach would have to be em- 
ployed to account for the multitude of minute reflecting 
surface parts whose distribution, if taken at random, 
would make the process largely independent of inci- 
dence (aside from shadow profile changes). 

H. THE FORCES EXERTED BY PARALLEL 
LIGHT ON A FINITE SURFACE 

The light should radiate through an isotropic 
homogeneous medium (not necessarily through vac- 
uum; expression (9) has been obtained with arbitrary 
values of e and /i) . Photoelectric effects must be 
disregarded. The surface is allowed to contain per- 
fectly absorbing and perfectly reflecting portions 
which may be shaped and distributed at will. 

For an easy description of optically different sur- 
face zones, it is advisable to replace the ordinarily 
Cartesian surface equation by a parametric repre- 
sentation of the surface points : 

x = x(<r, r ) 

y = y(<r. t) < 14 ) 

z = z(<r f t) . 

Each of the variables cr and r w ill move within a 
certain "natural" interval for the system (14) to de- 
fine either the whole surface or any of its parts in the 
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Cartesian system. In the (a, r) -plane the surface 
will usually appear as a rectangle; its parts may be 
odd -shaped areas, in which case the limits of one 
natural interval will vary when the other variable 
moves in its own interval. 

The local interior surface normal, a basic req- 
uisite, can be introduced through the unit vector 

n = in A + jn 2 + kn 3 

in its direction. The components n. follow from the 
system (i4) through the relations 


With perfect reflection, equality (13) applies to 

give 

d 2 F^= 2p^(n- s) 2 N dadr n . (20) 

While the vector _s_ can be taken here as fixed rela- 
tive to the surface , n and N depend on a and r . 
(On occasion, N may be found constant ; _n is so for 
the plane only. ) 

With a fully illuminated surface, the integrations 
extend over the natural intervals of a and t . Sha- 
dow zones will narrow them down. Such zones fall 
into two categories: tangential and cast. 


n i = 


i 

* N 


y z 

< 7 a 

y z 

r t 


I '* = ± N 


Z X 

a a 


Z X 
T T 


n 3 = ± 


N 


x y 
a a 

x y 

T T 


( 15 ) 


where the subscripts denote partial derivation with 
respect to the pertaining variable, and 


N = 


r 

y z 

2 

y x 

2 

x y 

2 


a a 

+ 

a a 

+ 

a a 



y 


z X 


x y 


V 

T T 


r r 


T J T 



( 16 ) 


Tangential incidence means that ot - 90°. Rela- 
tion (17) tells that the tangential shadow boundary, 
s^ , is given by the equation 

n . s = 0. (21) 

The cast shadow cylinder will be rooted on a line 
on the surface, either a sharp rim or a tangential 
shadow boundary (as, e.g. , of an airplane fuselage 
casting a shadow on the wings ) . Let the line be given 
in the parametric form 

x c = f(£), y c = g(£), z c = h(£) . (22) 


Once the surface is known, it is an easy matter to de- 
cide on the sign pertinent to the interior normal. 

This vector and the unit vector 

S = iSi + js 2 + ks 3 

in light direction combine to define the angle of inci- 
dence, O', by their interior product: 

n • s = cos a . (17) 

The elemental surface dS, a second requisite, 
is proportional to N: 

dS = N d(7 dT . (18) 

The differentials da and dT should have the same sign 
here, so that both variables must be taken increasing 
(or decreasing) at the later integrations. 

B y expressions (11), (12), (17), and (18) the 
second order force differential acting on a perfectly 
absorbing elemental surface may be written as 


The relevant generating lines of the shadow cylinder, 


x-x y -y z-z 

c _ c __ c 

Si S 2 s 3 


(23) 


start at the line segment (or segments) grazed by 
rays. They may intersect with other portions of the 
body forming a cast shadow boundary, s^, also to be 

traced on the (a,T ) -diagram; it leads to further re- 
striction of the integration interval, unless it lies 
wholly within the tangential shadow zone. 

With bodies of complex description, the deter- 
mination and use of proper integration limits may re- 
quire painstaking detail work. Often the integrations 
will have to be numerical. 

To know the force alone is insufficient informa- 
tion. Its point of attack (centroid) is also needed. 

To find it, one determines the arm of the force con- 
necting a point, C, on the line of attack with a point 
of reference that may be chosen arbitrarily and will 
here be taken as the origin O. 


d 2 F = p (n • s) N da dr s . (19) 

-a o — — ' ~ 

(The light direction is now s_ rather than k) . 


Let r denote the vector leading from the origin 
to an elemental surface. The force attacking there 
has the second order moment 
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d 2 M* = r x d 2 F , 

where F stands for either F or F a , The moment 
— -a — t 

of the integrated force is equal to the integrated mo- 
ment: 

r* x F = J f [r xd 2 F] (24) 

a r 

where 

r* = ix* + jy* + kz* (25) 

is the arm of the force F and as such defines 
C(x* , y* , z* ) . The integration limits above apply 
again. Equation (24) for the vector r* leaves one 
of its components indeterminate, since a force can 
be moved along its line of attack without altering its 
moment. One may place the centroid at a convenient 
location on this line by introducing a suitable rela- 
tionship 

f(x* , y* , z*) = 0 . (26) 


of the force F is not permanently zero, the vehicle 
will begin turning about shifting axes through G, so 
that F will be varying. This motion can best be de- 
scribed when the axes of the principal moments of 
inertia intersecting at G are introduced as coordinate 
axes. If vector components in this second body -fixed 
system, 2 , are identified by the indices a, (3, y, the 
torque (27) and the vector co of the angular velocity 
which it causes to develop are related by the Eulerian 
equations 


T 

d 

to 

- m 

to 

<T U 

-T ) 

a 

dt 

a 

P 

y 

b 

c 

T 

d 

to 

- to 

to 

(T 

- T ) 

b 

dt 

(3 

y 

0> 

c 

a 

T 

d 

CO 

- CO 

OK 

(T 

- T, ) 

c 

dt 

r 

a 

P 

a 

b 


(28) 


where T , T , and T are the principal moments of 
a b c 

inertia. The torque components, M^, M^, fol- 

low from M , M and M by the transformation equa- 
x y z 

tions connecting the two body -fixed systems. 


IE. REACTION OF A BODY TO THE IMPACT 
OF LIGHT 


For a limited time, an undisturbed satellite in 
Earth orbit will move at practically constant velocity 
along practically rectilinear path. A (|, 17 , £) -sys- 
tem through its mass center, G, will be equivalent 
to the original space -fixed system by the Galileo 
principle. Let us assume that at t = 0, the vehicle is 
at rest in ( £, 77 , t) • 

The surface equation ( 14) has been set up in a 
body -fixed (x,y, z) -system ; the differentials (19) 
and (20) of the light force F must be integrated in it. 
The light direction s_ is constant in (£, 77 , £) , or 
can safely be taken so. In (x,y,z) , however, it is 
different with different attitudes of the body. More- 
over, the shadow boundaries, and therefore the inte- 
gration limits, will often be different. Each attitude 
requires a separate force computation. 

Let the mass center G have the directed dis- 
tance 

r = xi + yj + zk 

from the originO of the (x,y,z) -system. If the 
torque 

M = (r* - r) xF = Mi + Mj + Mk (27) 

— - x y z 


The integration of the set (28) can be carried out 
independently of the translatory motion of the mass 
center, yet it remains laborious nonetheless. ( 1 ) 
Finding the angular velocity co prevailing at the end of 
a step is not enough. A second set of equations re- 
lating the co -components to the first order time- 
derivatives of the Eulerian angles must be solved to 
obtain the instantaneous values of these angles and 
thus the body's attitude in ( £,rj,£) at the beginning 
of the next step. Then one has (2) to determine 
s(x,y ,z) , ( 3 ) to use s for finding the shadow bound- 

aries and integration limits, (4) to carry out the 
force quadrature, (5) to establish the instantaneous 
centroid (r*) , ( 6 ) to compute the vector M, and (7) 
to compute its components in the system 2 . After 
this, one is ready to proceed by At in the step calcu- 
lation as it was set up. 

Since control will not permit major deviations 
from an original attitude, the working out of this 
scheme does not seem to be warranted. It will suf- 
fice to consider an initial period in which the system 
2 barely is shifting and which might be quite long 
since p is extremely small. The vectors F and M 

will then essentially retain the values they had at 
t = 0. The subtractive terms in the set (28) , as 
small of second order, can be neglected; the solution 
follows as 
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where i f , j f , k 1 are unit vectors on the axes of the 
principal moments of inertia. This expression de- 
fines the direction of the initial axis of the rotation 
about G and shows that co increases linearly with time 
during the initial period. 

Untouched by the rotation process as such, the 
mass center G will move relative to the ( £ , 77 , f ) - 
system. Newton's translation equation 



(m: body mass; s: distance covered in force direc- 
tion) cannot be solved, however, without regr ess to 
the rotatory equations which are to provide the dis- 
crete values of F needed for the necessarily stepwise 
integration. With F = const, in the initial period 
the translational velocity, like the rotational, is a 
linear function of time. 


IV. PRESSURE FORCES ON A RIGHT 
CIRCULAR CYLINDER 


The cylinder , as an approximation to an orbiting 
telescope, may be considered a fitting body to illus- 
trate the computational methods. It is assumed that 
its surface material is either totally absorbing or 
totally reflecting and that the lens end always points 
away from the sun. If the top of the cylinder ( Fig. 

2 ) is fitted with the optical objective, the light will 
have to arrive from below ( S 3 > 0 ) . 

The curved part of the cylinder surface in Figure 
2 is described by the system 

x = R cos r 

y = R sin t (30) 

z = 2 La 

where the natural ranges of the variables are 
-i <(r<+l 
0s t S2tt . 

From expressions (15) and (16) 
n = -i cos r - j sin r (31) 

N = LR = const. 



Condition (21) then gives the tangential shadow bound- 
ary s^ at 

s l 

tan r = - — . (32) 

t S 2 

If the first of the solutions is t* , the second will 

* 

be T t= T - 7 T . Thus, the integration intervals are 

t* - 7 r St St* , -l^o^l (unaffected). 

There is evidently no shadow cast by either top or 
bottom rim on the otherwise illuminated surface half. 
The (double) line s^ is the only shadew boundary. 

With perfectly absorbing surface material, the 
differential ( 19) , on integrating over the proper in- 
tervals, gives 

F = 2p LR J 1 - s 3 2 s s C s , 

-a,c o J - a - 

where relation (32) has been used to convert trigono- 
metric functions of t* into expressions involving Sj 
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and s 2 . The index c indicates integration over the 
curved part of the surface. The force on the blunt 
end is 

F =p 7r R 2 s 3 s = B s . 

-a,b *o 6 - a - 

This result can be obtained with ease from d 2 F 

- a,b 

p^ (k ■ s) dS, the interior normal being n = k. 


The line of attack of F is normal to the cylinder 
— t,c 

axis and will pass through O for symmetry reasons.* 
That of F t k is the cylinder axis itself. Thus, the 

centroid of the combined force again can be made to 
coincide withO; as above, the torque is zero with a 
homogeneous cylinder. The force direction, however, 
is no longer that of the incident light. 


The elemental forces d 2 F have the arms 
- a,c 


r = ix + iy + kz 


relative to the origin 0, where x, y, z are the coor- 
dinates (30). By expressions (19) and (31) the in- 
tegrated moment becomes 


+ 1 r 

— e ’ -»o “ / / 

<J = ~ 1 T* - IT 


l j k 

R cos r R sin t <t 
s, s 2 s 3 


{s t cos t + Sj sin t) do dr 


= - f P Q lr2s s < s 2 i - »i J) - 


The blunt end force F . obviously attacks at the 
“a,b 

base circle's center (r = - \ Lk ) ; its moment is 
i j k 

-b = \ P o * r2 83 0 0 -1 

s 2 s 2 S3 


O lr2 s 3 (s 2 1 - Si j). 


Since the forces have the combined moment zero, the 
arm of their resultant is r 5 ^ = 0.* The centroid is at 
the origin (C = O) , i.e. , at the geometric midpoint 
of the cylinder, which, with a homogeneous mass 
distribution, is the mass center as well. A solid 
(right circular) cylinder will experience no torque 
under parallel radiation if its surface is perfectly 
black. 


If the surface is perfectly reflecting, integration 
of the differential (20) leads to 


F 

“ t,c 


P Q LR n/i - S 3 2 (is t + js 2 ) = C t (is! + js 2 ) , 


while the blunt end force becomes 

F = 2 p 7T R 2 k . 
t,b o 


Suppose now that the cylinder is subdivided into 
compartments of rotational symmetry with respect to 
its axis, each filled with a homogeneous material dif- 
ferent from that in the adjacent compartments. The 

T -axis will then be the z-axis. Since T_ = T here, 
c b a 

the remaining two principal axes can be taken as paral- 
lel to the x- and y-axes. The primed unit vectors are 
identical with i, j, k in these circumstances. Hence, 
the components of any vector have the same magnitude 
in the two body -fixed systems , so that the indices a, 

13, y and x, y, z are interchangeable. The mass cen- 
ter, G, will be located on the z-axis at some distance, 
z = z , from the origin (and centroid; O = C). Thus, 

A ° * 

r = z^k, while r = 0. Expression (27) for the 
acting torque becomes 

M = -z [k x FI = z (iF - jF ) = iM + jlVl + 0. 

— o 1 o y J x' a J (3 

The angular velocity (29) takes the form 

2 =#< iF y-5 F X )> 

a 

which result holds for the initial period only as will be 
recalled. The symbol F is used either for F^ or for 

F . The x- and y -components of both these forces are 

proportional to Sj and s 2 . Thus, we may write 
z z 

S2 = A (is 2 - js t ) = A[s x k] , 
a a 

where A stands for C + B or C , respectively. The 
. a a t 

initial axis of rotation is shown to be the normal through 
the mass center to the plane containing light direction 
and cylinder axis. From this position and attitude, it 
will slowly wander out in subsequent stages following 
the complex pattern set by the Euler ian equations. 


* 

This can also be concluded from equation ( 24) where 
The use of equation (24) is not needed here. the right side can be shown to vanish. 
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The magnitude of the angular velocity , at any small 
t, is zero if S3 =* 1 (s = k; light directly from below) ; 
otherwise it is the ratio 

C t 4 L*il - S3 2 

C + B 3 > ■ ’ 

a a L V 1 - s 3 2 + | R s 3 

which is largest (4/3) for s 3 = 0. Incidence normal 
to the cylinder axis causes somewhat faster initial 
rotation with a perfectly reflecting than with a per - 
fectly absorbing surface. The ratio, however, de- 
creases as s 3 increases, and soon the absorbing cyl- 
inder will exhibit the faster rotation. This trend re- 
flects the fact that the component s 3 , totally lost in 
curved-surface reflection, but effective in absorption, 
increases at the expense of s t and s 2 (s t 2 + s 2 2 = 1 - 
S3 2 ). 

V. CONCLUSION 


is used to define a "light force coefficient" (C ) by 
putting F 

n =C_p . 

F o 

This coefficient moves between C = 0 (ideal trans- 

F 

mission) and C F = 2 (ideal reflection at a = 0°) . It 

represents one way of introducing a correction to the 
basic pressure p^. Since p ! Q and p"^ depend on a = 

arc cos (n • s) , the force quadrature is rather 
laborious even with simple surface and favorable 
incidence [41. 

Other, and more fundamental, difficulties arise 
with light incident on rough, scattering materials or 
on metallic surfaces emitting electrons when struck. 
In both cases, the theoretical derivation of the force 
differential must be re-examined and will have to 
proceed on lines deviating in places from those 
followed so far. 


The force, its point of attack, and the initial mo- 
tion under its impact can be determined, in basically 
the same way, for all surfaces whose parts are given 
in analytic terms, although the task is bound to grow 
wearisome with intricate shadow contours or multiple 
reflections in concave recesses. 

A refinement would consider the losses in par- 
tially transparent and absorbing materials as far as 
experimental or theoretical data are available. In 
such cases , the total force element has to be written 
as 

d 2 F = d 2 F (l) + d 2 F (r) = (p' k + p" k') cos a dS, 
t ““ — o o 

where both p 1 and p" are smaller than p . The 
o o o 

compact expression ( 13) is no longer obtainable then. 

In [4] , the absolute value 


II = !p T k + p n k’ | = v p |2 + p ,t2 + 2p T p" cos 2a 
o o o o o o 
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VII. STRUCTURAL DYNAMICS 



A DYNAMIC APPROACH OF FUEL SLOSHING PROBLEMS IN A SPACE VEHICLE 


By 

Frank C. Liu 


SUMMARY 


A dynamic approach to determination of the 
natural frequency of fuel sloshing in a space vehicle 
in flight is developed. This method enables us to 
couple directly the motion of the space vehicle re- 
sulting from fuel sloshing to the motion of the fuel. 
Hence, it yields a more realistic representation of 
the dynamic system than that given by the existing 
approach which links the motion of the vehicle and 
its fuel sloshing by means of a mechanical analogue 
or feedback. 

L INTRODUCTION 

The force used to maneuver a space vehicle in 
flight is usually an impulse. To determine how the 
space vehicle responds to an impulse force, it is 
important to know precisely the natural frequency 
of the whole system as well as the dynamic coupling 
of the motion of the vehicle and the motion of the fuel 
in its container or containers. The existing approaches 
to this problem are as follows: 

1. Mechanical Analogy. First, we analyze the 
motion of the fuel in the container, which rests on a 
fixed support, and determine the natural frequency 
of the sloshing. Then, the fuel tank system is con- 
verted to a mass-spring system by means of a me- 
chanical analog to constitute a dynamic component of 
the space vehicle. 

2. Feedback Approach. We impose a motion to 
the fuel tank and analyze the response of the fuel. 

Then, we determine the force created by the fuel 
sloshing and feedback to the vehicle to determine the 
response of the vehicle. 

It is easy to see that neither of the two methods 
couples the space vehicle and the fuel sloshing di- 
rectly. 

IL THE DYNAMIC APPROACH 

This new approach of solving sloshing problems 
of a space vehicle treats the vehicle and its fuel as 


one unit. First, let us show two examples to illus- 
trate how the frequency of a simple pendulum is 
affected when it is attached to a free -moving support. 
We will have some idea about how the natural sloshing 
frequency of the fuel tank is affected by its free 
moving support. 

Example 1. Let us attach the simple pendulum 
with mass m and length L to a mass M which is free 
to move horizontally in the plane of oscillation, as 
shown in Figure 1. 

The natural frequency of the system can be de- 
termined easily by using Lagrange*s method. 

w- J (i + (g/L) rad. /sec 

This shows that the frequency is increased by a factor 
\/ 1 + m/M . 



FIGURE 1. DIAGRAM OF EXAMPLE 1. 

Example 2. To increase the freedom of the 
support further, let us make the support also tree to 
rotate about its static center of gravity, as shown in 
Figure 2. Notice that a mass m is put at the opposite 
side of the axis of rotation to give a static balance of 
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FIGURE 2. DIAGRAM OF EXAMPLE 2. 

the system. Here M denotes the total mass of the 
support (stand, disc, and counter balance mass m), 

J is the moment of inertia of the disc with the balance 
mass m included, and R is the distance from the 
point where the pendulum is attached to the center of 
rotation. Now, the frequency of the system becomes 



FIGURE 3. SPACE VEHICLE HAVING 
PARTIALLY- FILLED FUEL TANK. 


Let us return to the space vehicle which has a 
partially filled fuel tank as shown in Figure 3. Because 
of its axial symmetry, only the motions x and © are 
to be considered. During space flight the surface 
wave of the fuel becomes the only restoring force for 
the free vibration of the system. Usually the mass of 
the fuel forms a major portion of the total mass of a 
space vehicle. Thus, we may expect that the fre- 
quency of the vehicle is considerably higher than the 
frequency of sloshing in a fixed tank. 

m. GENERAL PRINCIPLE 

In an attempt to determine the velocity potential 
function of the fluid which is assumed incompressible 
and inviscous in a rigid, free moving cylindrical 
tank we found that under the linearized boundary 
conditions the coefficients of the potential function 
are directly related to x and £ . For small oscilla- 
tions we may assume that x= X t sin cot and 6 = X 2 
sin cot and use a one-term approximation for the 
potential function. Formulating the kinetic energy 
of the dry vehicle in terms of x and 5 f and the kinetic 
and potential energies of the fluid in terms of the 
potential function results in the Lagrangian function 
of the system in the form 


I = |x T wx 

where X is a column matrix with two elements X t 
and X 2 , T denotes matrix transposition, and W is a 
2x2 matrix. The elements of W are polynomials of 
co, the natural frequency of the whole system. By 
Hamilton's Principle, I is an extremum when the 
determinant of W is equal to zero. Consequently, 
we can solve for co numerically. 


IV. CONCLUSION 

This method has addea advantage when it is used 
for a system which has more than one fuel tank. The 
dynamic coupling of all the tanks and the vehicle be- 
comes a simple matter. 

Digital computer programs for a vehicle with 
cylindrical tanks have been completed. The frequency 
of a vehicle of the Saturn class with one tank sloshing 
may be 5% higher than the fixed tank frequency. Work 
on the digital computer program for tanks with ar- 
bitrary shape is in progress. A technical report on 
this work will be available shortly. 
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WIND PENETRATION EFFECTS ON FLIGHT SIMULATIONS 


By 

James G. Papadopoulos 


SUMMARY 


Symbol 


Definition 


Present six-degrees-of-freedom wind monitor 
Saturn Simulation programs assume an instantaneous 
constant wind velocity over the entire vehicle length. 
Contrary to this popular but simplified assumption, 
the study presented here indicates that wind penetra- 
tion effects are not only desirable but necessary for 
any serious analysis of vehicle response to a wind 
profile. 

Wind penetration effects are evaluated in the 
frequency domain by developing a solution for theo- 
retical transfer functions. Numerical calculations are 
made for a Saturn SA-9 vehicle. Analytical transfer 
functions for angle -of -attack and average engine 
gimbal deflection are compared to a spectral analysis 
of actual flight test data. 


C 0 In-phase component of aerodynamic 

moment due to a sinusoidal unit wind, 
Newton -sec 

C t Out-of-phase component of aerodynamic 

moment due to a sinusoidal unit wind, 
Newton -sec 

C 2 Unit gimbal deflection restoring mo- 

ment divided by vehicle’s mass moment 
of inertia, i/sec 2 

Aerodynamic lift gradient distribution 

D 0 In-phase component of generalized aero- 

dynamic force due to a sinusoidal unit 
wind, Newton -sec /meter 



LIST OF SYMBOLS 

Symbol 

Definition 

a 0 

Attitude control gain 

a l 

Angular rate control gain, sec 

a l 

Longitudinal acceleration, m/sec 


Dj Out-of-phase component of generalized 

aerodynamic force due to a sinusoidal 
unit wind, Newton -sec /meter 

D t Aerodynamic damping, 1/sec 

F Vehicle thrust, Newtons 

F AERO aerod y namic f° rce acting on the 

vehicle, Newtons 


A 


m 


Generalized net aerodynamic force at 
the m bending mode , Newtons 


F AWP Aerodynamic force due to wind pene- 

tration, Newtons 


A 


mwp 


A 


mfd 


B 0 


Bi 


Generalized aerodynamic force due to 
wind penetration for the m bending mode, 
Newtons 

Generalized force term, Newtons 

In-phase component of aerodynamic 
force due to a sinusoidal unit wind , 
Newton -sec/meter 

Out-of-phase component of aerodynamic 
force due to a sinusoidal unit wind, 
Newton-sec /meter 




g2 


G a 


G 




G 


Viscous damping of the (jl bending 
mode 

Control accelerometer gain, rad -sec 2 /m 

Gimbal actuator transfer function 

Attitude control filter transfer function 

Attitude rate control filter transfer 
function 

Control accelerometer transfer function 
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LIST OF SYMBOLS (Cont*d) 


LIST OF SYMBOLS (Cont'd) 


Definition 

Aerodynamic force term, Newtons 

Aerodynamic moment term, Newton - 
meter 

Aerodynamic force term, Newtons 
Aerodynamic moment term, Newton - 


Complex notation for 


Mass moment of Vehicle about center 
of gravity, Kg-m 2 


Mass moment of inertia of gimbal engine 
about its center of gravity, Kg-m 2 

Aerodynamic force term, Newton/rad 


Aerodynamic moment term, Newton - 
meters/rad 

Aerodynamic moment term, Newton - 
meter 

m generalized bending force term due 
to the n bending mode. Newton/meter 

Distance between gimbal point and en- 
gine center of gravity, meters 

Vehicle length, meters 

La Place operator 

Vehicle mass, kg 

Net aerodynamic moment acting at the 
vehicle's center of gravity, Newton- 


Aerodynamic moment at the vehicle 
c.g. due to wind penetration, Newton- 
meter 

Mass of gimbal engine , kg 

Generalized mass of the p bending 
mode, kg 

Slosh mass of ith propellant, kg 


Symbol 


Definition 


Dynamic pressure, Newtons/m 2 

The real part of a complex number 

Complex notation for iu>, rad/sec 

Reference surface area, meters 2 

Vehicle velocity, m/sec 

Wind velocity along vehicle, m/see 

Amplitude of wind penetration, m/ sec 

Displacement along vehicle, meters 

Longitudinal center of gravity, meters 

Longitudinal location of attitude sensor 

Longitudinal location of gimbal engine 
center of gravity, meters 

Longitudinal location of the ith slosh 
mass, meters 

Normalized p bending mode displace- 
ment along the vehicle 

Normalized ^ bending mode displace- 
ment at the gimbal station 

Normalized p bending mode displace- 
ment at the angle -of -attack sensor 

Normalized g bending mode displacement 
at the accelerometer sensor 


Normalized /x bending mode displace- 
ment at the attitude sensor 

Normalized \x bending mode displace- 
ment at the angular rate sensor 

Normalized /x bending mode slope along 
the vehicle, i /meter 

Normalized n bending mode slope at the 
gimbal station , 1 /meter 

Normalized g bending mode slope at the 
angle-of-attack sensor , 1/meter 

Normalized p bending mode slope at the 
accelerometer sensor, i/meter 
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LIST OF SYMBOLS (Cont’d) 


LIST OF SYMBOLS (Concluded) 


Symbol 

Definition 

Symbol 

Definition 

y 1 

\x<P 

Normalized / u bending mode slope at the 
attitude sensor , 1/meter 

<P 

Angular rate of undeformed missile axis 
rad/sec 

y W 

Normalized \x bending mode slope at the 
angular rate sensor, 1/meter 


Indicated angular rate at the rate gyro 
location, rad/ sec 

y Si 

Normalized / \x bending mode displace- 
ment at the ith slosh mass 

<p 

Phase angle of wind penetration with 
respect to the vehicle, radians 

y, Si 

Normalized ^ bending mode slope at the 
ith slosh mass, 1 /meter 

B/u 

Frequency of the /u bending mode, 
rad/ sec 

y 

Lateral rigid body velocity at its center 
of gravity, meter/sec 


I. INTRODUCTION 


y 


O'. 

1 


a 

x 



£2 

£3 

£4 

£5 

<P 


Lateral rigid body acceleration at its 
center of gravity, meter/sec 2 

Indicated angle -of -attack, radians 

Local angle -of -attack, radian 

Engine gimbal deflection, radians 

Engine gimbal command, radians 

Acceleration of the control acceler- 
ometer, m/sec 2 

Generalized bending coordinate of the n 
bending mode, meters 

Lateral displacement of outer fuel slosh 
mass, meters 


As a vehicle flies through a wind profile, resul- 
tant local wind components along the vehicle create 
changes in local angle of attack. This effect is re- 
ferred to as wind penetration. This report will de- 
termine the relative importance of wind penetration 
effects in flight simulation programs by developing 
an analytical method of evaluating theoretical trans- 
fer functions. 


The expression for the local angle of attack, of , 
of a flexible vehicle with a local wind component, 

V , at a particular vehicle station, x, is 



Lateral displacement of outer LOX 
slosh mass, meters 

Lateral displacement of center LOX 
slosh mass , meters 

Lateral displacement of upper LOX 
slosh mass, meters 

Lateral displacement of LH 2 slosh mass, 
meters 

Angle between undeformed missile axis 
and space fixed reference, radians 


The net aerodynamic force acting on the vehicle is 
obtained by summing the product of the local lift 
gradient, dynamic pressure, reference area and lo- 
cal angle of attack. 

L 

F a = q S f C, a dx 

AERO J W x 


9 s ['”- v] / J' C lx V WX dX 


+ V*/ °lx ,,! CG- X)dX 'V2% 

o 1 o 


Indicated vehicle attitude, radians 



( 2 ) 
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The net aerodynamic moment acting at the vehicle’s 
center of gravity is 

T ’ 1 L 

M AERO - qS l*’ - vJ / Vco' x) dx * 


* ^ J' C 1 *'W x CG- x,dx 


* v’ / °lx |x CG _x,!dx 


as 


/ C lxV (X CG- X) 

[L-l o 


dx 


/ C lx y V X CG- X) 

n=i o 


(3) 


The generalized aerodynamic bending force for the 
(fi=m) bending mode is 

A = qS \<p - 1 f C y dx 

m L V J J lx mx 


^ f C V y dx 
V J lx wxV 
o 

L 

C lx (X CG~ X) y mx ^ 


as 

v 


qS 

V ^ "n 
n=t 


11 1-i 

Jj T) f c y Y dx 
J lx jjx mx 


o 

L 


- qS Yj i 7 f C_ y ’ y dx. 
U A a J lx ux mx 
M~1 o 


(4) 


These relationships utilizing local wind distribution 
to define local angle of attack, aerodynamic forces, 


aerodynamic moments and generalized bending forces 
(equations 1-4) will be the means of inserting the 
wind penetration effects into the development of the 
theoretical transfer functions. 


H. THEORETICAL TRANSFER FUNCTION 
ANALYSIS OF WIND PENETRATIONS 


Theoretical transfer functions are analyzed as 
sinusoidal amplitude and phase response of flight 
mechanic parameters with respect to a given sinu- 
soidal input. In this study , a sinusoidal wind profile 
is used as the forcing excitation. Transfer functions 
between any two flight mechanic parameters can be 
obtained from transfer functions of a common base by 
dividing amplitudes and differencing phase angles. A 
numerical evaluation of transfer functions was made 
for the SA-9 vehicle. 

The SA-9 control system employed a fully active 
ST -124 platform which generated the attitude error 
signals, a body-fixed control accelerometer located in 
the Instrument Unit which provided partial load re- 
lief, and a three -axis control rate gyro package for 
angular rate control also located in the Instrument 
Unit. Angle -of -attack components were measured by 
a mode F 16 Q-ball angle -of -attack transducer mount- 
ed on the top of the Launch Escape System. The en- 
gine gimbal deflections were recorded from the indi- 
vidual actuator positions. 

The mathematical technique of computing trans- 
fer functions is to determine the La Place Transform 
of the equations of motion and to evaluate the variables 
as a function of frequency. A total of fourteen equa- 
tions (5-18) were used to represent the vehicle’s dy- 
namics and control in the yaw plane. 
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rotation equation 

a l<p 

<P + D x <p - —j— 


L 4 | *vv | 3 M si^j * 


M 


AERO . „ „ . F rj y — 
~T~ + ° 2/3 + I (y E _X E y E } 


M ^ x + 1 
E E E E 


+ T | VV*| **lV*T | 3 VVi +A i- V 


(5) 


translation equation 
F 


y - A L<P- M 2M'"E M 


M L 8 F v’ 

+ " 


iVERO F „ ”E E^E . y " J M E _ i V M V - — f M f = 0 ( 6 ) 

L M T) ,, M Li A1 S1 


M n n M Sj^j M 

H j=l J = 3 J 1 


angle -of -attack sensor equation 
1 


-cp + a + 


V 


y "" V W a (X CG- X a )cp +T) V 


+ v y T = 0 

pa 


(7) 


control equation 


p - g 2 G*» y - a G <jP-aiG*<]P — 0 
p c y o <p <P 


( 8 ) 


engine actuator equation 


s - g a p =o 

P E A p c 


(9) 


control accelerometer equation 


(x <P + A_ <P) - Z + y - Y* (y •• r\ + A y f •• r\ ) = 0 

T , 1 LJ ^ //v 'ii T , //v 'n 


py P L |iy > 


( 10 ) 


attitude sensor equation 

<p. - <p + y.y f rj =0 
i p(p V 


„ rj =0 
pep 'p 


attitude rate sensor equation 

ip. - <P + Tj y' 

1 ^ b 

P 

generalized bending equation (mth mode) 

FY, A 

E • 2 v m 

j3 + (r? + a; g n+w n ) + rr - 

2M_ P E ' m Bm g Bm 7 m Bm % ; M 


Bm 


__ M 

Bm Bm 


(M Sj y SJ*j + A L y, Sj l J l 


^ Z ^Sj y Sj + A L y *Sj ^ ~ y mE + 'e y 'mE^*E 

Bm j=3 


slosh equations 


(11) 


( 12 ) 


(13) 


- x si* + a l^ + 2 

+ V 

” A L y 'si ^ 

+ <*i + w si g siq + w Si v = 0 

(14) 

~V + A L ^ + Z 

+ y S2 ^ 

- A L y S2 TJ 

+ ( ^2 + W S2 g S2^2 + “S2*2> = ° 

(15) 
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(16) 


-x„„ <p 
S3 


' X S4* 


-v 


+ A l * + 2 + V ' 

+ A L «P + *z + T) 

+ A L + z + y g5 v 


A L y, S3 TJ + ( ^3 + W S3 g S3^3 + W S3 ^3* 
A L y S4 17 + ( '^4 + W S4 g S4 ^4 + W S4 ^4 * 
“ A L y, S5 TJ + ( ^5 + “SS g S5 ^5 + W S5 *5* 


0 

0 

0 


(17) 

(18) 


The vehicle was assumed elastic and was repre - 
sented by the first bending mode. The resulting re- 
sponse should be fairly accurate up to approximately 
3. 0 cps which is the upper limit of the Flight Test 
Values. The cluster of propellant tanks in the S-I 
stage was assumed to have the same bending deflec- 
tions as the center LOX tank. This is a good as- 
sumption for the first two bending modes, particu- 
larly since the bending deflections of all the S-I pro- 
pellants are nearly identical and small with respect 
to the upper stage deflections. These cluster effects 
on transfer functions would be significant at frequen- 
cies greater than 5.0 cps. Quasi-steady aerodynam- 
ics, which states that the force along the vehicle 
depends only on the instantaneous local angle of 
attack, was assumed. The amplitude and phase of 
each control filter (accelerometer, attitude error, 
and angular rate) were continuously evaluated at the 
wind excitation frequency. 


The La Place transform for most of the equations 
can easily be written. However, the transform of 
terms involving the integral of the wind penetration 
distribution over the vehicle require special analysis. 
These terms occur in the expression for the aerody- 
namic force, the aerodynamic moment, and the aero- 
dynamic component of generalized force. Consider, 
for example, the wind penetration term in 

the aerodynamic force equation: 


F = SL§ f c 

AWP V J lx 
o 


V„ r * • 
Wx 


( 2 ) 


At any given vehicle flight condition, this term is de- 
pendent in frequency, amplitude and phase with the 
exciting wind oscillations. 


The transform of the wind penetration integral 
was numerically evaluated for particular flight time 
conditions. A common vehicle reference, station 0, 
was established to define the sinusoidal wind distri- 
bution over the vehicle (Figure 1) . 



FIGURE 1. REFERENCE WIND PENETRATION 
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Equation 19 indicates that the aerodynamic force due 
to wind penetration varies sinusoidally with the phase 
angle ip . In-phase and out -of -phase components for 
La Place Transform are referenced as ip equals 0° 
and 90®, respectively. Thus,at a particular frequency 
if 

L 


B ° = ¥ / c i sin (f) * 

o x 7 

B - = ^ {£■ 8ln (v + t)J 


dx 


then 


L(F awp ) = <B ° + Bli) L( V- 


(20) 

( 21 ) 

( 22 ) 


In a similar manner, the wind penetration component 
of the aerodynamic moment, equation (3) , can be ex- 
pressed as 

M ATTyT . V / °1 V WX (X CG~ X) ^ 


AWP 

= y r 

v v wJ 

c 

letting 


\ sin {f + *) 


( X CG " X ) dx (23) 


C 0 = 


V 


(X CG' X) dx 


C, = 


_ gs 

V 


/ c , - (f ) 

C 1 sln (f 4 fj 


(24) 


/ 


then 


L(M AWP ) = ( c 0 + c l i) L(V W ). 


(x CG -x)dx (25) 


(26) 


Also, the wind penetration component of the general- 
ized aerodynamic bending force for the mth bending 


then 


L<A mwp> ' < D . 4 °.‘> L 'V 


(30) 


The transform coefficients B 0 , B t , Cq# c i» d o 
and Dj can be expressed as polynomials in co for any 
particular flight time condition by performing least 
square solutions of results from specific values of 

CO . 

n 


Bo = E B 


k=0 

n 


Ok 


B t = E B 


k=0 

n 


lk 


c o - E c ok 

k=0 0k 
n 

Cl = E c 

k=0 

D o = E D ok * 

k=0 
n 


lk 


Dj = E D 


k=0 


lk 


co 


(31) 


For ease of presentation, definite integrals are 
represented by the following set of constants: 


J 0 = qS ( C dx 

IX 


3, ■ «S f C b (x CG -x) 


dx 


mode equation ( 4) is 

A = ^ f C, V„_, Y dx 

mwp V J 1 WX mx 

0 



J 2 

= qS 

! C lx (X CG- X,!dX 

o 

L 


y dx . 

mx 

(27) 

G 

= qS 

f C. y dx 

J 1 X fl 

o 

L 

Letting 

- / \ 



G 

V 

= qS 

f C lxV*CG- X,dX 

o 

D o = \r~ J C sin (— j y dx 

u V J 1 \V / mx 

o x ' 

(28) 

H 

= qS 

L 

f C. y f dx 
J lx // 

0 

D i = ¥ / c i sin (f + i) : 

V dx 

mx * 

(29) 

h" 

= qS 

/ C lx J '„ ,X CG- X,dX 
0 
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u 

A = qS f C. y y m dx 
mix J lx (i m 

o 

L 

K = qS f C. y f y dx 
m/i M J lx |i m 


(32) 


The La Place Transform of the yaw flight me- 
chanic equations (5-18) can now be expressed as the 
following set of simultaneous linear equations in the 
S domain: 
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(x^ S 2 + A) L(<p) - S L(Z) + L(y) - (y~ S 2 + Ay'”) L (77) = 0 
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These equations rewritten in matrix notation are 



(47) 


Although both matrix "A" and column ”B” are complex 
coefficients of the transform elements, the solution of 
these simultaneous linear equations is generally well 
known. One approach would be to expand the deter- 
minants as an S polynomial according to Cramer’s 
rule. For example, 

L (<p) = [function of S] L (V„_) 

w 

or 

= function of S . (48) 

L(V W ) 

If this were evaluated at a particular frequency , the 
transfer function would result in a complex number , 
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The transfer function is more conveniently expressed 
as a harmonic relationship at any particular fre- 
quency with amplitude ratio and phase angle. 


Amplitude ratio = •si R 2 + I 2 


Phase angle 



(50) 

(51) 


III. ANALYTICAL RESULTS AND COMPARISONS 
WITH RECORDED FLIGHT TEST DATA 


The Flight Mechanic Transfer Functions includ- 
ing the wind penetration effects as indicated by the 
equations were evaluated for the SA-9 80 -second 
flight condition. The resulting analytical transfer 
functions for gimbal engine deflection and angle of 
attack are shown as solid lines on Figures 2 and 3 
respectively. For comparison, corresponding trans- 
fer functions were evaluated for a wind distribution 
instantaneously constant over the entire vehicle - 
that is, no wind penetration. These transfer functions 
are indicated as dashed lines in Figures 2 and 3. For 
this particular condition, the wind penetration effects 
are more pronounced in the gimbal engine transfer 
function than that of the angle of attack. In fact, the 
wind penetration effects for the gimbal engine deflec- 
tion are five times greater in amplitude at 1.7 cps 
than corresponding analytical value without wind 
penetration effect. 

Further inspection of the gimbal engine deflection 
and angle -of -attack transfer functions (Figures 2 and 
3) indicates that the wind penetration effects are 
strongly coupled to the first bending mode (peak at 
2. 3 cps) . A transfer function plot of the first bend- 
ing mode generalized coordinate, 171 , shows this to 
be true ( Fig. 4) . The wind penetration effect in- 
creases the amplitude of the first bending generalized 
coordinate transfer function by a factor of 2.6 at 2.3 cps 
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and 6.6 at 4.0 cps. This, in effect, states that ne- 
glecting wind penetration effects results in under- 
estimating the bending mode coupling by a consider - 




FIGURE 3. ANGLE-OF- ATTACK TO WIND 
EXCITATION TRANSFER FUNCTION 


FIGURE 4. WIND PENETRATION EFFECTS ON 
GENERALIZED BENDING MODE TRANSFER 
FUNCTION 

A spectral analysis of recorded SA-9 flight test 
data was made hopefully to corroborate the wind pene- 
tration effects. The corresponding amplitude transfer 
functions are shown plotted on Figures 2 and 3. Even 
with the uncertain amount of noise content associated 
with flight telemetry data, the results are surpris- 
ingly close. The comparison of engine gimbal de- 
flection indicates the flight test results to be much 
closer to the analytical values which include the wind 
penetration effects. The angle-of -attack transfer 
function comparison indicates the flight test ampli- 
tudes to be within 15 percent of analytical values and 
in general closer to those computed using wind pene- 
tration effects. With a ten -samples -per -second rate, 
flight test transfer functions were limited to low fre- 
quencies. 


107 


IV. CONCLUSIONS 


The effect of wind penetration was seen to be 
significant at frequencies greater than 1.5 cps for the 
SA-9, 80-second flight condition. These wind pene- 
tration effects are aerodynamically introduced in the 
equations of motion by analyzing the local angle -of- 
attack distributions. Some insight as to the relative 
importance of these effects during a particular tra- 
jectory can be made by investigating the difference 
between an average angle of attack and one obtained 
for an instantaneous constant wind distribution over 
the entire vehicle length. This resolves into an 
analysis of the vehicle’s length with respect to the 
penetrating wind wave length. As the vehicle increases 
velocity , a given exciting wind frequency will require 
a larger wave length. This larger wave length will 
reduce the variation in local angle of attack. 

The relationship between angle -of -attack error 
profiles as a function of wind penetration frequency 
and vehicle flight time is shown on Figure 5 for the 
Saturn 201 vehicle. Although these percentile curves 
were obtained by an open loop analysis, they are 
reasonably accurate. These curves are to be inter- 
preted in the following manner: For example, neglect 
of wind penetration of a wind excitation at 1.7 cps 
occurring at 80 -second flight time will result in a 50 
percent error in angle of attack. 

Wind penetration effects can and should be in- 
cluded in the flight simulation analysis of a wind pro- 


file. Such an analysis would compute vehicle angle 
of attack, not at a single point, but at many stations 
along the vehicle and numerically integrate for the 
overall aerodynamic effect. 


Frequency (cps) 



FIGURE 5. ANGLE-OF -ATTACK ERROR 
PROFILES RESULTING FROM NEGLECTING 
WIND PENETRATION FOR SATURN 201 FLIGHT 
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VI 1 1. PUBLICATIONS AND PRESENTATIONS 


A. PUBLICATIONS 


ABSTRACT 


TECHNICAL MEMORANDUM X-53249 
April 23, 1965 

NEAR-OPTIMUM GUIDANCE - AN ANALYSIS 
OF FUEL PENALTY 

By 

Lyle R. Dickey 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


An explicit solution of the linearized equations of 
motion is used to obtain the deviations in end condi- 
tions , Ar and A0, and the additional burning time, 

At, expressed as integrals of functions of thrust angle 
deviations, Ax- Under the constraint that Ar = A0 = 

0, the calculus of variations is applied to minimize 
At. The resulting Euler -Lagrange equation evaluated 
at Ax = 0 gives a simple relationship which is used to 
show that At is a second order function of AX. This 
function is evaluated numerically for the second stage 
of an early SA-6 design. Results show that, if the 
mission is accomplished, the thrust angle may differ 
by as much as 2 degrees throughout the second stage 
with a propellant penalty of only 34 pounds. It follows 
that the capability of meeting the mission is the prime 
requisite of a guidance function and optimality is then 
only a second order consideration. 


TECHNICAL MEMORANDUM X-53290 
July 9, 1965 

F PS-16 RADAR/JIMSPHERE WIND DATA 
MEASURED AT THE EASTERN 
TEST RANGE 

By 

James R. Scoggins and Michael Susko 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


This report presents wind measurements at 
Eastern Test Range employing the FPS-16 Radar/ 

Jim sphere method. One hundred and twelve wind 
velocity profiles are presented. These measurements 
were made from November 1964 to May 1965. These 
data are not intended as a recommended statistical 
sample for structural and control system design crite- 
ria. They are presented to document this limited group 
of detailed wind profiles, and to provide a source for 
these data for use by other research organizations. 

The wind data presented herein supersedes the wind 
data presented in MTP-AERO-62-86. 


TECHNICAL MEMORANDUM X-53292 
July 12, 1965 

PROGRESS REPORT NO. 7 ON STUDIES IN THE 
FIELDS OF SPACE FLIGHT AND GUIDANCE 
THEORY 

Sponsored by Aero-Astrodynamics Laboratory of 
Marshall Space Flight Center 


ABSTRACT 


The progress reports of NASA -sponsored studies 
in the areas of space flight and guidance theory are 
presented. The studies are carried on by several 
universities and industrial companies. This progress 
report covers the period from July 23, 1964, to April 
1, 1965. The contracts are technically supervised 
by personnel of the Astrodynamics and Guidance 
Theory Division, Aero-Astrodynamics Laboratory, 
Marshall Space Flight Center. 


TECHNICAL MEMORANDUM X-53293 
July 13, 1965 

SA-10 FLIGHT MECHANICAL SUMMARY 

George C. Marshall Space Flight Center 
Huntsville, Alabama 
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ABSTRACT 


This report presents the Flight Mechanical Sum- 
mary for Saturn I vehicle SA-iO to be launched from 
Pad 37 B, Eastern Test Range,on a flight azimuth of 
95.2 degrees East-of-North. A successful flight will 
insert the spent S-IV stage and payload consisting of 
an Apollo boilerplate (BP -9) and a Meteoroid Tech- 
nology Satellite ( Pegasus C) into a 535-km circular 
orbit. Included is a discussion of the operational 
predicted trajectory with its orbital mission objec- 
tives and constraints, guidance, sequencing, and in- 
sertion parameters. The three -sigma envelope of 
thrust, flow rate, liftoff mass, and wind speed vari- 
ations is given, along with impact data, in the Range 
Safety part of this report. A study of wind disturb- 
ances, three-sigma C t and C 2 variations, and ± 10 
percent control gain deviations shows the vehicle to 
be structurally capable of withstanding the expected 
winds during the launch period. Pertinent tracking 
and telemetry ground coverage data are presented 
for powered and orbital flight phases. 

TECHNICAL MEMORANDUM X-53298 
July 21, 1965 

A PRACTICAL APPROACH TO THE OPTIMIZATION 
OF THE SATURN V SPACE VEHICLE CONTROL 
SYSTEM UNDER AERODYNAMIC LOADS 

By 

Robert S. Ryan and Dieter Teuber 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


The equations for the bending moment of a launch 
vehicle are written including the effects of bending 
and sloshing dynamics. By making certain simplify- 
ing assumptions, the bending moment equation is 
transformed to be a function of rigid vehicle angle of 
attack, engine deflection, bending mode acceleration 
and sloshing mode acceleration. This form is partic- 
larly advantageous for use in optimization techniques 
by the control engineer. The response of a Saturn V 
class vehicle was solved and the control system opti- 
mized to bending moments using data covering a 
period of six years, amounting to 4, 384 measured 
wind profiles, and the GPS high speed repetitive an- 


alog computer. A comparison between the bending 
moment response envelope of the measured winds and 
the bending moment response of the MSFC synthetic 
wind profile is made. 


TECHNICAL MEMORANDUM X-53303 
July 26, 1965 

HOPI-DART AND CAJUN-DART ROCKET 
WIND MEASURING SYSTEMS 

By 

* 

Robert E. Turner and Luke P. Gilchrist 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


Three Hopi-Dart rockets were fired in February 
1964, and during February 1965 through June 1965 
seventeen Cajun-Dart rockets were fired to obtain 
wind flow data. These wind flow data are obtained by 
radar tracking of Mylar chaff expelled from these 
rocket systems. An altitude profile of the measured 
winds and a description of the two -stage rocket - 
delivery systems are given. 

In February 1964 the observed winds were 
westerly. In 1965 the observed winds were westerly 
during February and March, changing to easterly in 
late March and continuing through May, and reversing 
to northwesterly in June. These data are presented 
graphically. 

TECHNICAL MEMORANDUM X-53319 
August 16, 1965 

CONIC SOLUTIONS TO THE INTERPLANETARY 
TRANSFER PROBLEM WITH CONSTANT OR 
MINIMUM INJECTION ENERGY 

By 

Lamar E. Bullock 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


* 

Luke P. Gilchrist is associated with the Lockheed 
Corporation working under Contract No. NAS 8-20082. 
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ABSTRACT 


This report describes a method for calculating 
minimum and constant injection energy, C 3 , values 
for conic interplanetary trajectories. The planets are 
assumed to move on mutually inclined elliptical or- 
bits about the sun. Lambert’s theorem is used to 
solve the time-constraint on the two -point boundary 
value problem, and an iterative procedure, with ar- 
rival date as the isolation parameter, is used to de- 
termine the desired C 3 . 

Typical results are presented for Earth-to-Jupi- 
ter trajectories in the 1970-71 launch period. The 
comparisons made with previously used methods in- 
dicate the procedures described in the report produce 
results of comparable accuracy, with a savings in 
computation time. 

TECHNICAL MEMORANDUM X-53325 
September 2, 1965 

STATISTICAL ANALYSIS OF PHOTOGRAPHIC 
METEOR DATA - PART 1 - OPIK’S LUMINOUS 
EFFICIENCY AND SUPPLEMENTED WHIPPLE 
WEIGHTING 

By 

Charles C. Dalton 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


After adjustment of the estimate for the mass of 
a zero-absolute-visual -magnitude meteor, formulas 
for meteoroid mass are derived from E. J. Opik’s 
(1958) physical theory for dustball meteors. Masses 
are calculated for G. S. Hawkins and R. B. South - 
worth’s (1958) random sample of 285 sporadic photo- 
graphic meteors. Weighting functions are derived 
for the statistical minimization of the effects of sev- 
eral sources of bias. Means, standard deviations, 
and ordinary, multiple, and partial correlations are 
shown for a 5 -variate statistical analysis. Weighted 
cumulative distributions for parameters are plotted 
for two gradations with respect to mass. Data -point 
distributions are plotted which establish cumulative 
flux as functions of mass, momentum, kinetic energy, 


etc. , which show the effects of the earth’s motion on 
the distributions and correlation of parameters. 

The results support revisions in the author’s 
(Jan. 1965) interpretation of Explorer XVI data and 
predicted puncture rates for the Pegasus meteoroid 
measurement satellites. The apparent discrepancy 
between the logarithmically linear relations between 
cumulative flux and mass, as inferred from photo- 
graphic meteors by G. S. Hawkins and E. K. L. 
Upton ( 1958) and from zodical light and solar F- 
corona data by D. B. Beard ( 1963), is shown to be 
accountable by E. J. Opik’s (1951) probability that 
a particular meteoroid will encounter the earth in 
one revolution of the particle. 


TECHNICAL MEMORANDUM X-53326 
September 3, 1965 

SUMMARY OF BASE THERMAL ENVIRONMENT 
MEASUREMENTS ON THE SATURN I BLOCK I 
FLIGHT VEHICLES 

By 

Ira P. Jones, Jr. 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


Base thermal environment measurements were 
made on all of the Saturn I Block I vehicles (SA-i 
through SA-4) . Total and radiation calorimeter 
measurements were made in the base region along 
with base gas temperatures and pressures. Details 
of the instrumentation and problems encountered in 
the data evaluation are discussed, and the data are 
compared from flight to flight on the basis of altitude. 
Some representative correlations of total and radia- 
tion heating, gas temperatures, and base pressures 
with small scale model tests are included. Flight 
heat transfer coefficients are calculated and com- 
pared with the model results. 

TECHNICAL MEMORANDUM X-53327 
September 7, 1965 

DERIVATION OF PARAMETRIC VALUES FOR A 
GREENHOUSE MODEL OF THE CYTHEREAN 
ATMOSPHERE 
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By 

Robert B. Owen 

George C. Marshall Space Flight Center 
Huntsville, Alabama 

ABSTRACT 


Presented is a self-consistent model atmosphere 
for the planet Venus, with emphasis on parameters of 
general use in developing vehicle design criteria. 

Kern and Schilling's program [1] for generation of 
atmospheric parameter values is coupled with the 
most recent estimates of surface conditions in order 
to derive the atmospheric environment at different 
altitudes. A range of initial values is used to produce 
a parametric band most representative of actual con- 
ditions. Because of the extreme pressures and tem- 
peratures involved, it is felt that design for such an 
environment presents an engineering problem of the 
first magnitude. 

TECHNICAL MEMORANDUM X-53328 
September 9, 1965 

TERRESTRIAL ENVIRONMENT (CLIMATIC) 
CRITERIA GUIDELINES FOR USE IN SPACE 
VEHICLE DEVELOPMENT, 1965 REVISION 

By 

Glenn E. Daniels 
James R. Scoggins 
and 

Orvel E. Smith 

George C. Marshall Space Flight Center 
Huntsville, Alabama 

ABSTRACT 

This document provides guidelines on probable 
climatic extremes and probabilities -of -occurrence of 
terrestrial environment data specifically applicable 
for NASA space vehicles and associated equipment 
development. The geographic areas encompassed are 
the Eastern Test Range (Cape Kennedy, Florida); 
Huntsville, Alabama; New Orleans, Louisiana; the 
Western Test Range (Point Arguello, California); 
Sacramento, California; Wallops Test Range (Wallops 
Island, Virginia); White Sands Missile Range, New 
Mexico; and intermediate transportation areas. In 


addition, a section has been included to provide infor- 
mation on the general distribution of natural environ- 
ment extremes in the United States (excluding Alaska 
and Hawaii) that may be useful to specify design cri- 
teria in transportation of space vehicle components 
from subcontractors. This document omits climatic 
extremes for worldwide operational conditions. This 
is consistent with the existing philosophy regarding 
the employment of large space vehicles, since launch- 
ing and test areas are relatively restricted due to the 
availability of facilities and real estate. 

This document presents the latest available 
information on probable climatic extremes, and 
superseded information presented in TM X-53023. 

The information in this document is recommended for 
employment in the development of space vehicles and 
associated equipment, unless otherwise stated in con- 
tract work specifications. 

TECHNICAL MEMORANDUM X-53344 
October 7, 1965 

GROUND LEVEL ACOUSTICAL FOCI IN A 
THREE-LAYERED ATMOSPHERE 

By 

Willi H. Heybey 

George C. Marshall Space Flight Center 
HuntsviUe, Alabama 

ABSTRACT 


An infinitesimally slender bundle of sound rays 
emitted by a point source can, under certain meteo- 
rological conditions, be returned to the horizontal 
plane on which the source sits, and at the same time, 
may be narrowed down to enclose an area actuaUy 
equal to zero. In such cases the common point of 
arrival is called an acoustical focus. The intensity 
of returned energy there attains a very high degree 
capable of damaging effects. 

For determining such points the atmosphere is 
usually divided up into layers in each of which the 
change of the sound speed with height is considered 
constant. Three layers often suffice to approximate 
the many layers arising from observation. 

The report studies focus formation engendered 
by the third layer in coUaboration with the two lower 
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layers. A focus may appear in variegated circum- 
stances. For each case a computational scheme is 
given which may be followed by slide rule or assigned 
to machines if higher accuracy or a systematic survey 
is desired. 


TECHNICAL MEMORANDUM X-53351 
October 21, 1965 

AN ENVIRONMENTAL MODEL FOR VAN ALLEN 
BELT PROTONS 

By 

William T. Roberts 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


The majority of the manned and unmanned 
missions in space for at least the next ten years will 
be composed mainly of earth orbiting satellites and 
space stations. Although there are many environ- 
mental aspects of space which must be explored and 
elaborated upon before the successful launching of 
such missions, few are so critical as the radiation 
environment which will be encountered. 

The radiation environment for earth orbiting 
vehicles will consist primarily of ( 1) Van Allen 
particles, (2) solar flare particles, and (3) galactic 
cosmic particles. Of these three, the Van Allen 
belt particles present the greatest hazard. 

This report summarizes the knowledge of the 
energy spectrum and distribution of the Van Allen 
proton belt as it is known today. A computer pro- 
gram has been developed for calculating Van Allen 
particles incident upon vehicles in earth orbits. 

Based upon the 1963 FLUX program of McDwain, we 
find that a vehicle in a highly elliptic orbit may 
encounter as many as 1. 1 x 10 8 protons/cm 2 / orbit 
with 40 Mev ^ F ^ 110 Mev. In the magnetic equa- 
torial plane, the flux of protons with E > . 5 Mev may 
be as high as 2. 2 x 10 7 protons/c m 2 /sec. 


TECHNICAL MEMORANDUM X-53354 
November 1, 1965 


LAUNCH WINDOWS FOR TWO TYPES OF ORBITS 
SYNCHRONOUS WITH THE LUNAR PERIOD 

By 

E. H. Bauer and L. D. Mullins 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 

This report presents launch windows for two 
types of orbits designed to be synchronous with the 
lunar orbital period. One class of conic, the 1/8- 
class orbit, maintains a period one-eighth that of the 
lunar period and an apogee of approximately one -half 
the lunar distance. The second type, 1/2 -class or- 
bit, consists of a conic with one-half the lunar period 
and apogee radius approximately equal to the lunar 
distance, an Arenstorf Orbit. Launch windows are 
presented which reflect the assumptions and con- 
straints applicable to the specific geometry of each 
class and to the S-IB/S-IVB/Centaur vehicle. The 
time period considered for the 1/2-class orbit is 
December 1967 through June 1968 where only discrete 
launch opportunities are found. The window for the 
1/8 -class is less constrained and is shown for each 
day of any year. 

TECHNICAL MEMORANDUM X-53360 
November 18. 1965 

STATISTICAL ANALYSIS OF PHOTOGRAPHIC 
METEOR DATA, PART E: VERNIANTS 
LUMINOUS EFFICIENCY AND 
SUPPLEMENTED WHIPPLE 
WEIGHTING 

By 

Charles C. Dalton 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


This report is Part E in a series of four reports 
on the statistical analysis of the Hawkins and South- 
worth random sample of 285 sporadic photographic 
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meteors. The four parts of the analysis comprise 
the combinations from two alternative formulations 
for meteoroid mass and two alternative formulations 
for data weighting. Parts I and II have the same 
weighting as a function of air -entry velocity and other 
parameters. But the surprisingly large disparity 
between the results with Opik’s luminous efficiency 
suggests the possibility of considerable bias from 
weighting inversely with the 1. 5 -power of velocity. 


TECHNICAL MEMORANDUM X-53363 
December 3, 1965 

MEASUREMENTS OF WINDS BY CHEMICAL 
RELEASES IN THE UPPER ATMOSPHERE 

By 

Hugh Wilson Morgan 

George C. MarshaU Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


Upper atmospheric winds are of importance to 
the designers of space vehicles as well as to geophys- 
icists. For this reason a number of organizations 
have begun to study winds above 80 kilometers by re- 
leasing into the upper atmosphere chemicals which 
luminesce either by the action of sunlight or by a re- 
action with naturally occurring molecules. Seventy - 
seven of these releases were studied in an effort to 
determine a statistical model for wind velocities 
between 80 and 200 kilometers at mid -latitude sites. 


TECHNICAL MEMORANDUM X-53365 
December 3, 1965 

BOUNDARY VALUE PROBLEMS ASSOCIATED 
WITH OPTIMIZATION THEORY 

By 

Hugo L. Ingram 

George C. Marshall Space Flight Center 
HuntsviUe, Alabama 


ABSTRACT 


Optimization theory is applied to the physics of 
trajectory problems to yield a boundary value formu- 
lation. Three practical numerical methods for ob- 
taining solutions to boundary value problems are dis- 
cussed, and a detailed explanation of the application 
of the most efficient method to a sample problem is 
included. Numerical methods are needed because the 
application of optimization theory to trajectory prob- 
lems generally results in nonlinear differential equa- 
tions connecting the boundary conditions that are not 
all specified at the same value of the independent 
variable. The methods discussed are adaptable to the 
solution of boundary value problems other than the 
ones associated with trajectory problems or optimi- 
zation theory. 


TECHNICAL MEMORANDUM X-53367 
December 10, 1965 

ELASTIC STABILITY OF A SLENDER BAR WITH 
FREE -FREE ENDS UNDER DYNAMIC LOADS 

By 

Frank C. Liu 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


Problems of elastic stability of a space rocket under 
dynamic, nonconservative loads are presented. The 
spacecraft is treated as a slender, uniform cylindrical 
bar with free-free ends. Analysis of the frequency of 
the transverse vibration coupled with the longitudinal 
vibration under the action of time -varying thrust force 
based on Galerkin's method leads to Mathieu ! s equa- 
tion. The elastic stability boundaries of the frequency 
ratio, longitudinal to transverse, versus the thrust are 
obtained. The dynamic buckling thrust of a uniform 
bar with various boundary conditions is presented. 

The aerodynamic force induced from the oscillation 
of the rocket in supersonic speed on the critical 
thrust is investigated. A numerical illustration is 
given to show the relationship of the critical thrust 
versus flight speed at various flight altitudes. 


TECHNICAL MEMORANDUM X-53368 
December 13, 1965 
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NUMERICAL PROCEDURES FOR ROLL 
STABILITY STUDIES 


ABSTRACT 


By 

Philip J. Hays 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


Practical simplified procedures are developed in 
this paper for calculating estimates of parameters of 
the negative binomial distribution with probability 
function 


f(x) = 


r (x + k) 

x: r<k) 


k /, x X f\ 

p ( i - p) ; X = 0, 


i, 2. . . 


ABSTRACT 


This report illustrates the numerical procedure 
used by the Aero-Astrodynamics Laboratory in per- 
forming roll stability analyses for Saturn type space 
vehicles. Vehicle control in the roll channel is pro- 
vided by a position gyro and a rate gyro. The dynamics 
of the rate gyro are included as a separate equation in the 
equations of motion. The modes of oscillation in this 
system are (1) torsion, (2) roll, (3) roll sloshing, 

(4) rate gyro roll, and (5) swivel engine. Two 
methods are available for solving the system’s 
eigenvalues: the matrix iteration method and the 
characteristic equation approach. In both cases the 
system is solved for the normalized eigenvectors, 
which provide additional information about the system. 

The frequency response of the system, at the 
control thrust vector point, can be found for a sinu- 
soidal forcing function at various sensor locations. 

The equations of motion are derived in the 
appendix. 


where 0 < p < 1 and k > 0. Moment estimators, max- 
imum likelihood estimators, and estimators based 
on moments and frequencies in selected classes are 
given both for the complete and for the truncated 
(with missing zero class) distribution. To facilitate 
calculation of the various estimators given, a table of 
the function -p In p/( 1 - p) with entries to six deci- 
mals at intervals of 0.001 for the argument p, is in- 
cluded. Illustrative examples are also included. 

TECHNICAL MEMORANDUM X-53362 
Nov. 30, 1965 

LINEAR FEEDBACK GUIDANCE 
By 

Lyle R. Dickey 

George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


TECHNICAL MEMORANDUM X-53372 
December 21, 1965 

ESTIMATION IN THE NEGATIVE BINOMIAL 
DISTRIBUTION 

By 

A. Clifford Cohen, Jr. 

George C. Marshall Space Flight Center 
Huntsville, Alabama 

Professor of Mathematics, University of Georgia, 
Athens, Georgia. The research reported in this 
paper was performed under NASA Contract 
NAS8-11175 with the Aerospace Environment Office, 
Aero-Astrodynamics Laboratory, Marshall Space 
Flight Center, Huntsville, Alabama. Mr. O. E. 
Smith is the NASA contract monitor. 


The Problem of determining the coefficients for 
a linear feedback guidance system with polynomial is 
investigated. The solution to the linearized equations 
of motion is employed to determine the first partial 
derivatives of the optimum thrust angle deviations with 
respect to position and velocity coordinates. The 
partial derivatives are approximated by quadratic time 
functions, which in turn are the first approximation 
to a time-variable set of gains. This system is then 
analyzed, and two of the weighting functions are 
linearized with respect to perturbations in the poly- 
nomial coefficients. Changes in these coefficients 
are determined to minimize the weighted sum of 
squares of these expressions over several time points. 
The resulting changes were made, and a satisfactory 
reduction in magnitude of these weighting functions 
resulted. This, coupled with a linear feedback of 
thrust acceleration and a small time function derived 
to cancel the effect of initial conditions, produced 
the final guidance function which performed exception- 
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ally well on a 100 n. m. orbital mission of an early 
SA-6 second stage vehicle. Of particular interest is 
a demonstration of the effectiveness of the numerical 
methods used which are generally applicable to a wide 
variety of guidance and control feedback problems. 


NASA TN D-3117 
November, 1965 

MONTE CARLO APPROACH TO TOUCHDOWN 
DYNAMICS FOR SOFT LUNAR LANDING 

By 

Robert E. Lavender 
ABSTRACT 


Results of analytical touchdown dynamics inves- 
tigations are presented which were conducted to 
obtain estimates of the probability of stable landing 
for configurations with various landing gear diam- 
eters. Both three-legged and four-legged vehicles 
are considered in the analysis. Dynamic scaling 
considerations are taken into account so that the re- 
sults are applicable to a wide range of vehicle size 
and mass. A Monte Carlo approach is taken in the 
determination of initial landing conditions. 

Results indicate that, for a given probability of 
stable landing, a three-legged vehicle requires a 
landing gear diameter only slightly larger than the 
diameter required for a four-legged vehicle. There- 
fore, the three-legged vehicle's landing gear should 
weigh less. 


B. PRESENTATIONS 


AN ANALYTICAL REDUCTION OF THE OPTIMAL 
TRAJECTORY PROBLEM 

By 

Rowland E. Burns 


ABSTRACT 


The problem of optimal thrust programming of 
a single stage rocket powered vehicle to achieve 
maximum payoff at specified end conditions is treated. 


The vehicle is assumed to operate in vacuo in an 
inverse square gravitational field. For the case of 
the so-called singular arc, the Lagrange multipliers 
are determined, and the mass flow rate to maintain 
a singular arc condition is derived. 


M. S. Thesis, Mathematics, August 1965 
University of Alabama 


STUDIES OF SATURN BOOSTER RECOVERY 
By 

Jesco von Puttkamer 

Aero-Astrodynamics Laboratory 
NASA -George C. Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


The paper deals with the discussion of booster 
recovery studies, as performed at MSFC, and their 
results. 

First, the rationale of booster recovery studies, 
the historic background of the subject, and the basic 
principles of booster recovery are presented in a 
cursory maimer. Then, the most recent thinking in 
the area of booster recovery is illustrated by sum- 
marizing the results of recent studies of SATURN-V 
first stage (S-IC) recovery at Marshall. The re- 
quirements and design criteria of suitable recovery 
systems are outlined by treating the subject step by 
step according to actual events encountered during the 
recovery sequence. 

By necessity, many of the points discussed are 
engineering considerations, directly related to "rou- 
tine” booster stage hardware. 


Presented at the University of Minnesota's course 
"Final Phase, AERODYNAMIC DECELERATION," 
University of Minnesota, July 6-16, 1965, 
Minneapolis, Minnesota. 

ON TOUCHDOWN DYNAMICS ANALYSIS FOR 
LUNAR LANDING 

By 

Robert E. Lavender 
Scientific Assistant to Director 
Aero-Astrodynamics Laboratory 
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George C. Marshall Space Flight Center 
Huntsville, Alabama 

ABSTRACT 


Analytical methods which may be used for ob- 
taining touchdown dynamics results vary from a 
simple rigid body -impulsive force procedure which 
accounts for only a few parameters to elaborate 
three-dimensional numerical procedures which at- 
tempt to account for practically all of the parameters 
affecting the motion. Landing stability boundaries 
obtained from six methods are compared to each 
other as well as with results from experimental drop 
tests. The results show considerable variation be- 
tween the various methods. A review of each method 
(parameters are considered in each) and comparison 
of the analytical results lead to the conclusion that the 
most detailed procedure considered gives results 
which most closely correspond to the results which 
would be obtained by a full scale vehicle landing on 
the moon. This conclusion is reached even though 
another method gives the best correlation with scale 
model experimental drop tests. It is further con- 
cluded that dynamically scaled models must be con- 
structed with great care in order to obtain the proper 
representation of the full scale vehicle. Results are 
given for both two-dimensional symmetric landings 
as well as three-dimensional landings. 


AIAA Symposium on Structural Dynamics and 
Aeroelasticity 

MIT, Cambridge, Massachusetts 
August 30-September 1, 1965 


LIQUID-FREE SURFACE INSTABILITY RESULTING 
FROM RANDOM VERTICAL ACCELERATION 

By 

L. L. Fontenot 
General Dynamic s/Convair 
Huntsville, Alabama 

G. F. McDonough 

Aero-Astrodynamics Laboratory, MSFC 
D. O. Lomen 

General Dynamics/Convair 
San Diego , California 


ABSTRACT 


In this paper the stability of the plane -free sur- 
face of a heavy liquid enclosed in a vertical cylinder .. 
undergoing random vertical accelerations is investi- 
gated. The cylinder studied has an arbitrary cross 
section and a flat bottom and stability is investigated 
in a mean square sense. The liquid is assumed fric- 
tionless, homogeneous and incompressible, and vis- 
cous and capillary forces between liquid and container 
are considered negligible. Fundamental hydrodynamic 
equations are used as the basis of the mathematical 
description of the fluid behavior. 

Random inputs in the form of white noise are 
applied and the correlation functions between princi- 
pal modes of oscillation are determined (considering 
the ensemble average) and, from the solution of the 
resulting integral equations stability of the surface is 
established by use of the Routh-Hurwitz criterion. 
This shows that the system is always unstable for 
white noise inputs if damping is not present or is 
sufficiently small. 


Presented at the Sixth International Symposium on 
Space Technology and Science, Tokyo, November 
29 -December 4, 1965. 


DYNAMIC STABILITY OF THIN ELASTIC PLATES 
UNDER THE ACTION OF NON-DETERMINISTIC 
LOADS 

By 

D. O. Lomen 

General Dynamics/Convair 
San Diego , California 

L. L. Fontenot 
General Dynamics/Convair 
Huntsville, Alabama 

G. F. McDonough 
Aero-Astrodynamics Lab, MSFC 
Huntsville, Alabama 

ABSTRACT 


This paper treats the problem of a simply 
supported rectangular plate under the action of ran- 
dom compressive forces applied at the edges. The 
theory of linear stochastic systems is applied to the 
problem for an input of pure white noise. This theory 
predicts the occurrence of lateral instability of the 
plate in a mean square sense for any amount of white 
noise in the case of very low damping. 
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Preliminary numerical results have been ob- 
^ tained by analog simulation for the case of zero 
damping (other than that inherent in the elements of 
the computer) , A first attempt at correlating the 
average time rate of change of amplitude with a fre- 
quency parameter , there appears to be a definite 
straight-line trend in the data plot (log -log scale) 
although considerable scatter exists in the data. 


Presented at the Sixth International Symposium on 
Space Technology and Science, Tokyo, November 29- 
Dec ember 4, 1961. 


STABILITY PROBLEMS IN THE CONTROL OF 
SATURN LAUNCH VEHICLES 

By 

G. F. McDonough 
Aero-Astrodynamics Laboratory 
Marshall Space Flight Center 
Huntsville, Alabama 


ABSTRACT 


The control system for attitude control of a 
launch vehicle during powered flight is designed to 
insure that the engine thrust vector is so oriented 
that attitude commands do not result in vehicle ro- 
tation instabilities. Stability problems arise because 
the vehicle is elastic rather than rigid, contains pro- 
pellants which are free to oscillate, and is at most 
times aerodynamically unstable. 

Elastic body bending modes arise from engine 
gimbal forces and aerodynamic forces. Control sen- 
sors must differentiate such modes from rigid body 
modes to prevent the occurrence of closed -loop os- 
cillations. Methods of stabilization include phase 
stabilization (achieved by sensor location, shaping 
control networks, adaptive control techniques, or 
sensor signal blending), gain stabilization (by at- 
tenuation) , or by change in configuration to eliminate 
instability of aerodynamic forces. 

Propellant oscillations can couple through the 
control system to produce closed-loop oscillations. 
This instability can be reduced or eliminated through 
choice of tank location, shaping of control system 
networks, or by mechanical means such as baffles. 

Engine compliance modes can couple with the 
vehicle modes; these can be eliminated through ac- 
tuator and thrust structure design. 


Control mode effects can be eliminated by pro- 
per choice of control frequency or by providing high 
damping in the control mode. 

The difficulty in handling these problems for 
the Saturn V vehicle lies in the close grouping of fre- 
quencies for the various modes ; stabilization of one ’ 
mode by present techniques generally worsens the 
effect of one or more other modes. For this reason, 
solutions of vehicle stability problems are compro- 
mises. 

Several criteria are available for examining 
vehicle stability; however, these require linear equa- 
tions with constant coefficients. To provide these 
coefficients time -slice investigations are used. Ve- 
hicle characteristics must be determined; this re- 
quires extensive testing as well as theoretical analysis 
and analog investigations. Where unknowns cannot 
be eliminated, latitude must be provided in the control 
system design. 


Presented at the 1T Internation Conference on Dynamic 
Stability of Structures" at Northwestern University, 
October 1965. Sponsored by Northwestern University 
and the Air Force Office of Scientific Research. 


DYNAMIC TESTING OF SATURN LAUNCH 
VEHICLES 

Bp 

George F. McDonough 


ABSTRACT 


Dynamic testing of full-scale launch vehicles 
is an important part of the Saturn program at Marshall 
Space Flight Center. Results of these tests are used 
to verify existing cfynamic analyses as well as to pro- 
vide a basis for the development of new theories 
where these are needed. In addition, necessary in- 
formation is obtained from these tests which either 
cannot be obtained elsewhere or, if there are other 
sources, these entail additional expense or time to 
achieve the same result. 

The Saturn V vehicle will serve as a suitable ex- 
ample of the application of cfynamic testing. When 
this vehicle was first conceived, there was no exist- 
ing knowledge of the dynamic behavior of clustered 
vehicles. Dynamic tests of scale models were con- 
ducted and theories of the behavior were developed; 
however, neither of these approaches was to prove 
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adequate in defining the required vehicle character- 
istics either for structural design or for control 
system design to insure stability during flight. Pri- 
marily through the application of the results of 
carefully controlled full-scale dynamic testing at 
MSFC, together with necessarily limited flight test 
results on Saturn I, the dynamic characteristics of 
clustered vehicles are reasonably well understood 
today. As a consequence, the amount of dynamic 
testing that has been performed on the second- 
generation Saturn IB is extremely limited. There 
remain certain vehicle characteristics which cannot 
as yet be determined theoretically with sufficient 
accuracy; chief among these are structural damping 
and local damping effects, which are particularly 
important in the definition of sensor response for the 
flight control system. In addition, much of the 
necessary knowledge on the complex elastic connec- 
tions of the clustered vehicles must be obtained 
through testing, and can best be gained through dy- 
namic testing. 

In the Saturn V program , once again we are 
faced with extremely important unknowns. As in the 
past, the ground rules under which we operate em- 
phasize that we do not have the luxury of development 
flights; in short, there can be no flight failures if the 
program is to meet its schedules. The only means 
by which we can reasonably insure the required re- 


liability is through a well conceived and carefully 
executed ground test program which includes dynamic * 
testing of a flight configuration vehicle. Among the 
more urgent problems for the Saturn V vehicle are 
the longitudinal vibration problem and local effects , 
which become even more important as the vehicle be- 
comes more flexible and as nonlinear behavior exerts 
a greater influence. Also, damping and three- 
dimensional effects are increasingly important. A 
well-balanced program of theoretical effort to define 
these effects is underway at MSFC, as well as 
throughout NASA and the aerospace industry; however, 
these efforts will not provide the critical information 
without verification by dynamic testing. Furthermore, 
the state-of-the-art of dynamic testing has had to be 
improved before the required results could be ob- 
tained. Several improvements have been made, a 
few examples of which are (i) an improved support 
system to eliminate damping in the support system to 
allow better measurement of structural damping, 

(2) use of scale model tests to aid in full-scale test 
planning and to enable earlier checks of theory to be 
made, and (3) improved instrumentation systems and 
data handling systems to allow a maximum amount of 
accurate data to be obtained in the shortest possible 
time to offset severely restricted test schedules. 
Presented at the Sixth International Symposium on 
Space Technology and Science, Tokyo, November 29- 
December 4, 1965. 
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